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1.0  PURPOSE 

The  work  undertaken  under  this  contract  and  presented  in  this  Report 
is  related  to  research  investigations  to  contribute  to  the  development 
of  air-breathing  fuel  cells  based  on  ion  exchange  membranes  as 
"solid"  electrolytes. 

Technical  guidelines  EPP  No.  58238  of  the  contract  were  defined  as 
follows: 

1.  These  Technical  Guidelines  cover  the  requirements  for  the 
investigation  of  hydrogen-oxygen  ion  exchange  membrane 
fuel  cells.  These  guidelines  further  cover  limited  investi¬ 
gations  of  fuels  other  than  hydrogen  which  show  practical 
and  economical  feasibility  with  respect  to  military  applica¬ 
tions. 

2.  The  objective  of  these  investigations  is  to  improve  the  power 
density  and  lifetime  characteristics  of  ion  exchange  membrane 
fuel  cell  power  supplies  for  ground  power  applications.  Cells 
utilizing  fuels  other  than  hydrogen  shall  also  be  investigated. 
Emphasis  shall  be  placed  on  the  improvement  of  ion  exchange 
membrane  fuel  cell  technology  and  a  systematic  understanding 

-  of  the  phenomena  observed.  This  effort  should  lead  to  an 

increase  in  the  state  of  the  art  which  will  allow  the  design  of 
improved  fuel  cell  power  sources  to  meet  groundpower  require¬ 
ments  with  a  known  reliability. 

3.  The  Contractor  shall  conduct  the  necessary  investigations 
leading  to  the  establishment  of  an  analytical  model  expressing 
factors  determining  fuel  cell  performances.  This  model  shall 
be  used  to  further  improve  the  overall  operating  performance 
of  the  ion  exchange  membrane  fuel  cell.  The  Contractor  shall 
perform  the  following  experimental  and  theoretical  studies: 

a.  Determine  the  factors  causing  ion  exchange  membrane 
oxygen  electrode  performance  to  decrease  when  oper¬ 
ating  on  convected  ambient  air  as  the  source  of  oxygen. 
Work  will  be  directed  toward  obtaining  optimum  air 
electrode  performance. 

b.  Conduct  fuel  cell  performance  and  life  studies,  using 
new  improved  cation  exchange  membrane  materials. 
Emphasis  shall  be  placed  on  structures  which  improve 
ionic  conductance,  physical  properties,  and  stability. 
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PURPOSE  (Cont'd) 

c.  Complete  a  study  of  factors  determining  the  membrane's 
water  balance. 

d.  Investigate  factors  influencing  electrode-electrolyte 
coupling  to  produce  low  electrode  polarization. 

This  work  shall  be  directed  toward  optimization  of 
ion  exchange  membrane  electrodes. 

e.  Investigate  fuels  other  than  hydrogen  in  ion  exchange 
membrane  fuel  cells.  These  fuels  shall  be  restricted 
to  those  found  by  other  investigators  to  work  such  as 
dissolved  liquid  fuels  in  acid  electrolytes.  Emphasis 
shall  be  given  to  determining  the  Faradaic  efficiency, 
main  reactions,  and  causes  of  anode  polarizations. 
Studies  of  the  methanol  ion  exchange  membrane  cell 
shall  be  included  in  these  investigations. 

The  initial  experime  ntal  program  started  in  October  1962,  and 
approved  on  November  1,  1962,  is  considered  compatible  with 
the  Technical  Guidelines. 


The  following  TASK  organization  was  adopted: 


TASK  I:  Interface  Structure  and  Interfacial  Properties  of 

IEM/Electrode  Systems 

1.  Isotropic  Characteristics  of  an  Ion  Exchange 
Membrane. 

2.  Ion  Exchange  Membrane/Electrode  Morphology. 

3.  Interfacial  Bonding  Forces. 


TASK  II:  Water  Transport  Across  Ion  Exchange  Membranes 

1.  Water  Permeability. 

2.  Electro-osmotic  Water  Transport. 

3.  Water  Gradients  in  Membranes. 

TASK  III:  Kinetic  Study  of  Electrochemical  Oxygen  Reduction  at 

Variable  Oxygen  Partial  and  Total  Pressures 

1.  Electrode  Reproducibility. 

Z.  Variables  Affecting  Reaction  Rates. 

3.  Transient  Rate  Measurements. 

4.  Steady-State  Rate  Measurements. 
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PURPOSE  (Cont'd) 


TASK  IV:  Heat  and  Mass  Transfer  Investigations  of  the  Air-Breathing 

I  EM  Cell  and  Battery 


1.  Over-all  Mass  Transfer  Processes. 

2.  Detailed  Mechanism. 

3.  Temperature  Distribution  and  Profile. 

4.  Property  Measurements. 

5.  Analytical  Model. 

6.  Life  Testing. 


TASK  V: 


Design  Parameters  Affecting 


oerated  IEM  Fuel  Cells 


1.  Performance  Data  as  Related  to  Water  Transport. 

2.  Systems  Evaluation  for  Product  Water  Conservation. 

3.  Air  Electrode  for  Constant  Water  Removal. 

4.  Air  Electrode  Configuration  and  Geometry. 

5.  Cell-Stack  Operation. 


TASK  VI: 


Ion  Exchange  Membranes 


1.  New  Ion  Exchange  Membrane  Materials. 

2.  Materials  Properties. 


Postponed  TASK  VII:  IEM  Methanol  Fuel  Cell 

1.  Methanol  Transport  through  IEM. 

2.  Rate  of  Methanol  Oxidation  in  Solution. 

3.  Reduction  of  Cathodic  Electrode-Poisoning. 

4.  Reaction  Rates  of  Intermediate  Products. 

This  Report  presents  experimental  work  regarding  TASKS  I,  II,  HI, 
IV  and  V. 
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ABSTRACT 


During  the  contract  period  from  January  1,  1963  to  June  30,  1963, 
experimental  work  has  been  conducted  on  membranes  and  inter¬ 
face  membrane/electrolyte  properties.  Reaction  rates  for  the 
electrochemical  oxygen  reduction  as  well  as  heat  and  mass  trans¬ 
fer  studies  have  yielded  interesting  results  regarding  air-breathing 
fuel  cell  electrodes.  The  reported  results  suggest  a  continued 
effort  of  the  proposed  program  throughout  the  remaining  contract 
period. 

1.  The  study  of  the  IEM-electrode  interface  by  optical 
microscopy,  capacity  and  adhesion  measurements  has 
yielded  results  regarding  type,  quality  and  nature  of 
the  electrolyte-electrode  bond.  Surface  contact  between 
the  two  media  rather  than  interfacial  penetration  has  been 
observed.  The  presence  of  a  fluorocarbon-type  wet¬ 
proofing  agent  has  shown  to  be  effective  as  a  water- 
rejecting  medium  and  generates  stronger  binding, 
resulting  practically  in  better  resistance  to  thermal  or 
mechanical  degradation  by  cycling  processes,  between 
membrane  and  electrode. 


2.  Hydraulic  transport  rates  in  IEM,  important  in  the  overall 
water  transport  processes  (heat  and  mass  transfer,  thermal 
convective  processes,  rates  of  evaporation  in  presence  and 
absence  of  forced  gas  flow,  etc. . .  )  have  yielded  values  for 
permeability  ranging  from  1.  3  to  25  •  10"  cm^-cm/ 
cm^-sec,  atm.  at  room  temperature.  These  transport 
rates  are  largely  dependent  on  ion  exchange  capacity,  and 
corrections  for  tortuosity  factors  have  allowed  the  deriva¬ 
tion  of  a  generalized  equation  for  all  IEC's.  Transport 
rates  are  temperature-dependent  and  yield  reduced  activation 
energies  for  diffusion  of  about  3-4000  cal/mole.  Diffusion  co¬ 
efficients  at  room  temperature,  as  determined  from  the 
corrected  permeabilities,  have  yielded  D~  2  *i  10  cm  /sec. 


3.  Limiting  current  densities  for  cathodic  oxygen  reduction  on 
IEM-Pt  black  type  electrolyte-electrode  systems  are  quasi- 
linearly  dependent  on  oxygen  partial  pressure,  for  the  binary 
system  O2/N2,  in  the  absence  of  convective  flow.  At  room 
temperature,  limiting  current  densities  for  air  are  about 
130  Ma/cm  .  These  current  densities  can  be  increased  by 
forced  gas  flow  over  the  electrode. 
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ABSTRACT  (Cont'd) 


In  air-breathing  channels  of  defined  geometry  current- 
potential  characteristics  are  position-dependent  and  can 
be  formulated  quantitatively  from  relationships  describing 
local  limiting  current  densities.  Channel-edge  currents 
can  become  appreciably  larger  than  center-channel  currents 
(2-3  times  for  channel-height/length  ratios  of  1/20  to  1/40). 
Diffusional  resistances  (defined  as  dE/dl)  can  be  predicted 
from  calculated  local  limiting  currents,  assuming  gas  dif¬ 
fusion  controlled  rates.  Dependence  on  channel  length, 
height  and  diffusion  coefficient  is  shown.  Experimental 
results  are  compared  to  derived  equations.  Under  defined 
flow  rate  conditions,  local  diffusion  polarization  is  appreciably 
reduced,  and  already  at  low  flow  rates,  quasi-uniform  current 
distributions  can  be  observed,  resulting  practically  in  cell 
operation  near  limiting  currents. 

4.  The  success  of  cell  operation  at  high  current  densities  depends 
essentially  on  the  study  and  design  of  the  air  electrode.  A  new 
design  of  an  air  electrode  current  collector  based  on  the  re¬ 
quirement  that  oxygen  and  water  partial  pressures  as  well  as 
temperatures  should  be  as  uniform  as  possible  in  order  to 
operate  at  uniform  current  densities  and  this  at  uniform  rates 
(oxygen. diffusion,  water  removal,  etc...)  has  resulted  in 
single  cell  operation  under  steady  state  conditions  of  60  amps/ 
ft2  at  D.  C.  voltages  about  0.  65  volt.  This  result  represents 
an  appreciable  improvement  over  experimental  data  reported 
under  earlier  contracts.  In  fact,  room  temperature  air- 
operated  fuel  cells  can  now  be  expected  to  extend  the  ran^e  of 
current  densities  up  to  100  amps/ft2  (at  about  50  watts/ft  ). 

5.  Overall  heat  and  mass  transfer  investigations,  including  de¬ 
termination  of  effective  thermal  conductivities  and  effective 
diffusion  coefficients  for  the  porous  catalytic  layer,  as  well 
as  measurements  of  rates  of  water  removal  separated  in 
thermal  contribution  by  Joule  effect  and  mass  transport  con¬ 
tribution  by  diffusion  to  bulk  gas  stream  will  allow  the  attempt 
for  quantitative  description  of  the  transport  processes. 
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ABSTRACT  (Cont'd) 


6.  The  realization  of  the  importance  of  major  parameters,  i.  e. 
oxygen  partial  pressure,  water  partial  pressure,  electrode 
temperature,  etc.,  resulting  in  surface  uniformity  regarding 
reaction  rates,  heat  and  mass  transfer  rates,  has  yielded 
various  conceptual  designs.  Consideration  has  been  given  to 
a  system  recycling  large  volumes  of  air  and  operating  at 
near -saturated  conditions.  A  considerable  fraction  of  the 
overall  heat  transferred  from  the  device  is  removed  through 
the  anode.  Such  a  design,  in  essence,  quite  different  from 
the  concept  ragarding  the  modified  air  electrode  current 
collector,  allows  for  a  broad  approach  to  battery  design 
aimed  at  the  multitude  of  applications  based  on  air -operation. 

7.  A  design  for  air-operated  fuel  cells  control  has  yielded  a 
simple  control  device  based  on  a  signal  input  corresponding 
to  positive  or  negative  rates  of  variation  of  cell  internal 
resistance.  Air  flow  modulation  allowing  for  changes  of 
rates  of  water  transport  from  the  electrode  surface  insure 
constancy  of  internal  cell  resistance  and  of  cell  voltage. 

Such  a  control  device  will  allow  cell  operation  under  largely 
variable  environmental  conditions,  without  appreciable 
variations  of  cell  power  density. 
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3.  0  PUBLICATIONS,  LECTURES,  REPORTS  AND  CONFERENCES 

3.  1  Work  planning  Conference  with  Contracting  Officer,  Fort  Monmouth, 

N.  J. 


A  review  meeting  of  the  work  accomplished  on  the  contract  up  to 
mid- January  1963  was  held  at  the  USAELRDL,  Fort  Monmouth, 

N.  J. ,  on  January  25,  1963.  The  attendance  included  Messrs. 

S.  Bartoch,  H.  Hunger  and  J.  Perry  of  the  USAELRDL  and 
Messrs.  A.  Glenn,  T.  K.  Johnson,  H.  Maget  and  E.  A.  Oster 
of  the  Direct  Energy  Conversion  Operation  of  General  Electric. 

3.  2  Work  Planning  Conference  at  the  DECO  Fuel  Cell  Laboratory, 

Lynn,  Mass. 

During  the  review  meeting  of  February  21,  1963,  held  at  the  DECO 
Fuel  Cell  Laboratory,  Lynn,  Mass.,  the  major  investigators 
associated  with  the  contract  program  presented  experimental 
results  and  planned  work  to  Mr.  S.  Bartoch  of  the  USAELRDL, 

Fort  Monmouth,  N.  J. 

Mr.  A.  Daniel  of  the  USAELRDL,  Fort  Monmouth,  New  Jersey, 
visited  the  DECO  Fuel  Cell  Laboratory  on  April  15,  1963.  DECO 
attendance  included  Messrs.  D,  W.  Craft,  T.  K.  Johnson, 

E.  A.  Oster  and  H.  Maget.  During  the  meeting,  the  status  of 
the  technical  program  as  well  as  planned  experimental  work  were 
discussed. 

3.  3  Power  Sources  Conference  Presentation 

At  the  Fuel  Cell  Batteries  Meeting  of  the  Seventeenth  Annual  Power 
Sources  Conference  held  on  May  22,  1963,  Dr.  H,  Maget  presented 
a  short,  prepared  note  regarding  some  experimental  results  obtained 
under  the  present  contract,  as  well  as  the  improvements  in  power 
capabilities  of  IEM  Air  Breathing  Fuel  Cells. 


-  7  - 


4.0 


FACTUAL  DATA 


4. 1  TASK  I 

INTERFACE  STRUCTURE  AND  INTERFACIAL 
PROPERTIES  OF  IEM/ ELECTRODE  SYSTEMS 

S.  J,  Krumbein 

4.  1.  1  Introduction 

Fabrication  of  an  ion-exchange  membrane  (IEM)  fuel  cell  results  in 
the  formation  of  a  new  interfacial  zone,  comprising  the  boundary 
region  between  the  electrode  and  the  IEM  electrolyte.  The  importance 
of  this  interface  is  known  from  the  dependence  of  the  cell's  operational 
characteristics  on  the  methods  used  for  combining  the  two  phases;  as 
well  as  by  the  great  effort  that  has  been  expended  in  developing  optimum 
manufacturing  procedures. 

As  pointed  out  in  the  first  semi-annual  progress  report,  the  amount 
of  work  accomplished  in  the  general  realm  of  solid-solid  interfaces 
has  been  quite  meager  when  compared  to  the  considerable  published 
literature  on  other  facets  of  surface  chemistry.  This  observation 
is  particularly  pertinent  to  the  complex  series  of  interfacial  regions 
inherent  in  IEM  fuel  cell  components;  and  is  due,  no  doubt,  to  their 
relative  experimental  inaccessibility  when  contrasted  with  classically 
well-defined,  partly  fluid  interfacial  systems.  Primary  emphasis 
was,  therefore,  placed  on  the  development  and  evaluation  of  represen¬ 
tative  experimental  techniques. 

Impedance  Measurements 

Introduction:  As  discussed  in  a  previous  report  ^  the  measured 
capacitance  of  an  IEM  fuel  cell  component  at  high  frequencies  is 
assumed  to  be  related  to  the  resultant  of  a  complex  parallel  circuit 
of  different  RC  ratios,  but  free  from  practically  all  faradaic  effects. 
Although  the  magnitude  of  these  limiting  capacitance  values  will  be 
far  below  the  "true"  ones  (®)  they  should,  nevertheless,  serve  as 
a  comparative  measure  of  the  relative  double  layer  capacities  for  the 
different  IEM/electrode  systems.  Hence,  if  it  were  possible  to 
maintain  all  other  cell  parameters  constant,  (viz.  catalyst  loading, 
fabrication  conditions,  electrode  thickness,  etc.)  much  could  be 
learned  regarding  relative  surface  areas  and  roughness  of  the 
interfacial  region  from  dependence  of  the  measured  impedance 
values  on  certain  of  these  properly-chosen  experimental  variables. 


4.  1.  2 
4.  1.  2.  1 
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4. 1.2.2  Experimental  Techniques  and  Procedures 

General:  The  program  in  progress  during  this  report  period  was 
pursued  with  an  eye  towards  evaluating  the  feasibility  of  obtaining 
meaningful  data  on  the  fuel  cell  interphases  by  the  A.  C.  impedance 
technique.  In  particular,  an  effort  was  made  to  characterize  the 
experiment  by  attempting  to  identify  the  particular  interface  or 
interfaces  being  measured  and  resolving  them  into  their  individual 
components. 

Two  "conducting"  interfaces  in  the  IEM  fuel  cell  can  be  visualized 
1)  the  metal/aqueous  solution  boundary  within  the  porous  electrode 
and  2)  the  metal/IEM  interfacial  region  with  the  measured  capa¬ 
citance  being  affected  by  both.  Since  these  two  components  are 
essentially  in  parallel,  the  relation 

C  =  f  (I)  +  g  (m)  (1-1) 

can  be  written,  where  I  refers  to  the  IEM/metal  interface  and  m 
is  the  amount  of  electrode  material  (e.  g.  platinum  black)  per  unit 
geometric  area.  For  electrodes  of  identical  relative  composition 
and  porosity,  g  (m)  might  be  identified  with  km  (i.e,  direct  propor¬ 
tionality),  as  a  first  order  of  approximation,  so  that  equation  (1-2) 
could  be  simplified  to: 

C  =  km  +  f  (I)  t1’2) 

In  the  course  of  this  program  two  experimental  techniques  for 
eliminating  the  g  (m)  --  or  km  --  term  were  tried.  These  were. 

1,  Varying  the  amount  of  liquid  electrolyte  (dilute  pH4  sulfonic  acid 
solution)  in  the  (porous)  electrode  phase  by  slowly  passing  dry  nitrogen 
through  the  cell's  gas  compartments.  As  drying  progresses,  the 
effective  metal/liquid  interfacial  area  will  decrease;  and  the  measured 
double  layer  capacitance  should  drop  accordingly.  It  was  postulated 
that  the  capacitance  would  finally  decrease  to  a  low,  but  finite,  limiting 
value  corresponding  to  the  IEM  interface,  since  this  region  was  not 
expected  to  dry  out  at  the  same  rate  as  the  bulk  electrode. 

2.  The  catalyst  loading  (i.e.  the  thickness  of  the  electrode)  was  varied 
but  all  other  conditions  were  held  constant.  By  this  means,  variations 
in  impedance  values  with  total  surface  area  of  the  platinum  could  be 
observed  and,  if  found  to  be  regular,  could  be  extrapolated  to  zero 
catalyst  loading  (Eqn.  1-1)  to  eliminate  the  bulk  electrode  term.  The 
success  of  the  method  hinges  on  the  possibility  of  fabricating  reprodu- 
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cible  cell  configurations  whose  only  variable  is  the  amount  of 
electrode  material  (of  constant  composition  and  porosity)  pressed 
onto  the  membrane.  A  correlary  condition  is  the  need  for  temporal 
reproducibility  so  that  cells  made  at  different  times  can  be  compared. 

Procedure:  -  In  order  to  avoid  some  of  the  ambiguities  of  assymetrical 
geometries  (2,  3)  when  using  the  bridge  technique,  all  measurements 
were  limited  to  back-to-back  parallel  configurations;  where  the 
measured  geometrical  areas  would  be  symmetrically  disposed  and 
(unless  otherwise  specified)  of  approximately  identical  character. 

For  this  purpose,  cells  were  fabricated  by  pressing  identical 
electrode  discs  onto  both  sides  of  3"  x  5"  ion-exchange  membranes  — 
each  of  which  was  designed  to  fit  into  the  smaller  half  of  a  4-compart¬ 
ment  test  fixture.  The  fixture  itself  is  described  in  detail  in  a  pre¬ 
vious  report  (Figure  2  3  of  Reference  (1)  and  was  found  particularly 
useful  due  to  the  ease  with  which  it  could  be  adapted  to  both  flooding 
and  drying  conditions. 

All  the  impedance  determinations  were  run  with  the  variable  ratio- 
transformer  bridge,  developed  by  the  Electrical  and  Physical 
Standards  Operation  of  General  Electric’s  Advanced  Technology 
Laboratory  (6).  Difficulties  were  encountered  at  the  outset  due  to 
large  inductive  residuals  of  the  circuit  components  at  the  high 
frequencies  used.  With  modifications  and  refinements  in  both 
measuring  procedure  and  circuitry,  however,  it  became  possible  to 
obtain  an  accuracy  of  10%  at  10,  000  CPS  when  encountering  capacitance 
values  up  to  3,  OOOytt-  F  (7). 

This  was  accomplished  by  the  use  of  the  substitution  technique  of 
Hillhouse  and  Kline,  provided  for  by  the  modified  circuitry  of 
Figure  I- 1.  In  practice,  the  experimental  fixture,  containing  the 
leads  and  current-collector  screens  (i.  e.  with  Cx  shorted  by 
switch,  Si),  was  placed  in  a  fixed  geometric  position  with  respect 
to  the  circuit;  and  resistance  and  tr.ansformer-ratio  values 
(Ap  or  1-A0)  were  obtained  at  all  frequencies  of  interest.  The 
experimental  cell  was  then  placed  in  the  fixture  and  another  set  of 
resistance  and  transformer-ratio  values  (Ai  or  1-Ai)  obtained. 

If  care  was  taken  to  keep  the  relative  geometric  relationships  of  all 
the  circuit  components  fixed  during  these  two  operations,  the  experi¬ 
mental  capacitance,  Cx,  could  be  calculated  by  the  relation 

C  B(  1  -A  i  }  ( l~ApC) _ 

(A1C  -  Ao>  [l-a2(Ls  +  Ld>J  t1'3* 

C  c 

where  Cg  is  the  "standard"  capacitor  (of  1/*F),  A1  and  AQ  are  the 
corrected  (for  inductance)  A^  and  AQ  values,  and  Lg  and  are  the 
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inductive  residuals  of  the  "standard  circuit"  and  of  the  balancing 
resistors  respectively.  One  should  also  note  that  neither  the 
exact  magnitude  of  nor  that  of  the  inductive  residuals  of  the 
test  circuit  (Lx)  need  to  be  known,  provided  they  remain  constant. 

In  most  cases,  the  high  resolution  of  the  method,  together  with 
the  low  Cg/Cx  ratios  generally  encountered  in  our  measurements 
(<0.01),  enables  the  use  of  the  approximate  relation 

a*=  ^  (1-a'0)2  (1-4) 

(where  Aq  is  an  average  of  the  Aq  values  from  2  to  10  Kilocycles), 
since  A/(l-A)  was  very  close  to  0.  1  at  practically  all  frequencies. 
One  could  then  be  justified  in  drawing  up  calibration  curves  for 
direct  "readout"  of  capacity  from  the  calculated  AA  values. 

4.  1.  2.  3  Results:  - 

1.  Drying  the  cell,  by  passing  "dry"  nitrogen  gas  at  low  pressure 
(2-5  psi)  through  the  small  gas  compartments  of  the  test  fixture, 
produced  a  continuous  drop  in  capacitance,  coupled  with  a  rise  in 
the  measured  cell  resistance.  In  the  initial  stages  of  the  drying 
procedure,  the  cell  was  observed  to  recover  (with  respect  to  its 
impedance  parameters)  when  the  gas  flow  was  stopped,  but 
continued  passage  of  N2  eventually  resulted  in  capacitance  values 
below  iytF/cm2  and  a  related  resistance  increase  of  several  orders 
of  magnitude.  During  this  latter  stage,  a  slight  tendency  toward 
recovery  was  also  observed  --  manifested  by  very  small  increases 
in  capacity- -due,  probably,  to  longitudinal  diffusion  of  water  from 
the  1  'electrode - le s s"  portion  of  the  ion-exchange  membrane. 

2.  Impedance  measurements  under  "normal"  wet  conditions  were 
generally  conducted  with  the  cell's  gas  compartments  flooded  with 
dilute  (pH  4)  sulfuric  acid  solution,  so  as  to  eliminate  ambiguities 
connected  with  the  position  of  the  liquid  level  in  the  electrode  phase. 
Under  these  conditions,  the  electrode  surfaces  were  at  a  potential 
of  approximately  1.  0V  (versus  NHE);  based  upon  measurements 
made  against  two  large  (3-5/8"  diameter)  back-to-back,  built-in 
hydrogen  reference  electrodes  (cf.  Figure  22  of  Reference  (1)) 

of  similar  composition. 
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The  first  set  of  experimental  determinations  --  all  of  which^were 
performed  with  cells  fabricated  from  0.005  inch  thick  GPS 
membranes  —  appeared  to  indicate  direct  variation  of  capacitance 
(at  constant  frequency)  with  electrode  loading.  That  this  result 
was  fortuitous  was  shown  by  later  measurements  on  the  same  cells. 


Tables  1-1  and  1-2  not  only  indicate  that  the  capacity  changed  with 
time,  but  that  these  changes  themselves  were  not  reproducible; 
although  they  were,  for  the  most  part,  in  the  general  direction  of 
increasing  capacitance  and  resistance. 

TABLE  1-1 

Constant-Freauencv  Impedance  Characteristics  of  IEM  Fuel  Cell  Configurations, 

Electrode 

Date 

C5Kc 

R5Kc 

Cell 

Loading 

Date 

Run 

F/crn 

Approximately 

No.  ** *** 

gm/crn 

Made 

1963 

,  ohms 

37 

0.025 

1/3/63 

1/31 

215 

0.  17 

30 

0.021 

12/21/62 

1/17 

155 

0.73* 

36 

0.  013 

12/24/62 

1/31 

85 

0.  15 

4/26 

92 

0.  49 

32 

0. 0074 

12/21/62 

1/16 

78 

1.  39* 

1/31 

111 

0.97 

4/25 

311 

1.05 

*  Not  flooded  during  these  two  runs.  The  cells  had  been  soaked  in  pH4  (H2S04) 
just  prior  to  being  placed  in  the  fixture,  however. 

##  All  cells  in  table  are  based  on  GPS  membranes  of  0.  005  inch  thickness. 


***In  this  text  GPS  will  represent  grafted  polystyrene  sulfonic  acid. 
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TABLE  1-2 

Constant-Frequency  Impedance  Characteristics  of  IEM  Fuel  Cell  Configurations 


Electrode 

Date 

Date 

C 

R 

Cell 

No.* 

Loading 
gm/ cm^ 

Made 

1963 

Run 

1963 

F/cm2 

5Kc 

8Kc 

5-10Kc 

57 

0.029 

3/18 

3/22 

+ 

680'3 

495 

0.20 

51-1 

0.0165 

3/11 

3/21 

686*1 

485 

0.  14 

3/22 

+ 

800"4 

630 

0.  19 

53 

** 

0.015 

3/12 

3/21 

+ 

845‘3 

610 

0.  14 

50 

0.0095 

3/8 

4/26 

1168*38 

940 

0.21 

54 

0. 008c 
b 

3/13 

3/21 

862*0 

675 

0.  18 

3/21 

+ 

788“0 

600 

0.  20 

4/25 

+ 

1 190” 18 

1030 

0.26 

56 

0.005 

3.  14 

3/22 

385*2 

278 

0.  32 

4/26 

+ 

500” 10 

412 

0.62 

*** 

59 

0.  009 

3/19 

3/22 

+ 

322'3 

249 

■» 

o 

OJ 

00 

*  All  cells  (except  No.  59)  in  this  table  are  based  on  GPS  membranes  of 
0.010  inch  thickness. 

**  Electrodes  for  No.  51-1  and  53  were  made  with  7  mil  and  10  mil  decals 
respectively;  No.  51-1  was,  therefore,  of  greater  compactness. 

PSS  membrane 
Average  deviation 
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Measurements  were  also  obtained  with  pre-pressed  electrode  configur¬ 
ations.  These  were  either  applied  with  heat  and  pressure  to  the 
membrane,  so  that  the  resulting  cell  was  outwardly  no  different  from 
those  of  Table  1-1  and  1-2,  or  else  consisted  of  individual  discs  (^-5/8" 
diameter  and  0.  005  inches  thick)  containing  about  0.  065  gms/cm  of 
catalyst;  and  held  against  the  IEM  by  simple  mechanical  pressure  of 
the  current  collectors. 

The  former  (with  a  0.  005  inch  thick  GPS  membrane)  yielded  capacitance 
values  of  approximately  the  same  magnitude  as  those  of  Table  1-1,  but 
showed  definite  tendencies  towards  increasing  capacitance  and  decreasing 
resistance  over  the  half-hour  period  of  the  actual  determinations  (from 
low  to  high  frequencies  and  back  --  cf.  Figure  1-2).  The  latter,  how¬ 
ever,  gave  lower  capacity  values  but  a  greater  tendency  for  change 
both  in  between  runs  (as  in  Tables  1-1  and  1-2)  and  during  the  course 
of  the  experimental  series.  These  changes  were  always  toward  higher 
capacities  and  lower  resistance,  and  are  illustrated,  to  a  certain 
extent,  by  the  broken  lines  of  figure  1-2. 

Auxiliary  determinations  were  conducted  with  the  same  two  GPS 
membrane  types,  but  with  platinum  foil  discs  (1-5/8  inches  diameter 
and  0.001  inch  thickness)  in  place  of  the  pre-pressed  electrodes. 

Although  time -dependence  during  the  measurements  was  similar, 
the  measured  capacitance  values  were  found  to  be  between  6  and  8 
/tF/cm^,  for  both  types  of  membrane,  over  the  same  frequency  range 

as  above. 

3.  Variations  of  capacity  with  frequency  were  virtually  independent 
of  all  experimental  parameters;  and  was  shown  to  be  so  even  during 
early  stages  of  drying  (Figure  1-3).  The  logarithmic  equivalents  of 
this  typical  hyperbolic  dependence  is  illustrated  by  the  plots  of 
Figures  1-2  and  1-3. 

4,  l,  2.  4  Discussion  and  Conclusions: 

1.  The  capacitance  values  are  apparently  a  relative  measure  of  the 
"active",  non-wetproofed,  interfacial  area  between  the  platinum 
black  and  aqueous  phases  (Figure  1-4).  This  would  include  only 
those  regions  of  free  platinum  surface  which  are  in  electrolytic 
contact,  through  the  aqueous  phase,  with  the  ion  exchange  mem¬ 
brane  (i.e.  not  in  blocked  pores).  As  such,  the  measurements 
could  not  differentiate  between  the  platinum  surfaces  within  the 
bulk  porous  electrode  phase  and  those  directly  adjoining  the  membrane 
material.  This  was  experimentally  corroborated  by  the  comparatively 
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fast  attainment  of  water  equilibrium  between  IEM  and  electrode  --as 
shown  by  the  tendency  toward  recovery  of  the  impedance  parameters 
during  the  drying  experiment  -  -so  that  the  assumptions  of  Section 
4.  1.  2.  2  would  appear  to  have  been  incorrect. 

2.  Capacitance  values  were  time-dependent,  but  the  magnitude  of 
the  changes  could  not  be  reproducibly  related  to  known  experimental 
parameters.  Considerable  change  occurred  even  when  the  cells  were 
stored  under  conditions  which  would  preclude  drying  out.  Judging 
from  the  odor  released  when  the  water-tight  bags  were  opened  (after 
a  month  or  more),  these  changes  may  be  partly  associated  with 
membrane  degradation.  Since  the  observed  time-dependence  was 
generally  in  the  direction  of  increasing  capacitance,  any  effect  due 
to  degradation  products  adsorbed  onto  the  platinum  surface  (a  process 
that  has  long  been  known  (18)  to  lower  the  double-layer  capacity) 
would  appear  to  have  been  masked  by  the  directionally  opposite  effects 
of  structural  changes  in  the  experimental  cell. 

Degradation  of  the  membrane,  as  well  as  shrinkage  and  swelling 
(see  Task  II),  would, however,  tend  to  produce  distortions  and/or 
separations  between  the  IEM  surface  and  parts  of  the  electrode 
material,  with  resultant  changes  in  the  interfacial  area.  Of  equal 
importance  would  be  the  gradual  release  of  the  original  compressional 
stresses  (probably  aided  by  water  absorption)  in  the  porous  electrodes  -- 
especially  those  that  had  been  applied,  in  the  powder  form,  directly 
to  the  IEM  (Tables  I  -  1  and  1-2).  This  would  lead  to  an  increase  in 
the  open  pore  network  connecting  with  the  membrane  surface  (Figure 
I  -  4),  with  a  concommitant  increase  in  the  total  area  of  active  plati¬ 
num/water  interface. 

Some  degree  of  evidence  for  this  hypothesis  can  be  seen  from  the 
capacity  values  of  cells  No.  51-1  and  53  (Table  1-2).  The  electrodes 
of  both  had  approximately  the  same  loading  factor,  but  differed  con¬ 
siderably  in  their  compactness  during  the  pressing  stage  of  cell  fabri¬ 
cation  --  so  that  initial  capacitance  values  for  cell  No.  51-1  were 
considerably  below  those  for  No.  53.  In  the  next  day's  run  (with  soaking 
in-between),  however,  values  quite  close  to  the  latter  were  produced, 


*  Figure  1-4  represents  one  opinion  of  the  electrode  structure  in  the 
region  adjacent  to  the  IEM.  It  is  intended  only  as  a  pictorial  aid  to  the 
discussion,  and  should  not  be  considered  a  scale  model  representation 
of  the  true  cross-section. 
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indicating  a  comparatively  rapid  opening  up  of  the  pore  structure 
(cf.  the  two  separate  runs  of  cell  No.  54). 


Figure  1-2  indicates  that  some  such  effect  may  have  also  occurred 
with  the  pre-pressed  electrode  discs,  after  having  been  soaked  in 
the  aqueous  solution  for  some  time.  It  should  be  noted,  in  this 
context,  that  different  0.  010-inch  thick  membrane  samples  were 
used  in  May  and  April,  so  that  the  great  differences  between  the 
two  months'  runs  --  and  the  small  effect  of  substituting  a  0.005-inch 
thick  membrane  --  was  chiefly  due  to  differences  in  the  electrodes 
themselves.  The  effect  of  water  on  the  porosity  is  also  suggested  by 
the  short-tirne  changes  in  this  figure,  the  points  on  the  dashed  lines 
being  those  obtained  at  the  end  of  each  series  of  rune. 


3.  The  results  of  this  phase  suggest  the  use  of  the  impedance  tech¬ 
nique  in  determining  aging  and/or  time -dependent  operational  effects 
occurring  at  the  IEM/electrode  interface.  Determination  of  relative 
interfacial  areas  and  other  structural  parameters  of  different  IEM 
fuel  cells,  however,  can,  at  best,  only  be  of  a  qualitative  nature; 
unless  improved  fabrication  and  storage  procedures  can  be  developed 
for  assuring  better  microscopic  reproducibility. 


4,  1.  3  Adhesion  Measurements 


4.  1.  3.  1  Theory: 

A  program  has  been  undertaken  to  determine  the  relative  adhesion 
existing  between  the  catalyst  and  membrane  phases  of  IEM  fuel  cells. 

The  method  adopted  for  initial  evaluation  was  that  developed  by 
Heavens  (9)  and,  later,  put  on  a  semi-theoretical  basis  by  Benjamin 
and  Weaver  (10).  Photograph  1-1  illustrates  the  complete  experi¬ 
mental  setup,  including  a  modified  form  of  Heavens-type  apparatus  built 
in  our  laboratory.  This  is  seen  to  consist  essentially  of  a  horizontal 
beam  with  a  weight- pan  in  its  center  and  a  smoothly  rounded  point  at 
its  end.  Samples  are  drawn  across  a  flat  table,  directly  under  the 
"point",  so  that  the  latter  travels  across  the  metal  (i.e.  electrode) 
surface.  Weights  are  added  to  the  pan  after  each  run,  until  --  at  a 
certain  critical  weight  --  the  film  should  be  cleanly  stripped  from  the 
surface. 

The  theory  of  the  experiment  is  based  upon  the  work  of  Bowden  and 
Tabor  (11,  12)  on  the  adhesional  theory  of  friction;  under  which 
Amonton's  second  law  acquires  the  form 

F  r  =  Ws/p  f1'5) 
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where  Fria  the  frictional  force,  W  the  normal  load,  and  p  the  yield 
pressure  of  the  softer  material,  while  s  is  the  force  per  unit  area, 
acting  in  a  direction  tangential  to  the  interface,  which  is  required 
to  shear  the  adhesive  bonds. 

It  has  been  shown  (10),  however,  that  the  "critical"  vertical  forces 
obtained  with  the  Heavens  apparatus  are  determined,  primarily, 
by  the  properties  of  the  film/substrate  interface.  The  hard  point 
(or  ball)  is  assumed  to  slide  over  the  metal  film  for  increasing  loads, 
up  to  the  point  where  deformation  of  the  substrate  causes  a  loosening 
of  the  film.  The  load  at  which  the  deformation  occurs  can  be  related 
to  the  yield  pressure  of  the  substrate  (11)  by 

W  =  Ap  (1-6) 

where  A'  is  the  area  of  contact  between  the  (hard)  ball  or  point  and 
the  softer  substrate. 


Any  effect  of  the  metal  film  is  assumed  to  be  negligible  provided  its 
thickness  is  much  less  than  the  ball's  radius,  so  that  the  former  can 
readily  be  deformed  to  the  contour  of  the  indentation.  The  maximum 
shearing  forces  between  this  film  and  the  substrate  surface  would 
then  be  expected  to  occur  at  the  lip  o i  the  identation.  It  is  here 
that  the  adhesive  bonds  would  be  initially  broken  by  the  critical 
vertical  load;  and  the  horizontal  motion  of  the  indentor  would 
merely  serve  to  push  the  sheared  film  aside. 

From  equation  (1-6)  and  analysis  of  forces  (including  these  two 
assumptions)  Benjamin  and  Weaver  (10)  have  shown  that  the  shearing 
force  (per  unit  area)  could  be  expressed  as: 


(trr2P  -  W)  l'Z 

where  r  is  the  radius  of  curvature  of  the  indentor  and  P  is  now  the 
deformation  or  indentation  hardness  of  the  substrate. 

It  is  seen  that,  for  identical  substrate  materials  (identical  ion-exchange 
membranes,  in  our  case),  comparative  measurements  of  adhesion  can 
be  obtained  from  the  vertical  load,  W,  alone;  although  this  parameter 
will  also  be  a  function  of  the  fraction  of  the  interfacial  area  on  which 
actual  IEM-electrode  bonding  has  taken  place.  The  hardness  of  the 
substrate,  P,  must  be  known,  however,  before  one  can  calculate  the 
shearing  forces.  This  was  determined  for  both  thick  GPS  and  PSS 

)j( 

non  cross-linked  sulfonated  polystyrene 
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membranes  by  measuring  the  apparent  compressive  yield  stress, 

Y,  by  standard  technique  (15),  and  applying  the  relation  (13) 

P  =  cY  (1-8) 


where  c  is  between  2.  5  and  3.  The  lower  value  was  chosen,  as 
proposed  by  Benjamin  and  Weaver,  and  the  experimental  value  of 
P  was  found  to  be  approximately  11  Kg/mm^  for  the  fluorocarbon 


based  GPS  membrane,  which  can  be  compared  to  the  value  of 

11.  5  Kg/mm^  measured  directly  by  King  and  Tabor  on  pure  Kel-F  (12). 


Experimental  compressive  yield  stress  values  for  the  PSS  membrane 
were  not  reproducible,  but  appeared  to  average  about  3  Kg/mm^  -- 
for  an  apparent  indentation  hardness  of  7.  5  Kg/mm^. 


4.  1.  3.  2  Experimental: 


In  addition  to  the  assumptions  that  were  specified  in  Section  4.  1.3.  1, 
certain  other  assumptions  and  experimental  conditions  must  be  met 
for  the  theory  to  be  tenable. 


1.  The  indentation  hardness,  P,  of  the  substrate  must  not  change 
significantly  with  W  or  the  time  of  actual  contact  with  the  indentor. 
Approximate  constancy  of  P  has  been  found  by  King  and  Tabor  (14) 
over  a  10  to  40- second  range  of  "loading"  time  --  with  the  least 
changes  for  the  shorter  times.  Nevertheless,  it  was  found  advisable 
to  make  all  measurements  at  constant  horizontal  pulling  speed  with  a 
Dillon  Universal  Tester  (Photograph  I  -  1),  so  as  to  eliminate  any 
additional  unknown  parameters  from  the  comparative  measurements. 
A  speed  of  0.075  cm/ sec  was  chosen,  corresponding  to  a  far  shorter 
contact  time  than  the  lowest  --  and  most  accurate  --  loading  time  of 
these  authors. 


2.  The  mechanical  integrity  of  the  bulk  substrate  material  must  not 
change  during  the  experiment.  The  substrate  must  remain  flat  so  that 
the  applied  horizontal  and  vertical  forces  will  subtend  a  constant, 
known  angle  (preferably  one  of  90  degrees).  Likewise,  the  other  bulk 
properties  of  the  material  should  also  be  unaffected  by  what  is  happening 
at  its  upper  surface.  The  thickness  of  the  piece  must,  therefore,  be 
considerably  greater  than  the  ball  radius,  especially  for  our  fairly 
soft  IEM  materials.  This  factor,  together  with  the  necessity  for 
examining  the  scratches  by  transmitted  light,  required  rather  artifi¬ 
cial  fuel  cell  configurations,  consisting  of  a  thick  membrane  with 
the  electrode  material  on  only  one  side. 
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3.  Equation  (I  -  7),  in  the  form 

W  =  lTr2  PF2  (1-9) 

pZ  +  F^ 

indicates  that  the  critical  load,  W,  is  directly  proportional  to  r2. 
Moreover,  the  variation  of  W  with  F  (assuming  constant  P) 

dW  _  2Tr2P3F  (I- 10) 

dF  (p2  +  F^T2 

shows  that,  with  soft  substrate  materials,  reproducible  points  of 
comparatively  large  radius  of  curvature  will  be  required  in  order 
to  obtain  sufficient  precision  and  ease  of  measurement.  Experi¬ 
mental  indentors  were  consequently  made  from  precision  bearing- 
balls,  rigidly  fixed  to  threaded  standard  stock  (see  Photograph 
I  -  8  in  Appendix  I  -  1). 

4.  1.3.  3  Results  and  Discussion: 

1.  Significant  differences  were  found  in  the  vertical  load  required 
to  strip  electrode  materials  of  different  formulations  from  a 
specially  prepared,  air-dried,  1 /8-inch  GPS  membrane.  Pure 
platinum  black  required  a  load  of  about  900  grams  (for  a  1 /64-inch 
diameter  tungsten-carbide,  ball-bearing  "point"),  while  a  platinum 
black  -  TFE  (Tetrafluoro  Ethylene)  mixture  remained  adherant,  over 
most  of  the  test  area,  at  loads  exceeding  1200  grams.  When  the 
pieces  were  soaked  in  water  just  prior  to  the  measurement,  the 
required  load  was  considerably  reduced  for  the  pure  metal  electrode, 
whereas  no  significant  changes  were  observed  for  the  fluorocarbon 
mixture.  On  the  other  hand,  PSS  substrates,  with  the  same  platinum/ 
TFE  electrodes,  yielded  lower  critical  force  values  --  which  decreased 
still  further  under  wet  conditions. 

These  results  were  obtained  with  fairly  uniform  electrode  coatings; 
and  are  summarized  in  Table  1-3. 
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TABLE  I  -  3 

Vertical  Loads  and  Shearing  (Adhesional)  Forces 
For  IEM/Electrode  Systems. 


IEM 

Electrode 

Condition 

W 

F 

Substrate 

Composition 

of  piece 

iKgL 

(10°  dynes/cm‘) 

GPS 

Pure  Pt 

Dry 

0.9 

15 

Wet 

0.  6 

10 

GPS 

Pt/TFE 

Dry 

1.2 

28 

Wet 

1.2 

28 

PSS 

Pt/TFE 

Dry 

0.7 

13* 

Wet 

0.4 

6  * 

*  Based  upon  "average" 

values  of  poorly- 

-defined  indentation  hardness. 

2.  Where  the  electrode  material  was  applied  unevenly  to  the  membrane 
due  to  local  differences  in  preparation  and  pressing --  the  "critical" 
vertical  weight  varied  with  apparent  thickness  of  the  layer.  In  fact, 
by  applying  a  900  gram  weight  over  various  portions  of  this  electrode 
it  was  possible  to  obtain  a  set  of  "contour"  lines  ranging  from  strong 
to  weak  (completely  stripped)  regions.  Since  previous  theoretical 
and  experimental  work  (10)  has  shown  that  the  critical  load  --  for  an 
ideally  coherant  coating--  should  be  independent  of  film  thickness, 
the  observed  adhesion  versus  thickness  dependence  would  appear  to 
reflect  an  increase  in  adhesive  bond  density  (between  the  1EM  and  the 
aggregated  electrode  material)  with  thickness. 

3.  Micrographs  of  "stripped"  grooves  from  this  type  of  sample  are 
shown  in  Photograph  1-2  for  the  900-gram  load,  as  well  as  in 
Photograph  1-3  for  a  600-gram  weight  on  a  region  of  weak  bonding. 

In  addition,  Photograph  1-9  (see  Appendix  1-1)  shown  the  furrow 
resulting  from  the  1200- gram  critical  load  that  was  required  for 
a  uniformly-pressed  Pt/TFE  electrode  (Table  I  -  3).  Under  reflected 
light,  most  of  the  stripped  material  was  found  to  have  been  removed 
from  the  inner  lips  of  the  furrow  --  thereby  pointing  up  the  applicability 
of  the  underlying  theory  (10)--  and  pushed  to  the  sides  as  fairly  large 
pieces  (photograph  I- 10;  Appendix  1-1). 
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As  shown  by  the  three  previous  photographs,  not  all  the  electrode 
material  was  ploughed  to  the  side  after  the  initial  adhesive  bonds 
were  broken.  Instead,  some  remained  imbedded  in  the  bottom 
of  the  groove,  the  quantity  of  such  "debris"  increasing  with  the 
critical  vertical  load  and  its  associated  adhesive  forces.  That 
this  effect  is  related  to  the  nature  of  the  adhesive  bonding  can 
be  seen  by  comparing  these  pictures  to  those  for  the  pure  plati¬ 
num  black  electrode.  Photograph  1-4  {  a  close-up  to  Photograph 
I- 11  of  the  Appendix)  shows  the  removal  of  a  much  greater  amount 
of  material--  by  an  800-gram  weight--  from  the  bottom  of  the 
furrow,  with  whole  sections  of  electrode  cleanly  stripped  away. 
Moreover,  prior  soaking  of  the  piece  in  water  served  to  accentuate 
this  effect,  resulting  in  an, almost  completely  clean  scratch  {with 
700  grams;  Photograph  1-5). 

4.  The  motion  of  the  indentor,  once  it  had  penetrated  to  the 
membrane  surface,  was  of  the  "stick-slip"  (13)  variety.  This 
is  clearly  revealed  in  Photographs  1-4,  I- 11,  and  (to  a  lesser 
extent)  1-2;  and  is  to  be  expected  for  the  fluorocarbon  based  GPS 
membrane  used  (14).  This  behavior  tended  to  give  way  to  a 
partial  penetration  of  the  "point"  into  the  material  even  when 
dry  (Photograph  I-  12,  showing  another  part  of  the  same  groove 
as  I- 11);  and  led,  in  some  cases,  to  actual  jamming  of  the  point 
when  the  substrate  was  wet.  The  latter  effect  is  illustrated  in 
Photograph  1-13  of  Appendix  1-1  for  another  section  of  the  same 
scratch  as  1-5.  Under  reflected  light  (Photograph  I- 14),  this 
path  of  the  indenting  ball  into  the  membrane  can  be  easily  followed 
as  a  black  "finger",  with  the  ploughed-up  catalyst  at  the  furrow's 
edges  revealing  a  comparatively  non-aggregated  structure 
(Photograph  I-  10). 

5.  The  problem  of  jamming  under  wet  conditions  was  especially 
critical  with  the  non-crosslinked  PSS  substrates,  even  when 
coated  with  the  same  platinum  black/T^E  mixture.  In  many  cases, 
the  hole  became  so  deep  that  the  "Dillon  Tester"  had  to  be  stopped 
and  reversed  to  avoid  snapping  the  pulley-cord;  and  the  bail-point 
pulled  out  by  hand. 

The  low  cohesional  strength  of  the  membrane  and  its  extreme  soft¬ 
ness  under  load  was  also  responsible  for  the  tendency  of  the  indentor 
to  strip  some  of  the  underlying  substrate  material  along  with  the 
electrode  layer.  This  took  place  even  when  the  test  piece  was  not^ 
pre- soaked;  and  would  appear  to  cast  doubt  on  the  applicability  of  the 
method  (along  with  its  underlying  assumptions)  to  PSS  membranes. 
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It  might  also  be  mentioned  that  these  pieces  tended  to  warp  when  dry; 
and  that  relatively  prolonged  air-drying  (i.e.  during  the  measure¬ 
ment)  tended  to  produce  local  blistering  and  loosening  of  the  electrode 
layer. 

6.  The  effect  of  polymeric  binders  on  the  cohesional  characteristics 
of  the  electrode  is  illustrated  by  the  reflected-light  micrographs  of 
some  grooves  made  with  less  than  the  critical  load.  Photograph  1-6 
was  made  by  950  grams  on  a  platinum  black/ TFE  electrode,  and 
shows  a  groove  whose  surface  is  little  different  in  appearance  from 
that  of  an  unscratched  electrode.  On  the  other  hand,  a  scratch  made 
on  a  platinum  black  electrode  (Photograph  1-7)  by  500  grams  reveals 
a  series  of  deep  cracks,  at  right  angles  to  the  groove,  even  though 
the  unscratched  surface  is  no  different  from  Photograph  1-6.  It 
should  also  be  noted  that  the  "unstripped"  grooves  in  Pt/TFE  material 
on  PSS  substrates  gave  the  same  appearance  as  that  for  the  GPS 
membrane  (Photograph  1-6). 

4.  1.3.  4  Conclusions: 

The  importance  of  polymeric  binders,  such  as  the  polyolefines  and 
their  halogenated  homologues,  as  wetproofing  agents  for  porous 
electrodes  has  long  been  recognized  (4,5).  More  recently,  a 
definite  parallel  was  demonstrated  in  this  laboratory  (2)  between 
optimum  fuel  cell  characteristics  and  maximum  contact  angle 
as  a  function  of  the  relative  composition  of  platinum  black  and 
Teflon  TFE  binder.  The  results  just  reported  --  using  these  same 
model  materials-  -  have  now  shown  that  hydrophobic  binders  can 
serve  the  additional  function  of  acting  as  adhesive  agents  between 
the  catalyst  and  ion-exchange  membrane. 

This  empirical  view  of  the  triply-functional  character  of  "electrode 
binders"  or  "wetproofing  agents"  can  be  best  understood  in  the  light 
of  the  very  recent  work  of  Sharpe  and  Schonhorn  (21)  on  the  adsorption 
theory  of  adhesion.  According  to  their  proposal  (which  is  outlined 
briefly  in  Appendix  1-2),  there  exists  a  one-to-one  correspondence 
between  thermodynamic  spreading  and  strong  adhesion,  so  that  the 
operational  condition  for  strong  adhesion  can  be  shown  to  be  identical 
with  Zisman's  wetting  relation  (20), 


r*  *  & 


(I- 11) 
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The  critical  surface  tension,  f  c,  of  a  solid  has  been  defined  by 
Zisman  (19)  as  the  maximum  surface  tension  that  a  liquid,  L, 
can  possess  in  order  to  just  wet  that  particular  solid.  Materials 
with  lower  surface  tensions  would  then  be  expected  to  make  good 
adhesives  for  such  solids  --  provided  they  could  be  applied  in 
fluid  form.  This  latter  condition  arises  from  the  fundamental 
mechanism  of  simple  adhesion,  which  is  postulated  (21)  as  being 
almost  exclusively  due  to  Van  der  Waals  forces,  and,  therefore, 
only  operative  over  very  small  distances. 


It  can  now  be  seen  that  the  very  same  property  that  makes  a 
material  a  good  "wetproofing"  agent  while  solid  (i.  e.  ita  low 
surface  free  energy)  also  serves  to  make  it  a  good  adhesive 
when  it  is  caused  to  undergo  a  certain  amount  of  flow.  In 
practice,  this  is  accomplished  by  the  high  pressures  and  moder¬ 
ately  high  temperatures  applied  during  fabrication,  so  that 
Teflon  TFE  t^c  =  18.  5)  ,  for  example,  would  not  only  act  as  an 

adhesive  between  the  high  surface-energy,  platinum  black  particles 
(cohesional  or  "cementing"  agent)  but  also  between  the  latter  and 
the  low  surface-energy  (}*c  =  33)  ,  fluorocarbon  based,  solid  substrate. 
In  this  context,  one  might  also  note  that,  whereas  the  softening  point 
of  Teflon  TFE  is  almost  twice  that  of  the  base,  fluorocarbon,  its  static 
yield  pressure,  at  20°C,  (5.  3Kg/mm^)  is  less  than  one-half  that  of 
the  latter  (14).  It  is,  therefore,  not  surprising  that  the  large  vertical  ^ 
forces  encountered  in  our  measurements  --  of  the  order  of  TO  Kg/mm  -- 
should  be  expected  to  cause  individual  platinum/Teflon  aggregates  that 
are  directly  under  the  ball-point  (i.e.  at  the  bottom  of  the  groove)  to 
adhere  even  more  strongly  to  the  substrate  (compare  Photographs  1-2, 
1-3,  and  I-  1 1  with  Photographs  1-4  and  1-5). 


4.  1.3.  5  Future  Work 

Further  work  is  anticipated  on  this  phase  of  the  Program,  for  the 
purpose  of  testing  our  conclusions  with  hydrophobic  binders  possessing 
critical  surface  tension  values  midway  between  those  of  TFE  and  the 
base  fluorocarbon.  Determinations  might  also  be  attempted  with 
similar  electrodes  containing  Teflon  FEP  (^ c  =  16)  as  binders. 


*  All  values  are  based  on  those  of  Zisman  and  co-workers  (19,  20). 
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4.2  TASK  II 

WATER  TRANSPORT  IN  ION  EXCHANGE  MEMBRANES 

J.  R.  Boyack 


4.2.1  Introduction 

Water  transport  phenomena  are  important  in  ion  exchange  membrane 
fuel  cell  technology  because  of  the  strong  dependence  of  membrane 
conductivity  on  water  content  and  also  from  the  standpoint  of 
controlling  the  water  balance  in  the  entire  operating  unit. 

The  driving  forces  which  affect  water  transport  are  electric 
potential  gradients,  water  activity  gradients,  hydrodynamic 
pressure  gradients,  and  temperature  gradients.  The  first  of  these 
has  been  treated  quite  extensively  in  the  literature,  (231,  (24),  (25), 
(26),  (27)  but  there  has  been  little  done  concerning  the  others,  (28) 

The  objective  of  this  task  is  to  obtain  additional  knowledge  in  the 
areas  mentioned  above.  Reported  here  are  experimental  results 
from  the  measurement  of  the  hydrodynamic  permeabilities  of  a 
series  of  polystyrene  sulfonic  acid  membranes  of  various  ion  ex¬ 
change  capacities. 

4.  2.  2  Experimental  Procedure 

The  flux  of  water  through  the  membrane  under  a  pressure  gradient  was 
determined  using  a  Stabin  osmometer,,  (29).  The  cell  was  fitted  with 
ion  exchange  membranes,  filled  with  water  and  placed  in  a  bath  of 
water.  Pressure  was  applied  through  the  capillary  tube  and  the  flux 
determined  by  observing  the  movement  of  the  meniscus  in  the  capi  ary 
as  a  function  of  time.  Membranes  used  were  polystyrene  sulfonic 
acid  having  ion  exchange  capacities  of  1. 65,  1. 80,  and  2.  15  meq/g. 

The  membranes  were  allowed  to  equilibrate  with  water  at  the  temp¬ 
erature  used  for  at  least  sixteen  hours  before  a  determination  was 

made. 

Membrane  water  content  was  measured  by  weighing  after  blotting 
off  surface  water  of  a  membrane  strip  which  had  been  equilibrated 
with  water  at  the  temperature  in  question.  The  dry  weight  was 
measured  after  drying  the  strip  in  an  oven  at  110°C. 
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4.  2.  3  Experimental  Results 

The  curves  obtained  by  plotting  the  volume  of  water  which  had  passed 
through  the  membrane  vs.  the  time  measured  from  application  of 
the  pressure  generally  had  the  form  depicted  in  Figure  II- 1.  After 
the  induction  period  of  from  five  to  thirty  minutes,  the  plots  became 
linear.  The  water  flux  was  calculated  from  the  slope  of  the  linear 
portion  of  the  curve.  The  induction  period  shown  in  Figure  II- 1 
did  not  always  occur.  A  few  of  the  plots  were  linear  right  from  the 
beginning.  The  time  of  the  induction  period  was  characteristic  in  a 
general  way  of  the  membrane  and  the  temperature  used,  but  varied 
sometimes  over  wide  limits. 

The  water  flux  through  the  membranes  was  in  all  cases  a  linear 
function  of  the  applied  pressure  as  can  be  Been  from  Figure  II-2. 

Table  II- 1  summarizes  the  permeabilities  obtained,  along  with 
other  calculated  parameters  to  be  discussed  below.  The  perme¬ 
abilities  were  calculated  from  the  least  squares  slopes  of  the 
pin*:,  in  Figure  II- 2. 

4.  2.  4  Discussion 

There  are  two  mechanisms  which  may  contribute  to  flow  of  water 
through  a  membrane  under  a  pressure  gradient:  (1)  diffuBional 
transport,  in  which  molecules  move  as  Bingle  units  and  the  influence 
of  the  pressure  is  to  raise  the  water  activity  on  one  side  of  the 
membrane,  and  (2)  convection,  in  which  bulk  flow  of  groups  of 
molecules  occurs  in  the  membrane  pores.  In  most  membrane 
systems,  both  mechanism  contribute  to  the  total  flux. 

The  flux,  J,  of  water  through  a  membrane  is  given  by 

J  =  D  dc/dx  +  Ph  Ap/t  (II_1) 

where  D  is  the  diffusion  coefficient,  dc/dx  the  concentration  gradient, 
P,  the  hydraulic  permeability,  Ap  the  pressure  difference  between 
the  two  sides  of  the  membrane,  and  t  the  membrane  thickness.  The 
ratio  of  the 
brane  is 

al/ao 
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where  V  ie  the  molal  volume  of  water,  R  is  the  gas  constant, 
and  T  the  absolute  temperature.  Assuming  equal  water  activity 
coefficients  at  the  two  faces  gives 

crc0  .  cQ  [\«p  liAEl  -IL  CoVAP 

10  L  RT  RT 

Replacement  of  dc/dx  by  (Cj-CQ)/t  yields 
C0  V  Ap 

J  =  (°  — -  +Ph)  —  <n‘4> 

RT  t 

or 

P  =  D  C°  V  +  Ph  (II-5) 

a  RT 


where  Pa  is  the  measured  apparent  permeability. 

The  induction  period  shown  in  Figure  II- 1  which  occurred  in  the 
majority  of  cases  may  be  associated  with  the  diffusion  coefficient 
of  water  in  the  membrane.  The  diffusion  coefficient  can  be  calcu¬ 
lated  approximately  from  (30) 

D  =  0.42  t Zlj  (II- 6) 

where  T  is  the  time  required  to  reach  steady  state  conditions. 

As  mentioned  before,  T  varied  widely  in  a  range  of  *30%  from  the 
average  value  for  a  given  set  of  conditions.  Nevertheless,  Eq.  (II - o ) 
yields  order  of  magnitude  values  for  D,  which  are  tabulated  in 
Table  II- 1. 

The  first  term  in  Eq.  (II - 5)  represents  the  contribution  of  the  diffu- 
sional  mechanism  to  the  total  observed  permeability.  Table  II- 1 
shows  values  for  this  term  under  the  heading  "Diffusion  Permeability." 
Comparison  with  the  measured  permeability  shows  that  only  about 
1%  of  the  total  flow  iB  by  this  mechanism. 

In  migrating  through  the  polymer  matrix,  the  water  is  forced  to  follow 
a  tortuous  path.  One  may  define  a  tortuosity  factor, 'T  ,  as 

•  T  =  t'/t  (H-7) 
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where  tf  is  the  actual  average  distance  traveled  in  migrating  through 
the  thickness  of  the  membrane.  A  simple  theory  presented  by 
Mackay  and  Meares  (31)  gives  the  tortuosity  as 

'T  =  (1  +/  )  /  (1  -0 )  (n-®) 

where  0  -  1-CQ  is  the  fraction  of  the  space  filled  by  polymer.  A 
normalized  permeability,  Pjg,  and  a  normalized  diffusion  coefficient, 
Dj,j,  are  given  by 

PM  =  p  (7,2/c  )  (n-9) 

N  a  o 

and 

dn  =  dT2  (ii-io) 

PNand  Djsj  are  the  values  which  would  be  expected  if  all  purely  geo¬ 
metric  effects  of  the  polymer  matrix  were  absent.  These  quantities 
also  appear  in  Table  II- 1.  One  should  note  that  the  DNare  all  very 
close  to  the  self  diffusion  coefficient  of  water,  which  is  2  to  5  x  10” 
cm2/sec.  (32)  over  the  temperature  range  considered.  The 
geometric  effect  of  the  matrix  alBO  accounts  for  most  of  the  variation 
in  permeability  of  the  three  membranes.  Since  hydraulic  flow 
•  resistance  is  a  function  of  the  square  of  the  effective  average  pore 
radius,  this  would  seem  to  indicate  that  the  pore  size  in  the  three 
membranes  is  about  the  same  with  differences  in  swelling  reflecting 
differences  in  the  number  of  pores. 

If  Pjvjis  assumed  to  obey  the  relationship 

p  _  P°  AE*/RT  (II- 1 1) 

*N  N 

A  plot  of  log  PN  vs.  1/T  yields  a  straight  line,  the  slope  of  which 
determines  the  activation  energy,  &E  “ .  For  the  three  membranes 
studied  here,  these  plots  are  linear  except  for  the  points  at  55  C. 
The  activation  energies  obtained  are  2700,  3900,  and  3100  cal.  /mole 
for  I.  E.  C. 's  of  1.65,  1.  80  and  2.  15  meq. /g.  1,  respectively. 

These  values  are  close  to  the  activation  energy  for  viscous  flow  of 
water,  which  is  3700  cal /mole,  (33) 

The  above  considerations  show  that  the  concept  of  tortuosity  is  a 
valid  one  in  dealing  with  these  systems  and  that  Eq.  (II-8)  is  useful 
in  correlating  membrane  transport  phenomena. 
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4.  2.  5  Future  Work 

Data  on  electro-osmotic  water  transport  is  insufficient  to  report 
at  this  time.  Work  on  this  part  of  the  program  will  be  accelerated 
soon.  In  addition,  vapor  permeabilities  and  hydrodynamic  perme¬ 
abilities  of  other  membrane  systems  will  be  determined. 
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4.  3  TASK  III 

KINETIC  STUDY  OF  THE  ELECTROCHEMICAL  OXYGEN  (AIR)  REDUCTION 

H.  J.R.  Maget 


4.  3.  1  Introduction 

Rates  of  cathodic  oxygen  reduction,  that  is  cell  currents,  are  dependent 
on  partial  pressures  of  the  oxidant  if  rate-controlling  steps  involve 
the  concentration  or  partial  pressure  of  oxygen.  This  is  likely  to 
be  observed  since  reaction  rates  will  be  either  liquid  film  or  gas 
diffusion  controlled.  However,  cases  can  arise  where  removal  of 
the  reaction  product  may  be  hindered  by  slow  transport  processes, 
thus  resulting  in,  possibly,  appreciable  lower  rates.  This  could 
be  the  case  of  water  removal  from  the  catalyst  surface  of  an  oxygen 
electrode.  If  local  current  densities  are  either  dependent  on  partial 
pressures  of  oxygen  or  water,  it  will  become  necessary  to  establish 
relationship  predicting  such  local  current  densities  and  to  design 
electrode  geometries  favorable  to  uniform  current  distribution,  and 
as  a  result  uniform  distribution  of  the  main  influential  variables 
affecting  oxygen  electrode  performance. 

The  ultimate  goal,  that  is,  quantitative  description  of  current- 
voltage  relationships  for  an  IEM/electrode  structure  as  a  function 
of  main  variables,  can  be  attained  from  experimental  studies, 
including  at  least: 

-  the  determination  of  limiting  current  densities 

-  an  attempt  to  establish  the  rate-controlling  process  for  known 
electrode  structures  in  the  current  density  range  corresponding 
to  0.  85  -  0.  5  volts 

-  the  derivation  of  relationships  describing  the  current-voltage 
behavior  over  practical  operational  ranges 

Experimental  results,  as  well  as  the  equations  derived  to  predict 
limiting  local  current  densities  for  purely  gas-diffusion  controlled 
processes  have  been  obtained  during  this  period  and  will  be  discussed 
in  the  following  three  sections. 
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4,  3,  2  Influence  of  oxygen  partial  pressure  on  limiting  current 

densities  of  the  oxygen/IEM  electrode-electrolyte  system. 

4.  3.2.  1  Introduction 

In  establishing  over-all  oxygen  electrode  capabilities  for  fuel  cell 
application,  one  of  the  important  variables  is  the  limiting  current 
density.  Although  these  limiting  currents  are  not  simply  dependent 
on  some  limited  parameters,  i.e.  electrode  activity,  electrode 
structure,  electrolyte  properties,  local  temperatures,  anisotropic 
current  density  distribution,  etc  ....  actual  measurements  are 
valuable  if  some  reproducible  characteristics  can  be  controlled, 
i.e.  catalytic  activity,  electrode  structure,  electrolyte  properties. 
Thus,  measurements  are  valuable  for  a  specifically  designed  system. 
Such  a  system  has  been  chose  to  be  a  Pt  black  electrode  associated 
with  a  cationic  exchange  membrane  as  electrolyte.  Measurements 
have  been  restricted  essentially  to  room  temperature  operation, 
and  to  some  specific  humidity  conditions. 

4  3  2.  2  Experimental  Equipment  and  Procedures 

Partial  as  well  as  limiting  currents  were  obtained  by  applying 
manually  set  voltages  by  means  of  300  and  2000  ma  potentio stats. 

A  Luggin  capillary /SCE  system  is  used  as  the  reference  electrode, 
the  capillary  tip  being  tightly  held  against  the  membrane.  Such 
reference  electrodes  have  been  described  and  discussed  elsewhere,  (1). 
Due  to  the  limited  current  output  of  the  initially  available  potentiostat, 
and  to  expected  high  limiting  currents  at  high  partial  pressures  of 
oxygen,  the  circular  wet-proofed  electrodes  had  apparent  areas  of 
either  0.33  or  0.70  cm2  (Figure  III- 1).  These  electrodes  could  be 
operated  for  forced  flow  conditions  (Figure  III  1  B-l)  or  as  self- 
breathine  device  (Figure  III  1  B-2).  Larger  wet-proofed  electrodes 
(3.88  cnr)  were  operated  as  self  breathing  device  in  an  enclosed 
environment  in  absence  of  forced  flow  and  currents  were  measured  by 
means  of  an  2000  mA  Wenking  potentiostat.  Currents  were  recorded 
for  short  periods  of  time  on  a  Leeds  and  Northrup  Co.  Speedomax 
G  Recorder  at  chart  speeds  of  about  30  cm/min.  It  was  hoped  that 
surface  conditions  regarding  water  partial  pressure  and  accumulation 
would  be  unchanged  during  such  short  time  intervals.  Conditions 
surrounding  the  electrode  were  kept  constant  by  placing  the  device 
into  a  constant  temperature  and  humidity  chamber  (Aminco  Climate- 
Lab).  Electrode  surface  temperatures  were  measured  by  a  fine 
Pt-Pt/Rh  thermocouple,  inserted  into  the  electrode  catalyst  layer. 
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Complete  polarization  curves  could  be  obtained  by  this  procedure 
and  steady-state  currents  evaluated  by  extrapolation  of  current 
time  curves  to  t  =  0.  Steady  or  quasi-steady  state  currents  could 
be  obtained  approximately  15-25  seconds  after  the  applied  potential. 


4,  3.  2.  3  Experimental  Results 
Current  Decay 

All  experimental  results  obtained  by  applying  fixed  potentials  displayed 
initially  large  "surge"  currents.  However,  currents  would  decay 
rather  rapidly,  to  steady  state  or  quasi-steady  state  values.  Steady 
state  values  reported  in  the  following  figures  represent  extrapolated 
values  to  t  =  0  of  the  linear  current-time  functions.  As  expected, 
currents  would  increase  with  increased  applied  potentials  up  to 
certain  limiting  currents. 

Attempts  to  correlate  experimental  results  on  the  basis  of  classical 
equations  for  diffusional  processes  have  not  been  successful.  How¬ 
ever,  as  represented  in  Figure  III-2,  rates  of  current  decay  seem 
to  be  dependent  on  electrode  size,  suggesting  possibly  contribution 
of  radial  diffusional  processes  over  the  electrode  surface  or  in  other 
words  possibly  variable  diffusion  layer  thicknesses  dependent  on 
electrode  geometry.  An  attempt  to  obtain  analytical  solution  for 
purely  gas  diffusional  processes  will  be  discussed  later.  Current 
decays  for  open-end  tube  as  represented  in  Figure  III-2  suggest  an 
induction  period  before  steady  state  conditions  are  observed.  This 
result  could  be  obtained  if  the  open-end  tube  creates  initially  a  barrier 
to  high  diffusion  rates  followed  by  establishment  of  steady  state  conditions. 

Electrode  Polarization  and  Limiting  Currents. 

A.  Forced  Gas  Flow  Conditions 

Limited  experimental  data  are  available  for  large  flow  rates  over 
electrodes  according  to  configuration  IB-1  (Figure  HI-1).  The®e 
results  indicate  influence  of  gas  flow  as  compared  to  self-breathing 
device.  Such  differences  are  associated  with  diffusional  processes 
at  the  electrode  surface.  This  phenomenon  will  be  discussed  in 
greater  detail  in  the  study  of  channel  geometry  influences.  Limiting 
currents  for  air ,  50  O2/50N2  and  oxygen  at  27  C  and  1  atm.  total 
pressure  can  be  determined  from  Figure  III-3.  The  following  results 

were  obtained: 
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Partial  Pressure  of  Oxygen 
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These  results  allow  the  following  observations:  1)  up  to  about 
70  mA/cm^,  the  polarization  curves  behave  almost  identically 
and  independently  of  the  oxidant  partial  pressure  and  2)  limiting 
current  densities  are  nearly  proportional  to  oxygen  partial 
pressures. 

B.  Self-Breathing  Electrodes: 

Since  gas  flow  conditions  are  doubtlessly  very  influential  on 
limiting  current  densities,  most  experimental  results  have 
been  associated  with  non-flow  conditions,  on  so-called  self¬ 
breathing  electrodes,  configuration  IB- 2  (Figure  III  - 1 ) . 

Results  for  air  electrode  polarization  are  reported  in  Figure 
I1I-4,  where  initially  appreciable  polarization  can  be  observed 
only  for  low  oxygen  partial  pressure  and  where  air  operated 
electrode  under  flow  conditions  yield  results  comparable  to 
higher  oxygen  partial  pressures  for  current  densities  up  to 
about  70  ma/cm2,  Figure  III- 5  represents  quasi-linearity 
between  limiting  current  densities  and  oxygen  partial  pressures. 
These  results  were  obtained  for  three  electrode  sizes:  0.  33, 

0.  70  and  3.  88  cm2  corresponding  to  electrode  radii  of  0,  33, 

0.47  and  1.  1  cm,  respectively. 

Recent  results  reported  by  M.  B.  Clark,  in  the  literature  (34) 
on  the  influence  of  oxygen  partial  pressure  on  limiting  current 
densities  for  carbon  electrodes,  indicate  large  deviation  from 
linear  relationships  between  i^  and  P02'  these  alkaline  electro¬ 
lytes  limiting  current  densities  are  about  2  times  larger  than  for 
I  EM,  viz.  1200  ma/cm2  against  550  ma/cm2  for  pure  oxygen. 
However,  limiting  current  densities  are  quasi-linear  with  p q2 
up  to  about  0,61  atm.  and  are  identical  to  results  obtained  for 
the  Pt  black  -  IEM  system.  Since  the  compared  systems  are 
quite  different,  similitude  in  experimental  results  suggest  a  gas- 
diffusion  controlled  process,  and  deviation  from  such  a  mechanism 
in  the  case  of  pure  oxygen. 
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Influence  of  Electrode  Temperature 

Temperature  effects  have  not  been  systematically  investigated. 

However,  for  all  measurements,  fine-wire  Pt/Pt-Rh  thermo¬ 
couples  have  been  placed  in  the  Pt-black  electrodes  and  temperatures 
have  been  recorded.  Qualitatively,  temperature  effects  are 
appreciable,  and  limiting  currents  can  double  for  temperature 
changes  from  about  25  to  60°C.  Practically,  this  may  suggest 
fuel  cell  operation  under  conditions  such  that  heat  transfer  would 
allow  electrode  surface  temperature  elevations  up  to  conditions 
compatible  with  membrane  stability. 

4.  3.  2.4  Influence  of  Total  Pressure  on  Rates  of  Oxygen  Reduction 

Since  total  gas  pressure  is  expected  to  affect  limiting  currents  in  a 
gas-diffusion  controlled  process  (influence  on  the  diffusion  co¬ 
efficient,  partial  pressure  gradient)  as  well  as  in  a  liquid  film 
controlling  process  (influence  on  the  gas  solubility,  partial  pressure 
gradients)  experimental  results  were  obtained  to  describe  such 
affects  with  air  at  total  pressures  up  to  about  8  atm. 

4.  3.  2.  4.  1  Experimental  Equipment  and  Procedures. 

The  experimental  apparatus  represented  schematically  in  Figure  III-6, 
was  designed  for  operation  up  to  10  atm.  total  gas  pressure.  Gases 
were  allowed  to  enter  the  two  compartments  at  a  rate  such  that 
pressure  differences  across  the  membranes  were  always  less  than 
0.  3  atm  during  pressurizing.  Pressures  could  be  read  within  0.01 
atm.  By  proper  manipulation  of  exhaust  valves,  air  was  allowed  to 
flow  over  the  electrode  at  the  operating  pressure.  A  Luggin  Capillary/ 
SCE  system  was  used  as  reference  electrode.  Such  reference  electrodes 
have  been  described  elsewhere  (1).  The  circular  wet-proofed  electrodes 
had  an  apparent  area  of  0.  33  cm2.  Time -dependent  currents  were 
recorded  at  chart  speeds  up  to  30  cm/min.,  for  various  potentials 
applied  manually  by  means  of  a  2000  mA  Wenking  potentiostat. 

4.  3.  2.4.  2  Experimental  Results 
Open  Circuit  Voltage 

Open  circuit  voltages  recorded  on  the  same  complete  fuel  cell  as  a 
function  of  total  pressure  and  for  two  gas  composition  are  represented 
in  Figure  III-7.  A  logarithmic  function  of  total  pressure  vs.  OC/V 
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yields  linear  relations  for  both  systems,  suggesting  the  following 
pressure  dependence  in  the  case  of  air: 

E  =  Eo  +  0.056  log  P  (III-l) 

where  E  =  OC/V 

Eo  =  OC/V  at  P  =  1 

P  =  gas  pressure  in  atm.  at  the  air  cathode  and 
HJ  anode 

Limiting  Current  Densities 

On  Ag  in  acidic  and  neutral  media,  Krasilschchikov  (35)  has  shown 
the  relationship: 

i  =  ->.r ‘*rr  £  (in- 2) 


(where  4.  =  0.  5)  to  hold  between  0.  21-70  atm.  This  result  was 
supported  by  Winkelmann  (36)  for  bright  Pt.  M.  B.  Clark  (34) 
reports  linear  relationships  between  i^,  and  P  #  (inverse  of  total 
pressure).  Results  obtained  on  IEM/Pt  black  systems  have  yielded 
linear  relationship  between  ij^  and  \fp".  The  slope  for  21%  02-79%  N2, 
Figure  III -8,  correspond  to  0.011  Amp-atm'1'  .  In  order  to  show 
possible  effects  of  bulk  diffusion  coefficients,  a  binary  gas  system 
was  chosen  displaying  much  larger  diffusion  coefficient,  i.  e.  0.  75 
cm2/sec.  at  27°C  for  He-02  (37)  as  compared  to  0.25  cm2/sec.  (38) 
for  N2-O2  mixtures.  Since  diffusion  coefficients  are  inversely  propor¬ 
tional  to  total  pressure  and  nearly  independent  on  gas  composition,  the 
slope  of  0.030  Amp. -atm'1'2  for  02-He  could  be  expected  as  compared 
to  air.  Although  general  trends  in  variations  of  limiting  currents  can 
be  explained,  it  remains  mandatory  to  obtain  quantitative  interpretations. 

4.  3.2.  5  Discussion 

Equations  derived  in  Appendix  III- 1  for  rate-controlled  diffusional 
processes  in  the  gas  bulk,  i.e. 

i  =  47 mFDbc0  tfn-22> 
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for  gas-controlled  diffusion  in  the  pores 
I  (pore)  =  4 


.iusion  in  uio 

irnFC^bc0^- 

m*l 


cosh  1 

KmSinh(Xmh/b) 


(III -2  7) 


and  for  liquid  phase  controlled  diffusion: 


(III- 35) 


yield  results  too  large  to  justify  a  process  described  by  Eqn.  Ul-22. 
However,  rate-controlled  processes  in  the  pores  could  be  either  gas 
diffusion- controlled  (Eqn.  III-27)  and  will  depend  essentially  onJ>  e 
and  b,  if  (h/b)  is  large,  or  liquid  film  controlled,  to  become  dependent 
onfi  e  *  ^  and  the  gas  solubility  CQ  in  the  thin  film.  It  is  not 
possible,  presently,  to  make  a  definite  distinction  between  these  two 
transport  mechanisms  in  order  to  explain  experimental  results.  How¬ 
ever,  it  will  be  shown  in  the  case  of  oxygen  (air)  reduction  in  restricted 
channels  that  bulk  diffusional  processes  may  become  small  enough  to 
justify  observed  local  limiting  current  densities  as  compared  to  more 
complex  pore-diffusion  controlled  transport  rates. 

4.  3.  3  Self-Breathing  Air  Electrodes 

4.  3.  3.  1  Introduction: 

Since  it  may  be  considered  as  a  fact  that  ". . .  fuel  cells  must  be  made^ 
to  operate  on  air  before  commercial  applications  can  be  contemplated", 
(39)  and  since  it  has  been  shown  previously  that  limiting  current  are 
essentially  proportional  to  oxygen  partial  pressure  in  the  oxidant 
stream,  it  becomes  imperative  to  measure  and  possibly  understand 
polarization  losses  incured  in  practical  air  electrodes.  Such  losses 
could  be  minimized  by  compressing  the  air.  However,  compressors 
could  only  be  justified  for  "fuel  cells  with  power  output  of  the  order 
of  100-200KW"  (39).  For  smaller  unit  to  be  operated  at  atmospheric 
pressure,  careful  air  channel  and  air  electrode  design  will  be  required 
in  order  to  minimize  diffusion  polarization. 

The  purpose  of  the  present  work  was  to  obtain  experimental  results  and 
interpretation  to  explain  polarization  in  channels  of  defined  geometry 
and  to  establish  influential  parameters  which  would  affect  electrode 
performance.  In  order  to  establish  applicability  of  a  self-breathing 
air  electrode,  viz.  without  forced  convective  air-flow,  for  low 
current  densities,  experimental  investigations  were  started  on  straight 
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air  channels.  Furthermore,  since  such  a  self-breathing  electrode 
had  to  operate  away  from  limiting  currents,  in  order  to  minimize 
diffusion  polarization,  additional  work  for  forced  flow  (at  various 
air  flow  rates)  was  conducted  in  straight  air  channels. 

4.  3.  3.  2  Experimental  Equipment 

Partial  as  well  as  limiting  currents  were  obtained  by  applying 
manually  set  voltages  by  means  of  a  300  Ma  Wenking  potentiostat. 

A  Luggin  Capillary/SCE  system  is  used  as  the  reference  electrode, 
the  capillary  tip  being  tightly  held  against  the  membrane.  Such 
reference  electrodes  have  been  discussed  and  described  elsewhere(  1). 

Additional  investigations  were  undertaken  for  galvanostatic  measure¬ 
ments.  The  equipment  is  represented  in  Figure  III -9.  Identical 
reference  electrodes  were  used. 

All  experimental  work  was  conducted  on  air  channels  2-1/2"  long  and 
1/2"  wide.  Ch?~.,el  height  could  be  caried  by  changing  a  removable 
Lucite  bar  placed  on  the  channel  top.  Thus,  channel  heights  could 
be  1/16,  1/8,  1/4  and  1/2".  The  channel  was  mounted  on  an  air 
electrode  and  placed  in  a  large  Lucite  box  in  order  to  avoid  small 
air  flow  sweeps  over  the  air  electrode. 

In  order  to  determine  local  current  densities,  the  air  electrode  was 
manufactured  by  placing  ten  parallel  electrode/screen  strips  on  a 
GPS-Ion  Exchange  Membrane.  (Electrode  position  is  described  in 
Figure  III  - 9. )  Gaps  of  1/16"  between  electrodes  allowed  for  electrical 
insulation  of  the  various  electrodes.  Electrode  dimensions  were 
1/2"  x  3/ 16"  with  an  actual  area  of  0.  60  -  0.05  cm  .  The  IEM  exten  e 
out  of  the  channel  for  reference  potential  measurements  and  stainless 
rods  contacted  the  screens  for  current  pick-up. 

The  counter  electrode  (H2-electrode)  was  prepared  in  a  similar  manner 
Catalytic  electrodes  faced  each  other  across  the  electrolyte  (membrane). 

The  two  ends  of  the  self-breathing  channel  were  open,  thus  displaying 
planar  symmetry  on  either  side  of  the  channel  center  cross-section. 
Closing  one  channel  end  actually  corresponded  to  doubling  the  channel 

length. 
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4.  3.  3.  3  Experimental  Results 

{ 

1.  Individual  Electrode  Polarization 

— , - 

Experimental  electrode  polarizations  are  presented  in  Figures  III- 10 
and  III- 11  for  various  channel  heights.  The  following  general  obser¬ 
vation  regarding  these  figures  can  be  made: 

a.  polarization  characteristics  are  almost  identical  for  all 
electrodes  for  current  densities  up  to  about  25  ma/cm2  for  a 
channel  height  of  1/8",  although  the  initial  slope  of  the  I-V  curves 
are  not  all  identical. 

b.  limiting  currents  are  larger  (in  essence  the  limiting 
current  for  exposed  electrodes)  for  terminal  electrodes  (1  and 
10)  than  other  electrodes  and  limiting  currents  decrease  harmon¬ 
iously  from  terminal  to  central  electrodes. 

c.  central  electrodes  5  and  6  display  minima  in  current  den¬ 
sities  as  expected  from  consideration  of  central  planar  symmetry. 

d.  channel  height  affects  performance  to  a  large  degree,  i.  e. 
limiting  currents  for  inner  electrodes  are  22mA  and  35mA  for 
channel  heights  of  1/16"  and  1/8"  respectively.  These  current 
densities,  i.e.  37  and  60  Ma/cm2,  correspond  to  oxygen  partial 
pressures  of  about  0,  07  and  0.  1,  respectively.  No  appreciable 
polarization  is  observable  for  channel  height^  0.64  cm  (1/4"). 

e.  for  identical  polarization  potentials  (as  determined  for 
polarization  curves)  current  density  distributions  display  minima 

for  inner  electrodes  and  can  be  extrapolated  to  edge  values  of  limiting 
currents  as  determined  from  previous  measurements.  Strictly  speaking, 
these  current  densities  cannot  be  described  from  Figures  III- 10  and 
III- 11.  However,  the  partial  pressure  of  oxygen  in  the  channel  would 
not  be  too  greatly  affected  at  higher  current  corresponding  to  limit 
currents  of  the  edge  electrodes.  Results  would  be  affected  greatly 
at  the  edge  electrodes,  let's  say  electrodes  1  and  2  and  their  corres¬ 
ponding  symmetrical  position.  Current  density  distributions  for 
channel  heights  of  1/8  and  1/16"  for  various  applied  potentials  are 
represented  in  Figures  III- 12  and  III- 13.  Edge  currents  (about  80mA) 
can  become  2  to  4  times  larger  than  center  channel  currents,  but 
display  much  more  uniform  distribution  for  lower  electrode  polari¬ 
zation,  i.e.  0.5  volts  vs.  SCE  for  1/8"  channel. 
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Interpretation  of  Individual  Polarization  Curves 


Current  potential  behavior  for  the  investigated  systems  can  be 
represented  by  the  following  Figures: 

diffusion  limited  current 


/'duf fusion 
POLARIZATION 


POTENTIAL-  VOLTS 

FIGURE  A 


/  rl(jmEo~  ^cjce 
/  Ea.m  actual  voltage 

Eujce  m  voltage  in  absence  of 
cbonne7  effects 

POTENTIAL  ~  VOLTS 

FIGURE  B 


Figure  A  represents  polarization  characteristics  as  observed  in 
absence  of  diffusion  polarization  (probably  activation  polarization, 
in  absence  of  ohmic  contribution,  which  are  generally  eliminated). 
Figure  B  represents  strictly  diffusion  polarization  terms  in  absence 
of  other  possible  polarization. 

Polarization  characteristics  can  be  represented  by  the  equation: 


O'  ■*)  =e‘ 


(E°  -  E) 


(IU- 3) 


where  i  and  ij^  represent  actual  and  limiting  currents,  respectively, 
«^=  f  nF/RT,  E°  =  reference  potential  in  absence  of  channel  polari¬ 
zation,  E°  -  E  =1|j 


Ln(  1  -  i/iL)  =  4.  (E°-E)  =«^d 


By  differentiation: 


di  ~<ULU  -  i/1  l.) 


(Ill -4) 


(III -5) 
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(III- 6) 


Since  i,  (x)  displays  a  monotonous  decrease  up  to  channel  center, 
equation  III -6  is  expected  to  display  in  increase  from  channel 
edge  to  center.  Figure  III- 14  represents  smoothed  experimental 
results  for  1/4,  1/8  and  1/16"  channel  heights.  Since  no  additional 
position-dependent  polarization  was  observed  for  ye  =  1/4", 
the  value  of  Rd  =  4  was  chosen  as  reference  systems  resistance  to 
obtain  a  relationship  between  channel-induced  diffusion  resistance 
Rd(x)  and  ij^(x).  Now,  in  equation  III-6,  the  slope  (dE/di)  defined 
as  diffusional  resistance  becomes  dependent  of  and  ii_.  (If  iL 
can  be  determined  analytically  and  since  <  represents  defined 
constants,  (dE/di)  is  defined. )  Equation  III- 6,  if  represented  by 
/  /AM  should  yield  a  constant  value,  determining 

exp. 

AV/AI  representing  the  diffusional  resistance  should  display: 


consta 


a.  no  diffusional  resistance  at  the  edges,  at  least  as 
related  to  electrode  geometry 

b.  maximum  diffusional  resistance  for  center  electrodes 


c.  increased  resistance  for  reduced  channel  cross-section. 

Experimental  results  regarding  these  observations  are  represented 
in  figure  III-15  for  channel  heights  of  1/8  and  1/16",  and  different 
electrode  positions. 

Table  III- 1  represents  results  obtained  for  1/8"  and  1/16"  channel 
heights. 
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TABLE  HI-1 


n* termination  of  from  experimental  data 


Electrode 

Limiting 

1 

LM.) 

Ui_ 

-rf  1 

Position 

Current,  iL 

*L 

(Amps)- * 

UV  exp. 

U1 

/  -R 
(Volt)*1 

Channel  Height 

Amps 

1/16" 

(ohm) 

(ohm) 

1 

0.041 

- 

• 

2 

- 

- 

- 

- 

- 

3 

0.028 

35.8 

7.3 

3.  3 

10.  8 

4 

0.024 

41.7 

7.8 

3.8 

11.0 

5 

0.022 

45.5 

8.2 

4.2 

10.8 

6 

0.025 

40.0 

8.  1 

4.  1 

9.8 

7 

X).  027 

38.0 

7.9 

3.9 

9.7 

8 

0.030 

33.3 

7.3 

3.3 

10.  1 

9 

10 

Channel  Height 

0.050 

1/8" 

20.0 

6.1 

2.  1 

9.6 

Average  *-= 

1 

- 

- 

- 

• 

2 

0.045 

22.2 

6.0 

2.0 

11.  1 

3 

0.043 

23.2 

6.7 

2.7 

8.6 

4 

0.037 

27.0 

7.  1 

3.  1 

9.0 

5 

0.035 

28.6 

7.3 

3.3 

8.7 

6 

0.035 

28.6 

7.3 

3.3 

8.7 

7 

0.037 

27.0 

7.  1 

3.  1 

9.0 

8 

0.  044 

22.7 

6.8 

2.8 

8.  2 

9 

- 

- 

- 

- 

- 

10 

- 

- 

- 

- 

Average  = 
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Evaluation  of  the  Limiting  Current. 

The  analytical  evaluation  of  ij^  (K)  will  allow  to  describe  the  polari¬ 
zation  curves  means  of  Equation  III-6.  Local  limiting  currents 
can  be  evaluated  by  solving  Laplace's  equation,  providing  transport 
by  convection  is  considered  to  be  negligible.  Detailed  calculations 
are  presented  in  Appendix  IU-2. 

The  concentration  distribution  for  oxygen  and  water  in  the  channel 
are  represented  by: 

_  OO  -  fj 


/>*/ 

where  Cn  and  C  represent  local  concentration  of  oxygen  and  water. 
w 

C°  =  oxygen  partial  pressure  at  channel  inlets. 

g(T)  and  f(T)  represent  partial  pressures  of  water  under  equilibrium 
conditions  at  y  =  yQ  and  y  =  o,  respectively. 
yQ  =  channel  height 

(-/)n- 1  j  tjjr 

^  w 

Local  current  densities  can  be  evaluated  from: 

and  total  channel  current  from:  ^ 

Jtydx  (m'101 

I  =  total  channel  current 
J  (x)  =  local  current  density 
W  =  channel  width 

Equations  can  be  derived  for  diffusion  of  oxygen  or  water  vapor  through 
the  ternary  gas  system. 


(IH-$) 


(III- 10) 
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Calculated  local  current  densities,  thus  become: 

j(x)o=  ianf)(X„y6 


and: 


aw, 3.. 


j(x\t  = 


zn^nxhCr) 


2 


Smfrn  y 


where: 


'i-cosha„y0 

smh»ny0 


-] 


(III- 11) 


(III- 12) 


n'  =  number  of  equivalent/mole  (n'  =  4  for  and  n'  =  2  for  H^O) 
p  =  Faraday  constant  =  96,  500 

D  =  Diffusion  constant  of  species  i  in  the  gas  medium 
2b  =  Channel  length  (6.  35cm) 

y  =  Channel  height  0.  32  and  0.16  cm 


Equation  III- 12  has  been  obtainel'  for  g(T)  =f(T)  =h(T) 
Application  to  specific  environmental  conditions 
For  operation  with  air,  equation  III- 11  becomes: 


j(*)0 ■  0.2$S^sm  (nr&)  hnk  (wrj |) 

1  n-/,3 

,,  r-  /  x\  r  i-coshu.^ 
0.03/' )sm  (nr a)  ^ 

n=i,i  L  J 

Equation  III- 13  and  III- 14  are  based  on 


3)  0z  =  o.2i  cm2/sec 
j)  w  =  o.  29  cm2/se c 

h(T)  =  30mmHg  (vapor  pressure  of  water  at  30°C) 


(III- 13) 


(III- 14) 


Some  calculated  data  are  presented  in  Table  III  -  2  for  a  channel 
height  of  1/16"  =  0.16cm. 
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TABLE  III -2 

Local  Current  Density  for  1/16" 

Channel  Height 

Minimum 
Current  Density 
j(x)  min. 

Position 

Local  Current  Density  j(x)/j(x)  min. 
j (x)  in  mA/cm^ 

Oxygen  Diffusion  b 

117 

1.0 

b/2 

163 

1.  41 

117 

b/4 

298 

2.  5C 

Water  Diffusion  b 

20 

1.0 

b/2 

22 

1.  1 

20 

b/4 

35 

1.  7 

Measured  b 

38 

1.0 

b/2 

44 

1.  2 

38 

b/4 

61 

1.  6 

0 

110 

2.9 

Calculated  and  experimental  results  are  presented  in  Figure  III-16 
for  the  1/16"  channel  and  for  fixed  diffusion  coefficients.  Data  can 
be  represented  fairly  well  by  assuming  a  water  diffusion-controlled 
process. 
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4.  3.  3.  4  Discussion 


4.3.4 


Restriction  of  channel  cross-section,  i.  e.  of  yQ,  may  result  in 
appreciable  reduction  in  gas-phase  transport  rates.  In  fact, 
these  rates  may  be  calculated  from  purely  gas-diffusion  controlled 
processes.  In  channel  where  forced  air  flow  rates  may  be  small, 
unfavorable  current  density  distributions  may  also  be  encountered. 

Such  conditions  and  results  will  be  discussed  in  section  4.3.4. 

The  results  also  suggest  reduction  of  channel  length  (for  defined  yQ) 
in  an  attempt  to  maintain  quasi-uniform  oxygen  partial  pressure 
and  surface  temperature.  The  results  indicate  applicability  of  self¬ 
breathing  electrodes  up  to  about  25  mA/cm  (at  ambient  temp.) 
in  a  design  where  air-electrodes  would  be  disposed  face-to-face  and 
separated  by  air-channels  such  that  yQ>  1/4".  Better  current  distri¬ 
bution  may  also  be  expected  for  disc- shaped  electrodes,  by  careful 
selection  of  the  channel-gap/electrode-radius  ratio.  The  stability 
of  such  a  design  would  require  additional  investigations  related  to 
heat  and  mass  transfer  processes. 

Current  distribution  in  air-breathing  electrode  channels  under  forced  flow 
conditions. 


4.  3.  4.  1  Introduction 

Experimental  results  and  interpretation  presented  under  the  title: 
"Self-Breathing  Air  Electrodes"  apply  strictly  to  non-flow  conditions 
and  represent  limiting  conditions,  since  fuel  cells  operating  under 
high  current  drains,  will  be  subjected  to  air  flowing  over  the  electrodes 
Other  limiting  conditions  would  correspond  to  infinite  flow  rates  o 
air,  insuring  identical  oxygen  partial  pressure  over  the  electrode. 
However,  such  operating  conditions  can  not  be  realized,  although 
approached,  under  practical  flow  conditions.  In  fact,  in  the  case  of 
the  IEM  air-operated  cathode,  air  flow  rates  have  to  be  adjusted  in 
order  to  remove  product  water,  and  only  product  water,  lest  the 
membrane  water  would  decrease,  generating  non-uniform  surface 
and  bulk  conditions,  very  detrimental  to  proper  operation  of  the 
complete  system. 

In  order  to  measure,  analyze  and  possibly  interpret  effects  of  air 
flow  rates  on  local  current  densities  in  air-breathing  channels, 
additional  work  was  conducted  on  the  channels  discussed  earlier. 
Air-flow  rates  were  maintained  at  low  Reynolds  number,  i.  e.  in  a 
cases  less  than  20. 
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4.  3. 4.  2  Experimental  Equipment 

Partial  and  limiting  currents  were  measured  indirectly  by  galvano- 
static  measurements.  A  Luggin  Capillary/SCE  system  was  used 
as  reference  electrode,  the  capillary  tip  being  tightly  held  against 
the  membrane.  The  experimental  equipment  is  represented  in 
figures  UI-17  and  III- 18.  Air  channels  are  2-1/2"  long  and  1/2"  wide. 
Channel  height  can  be  varied  by  moving  a  Lucite  bar  placed  on  the 
channel  top.  Channel  heights  could  be  1/16,  1/8,  1/4  and  1/2". 

The  channel  was  mounted  on  an  air  electrode.  In  order  to  avoid 
appreciable  entrance  effects  due  to  possible  turbulence,  the  air 
velocity  profile  was  generated  in  a  channel  similar  to  the  reaction 
channel,  and  had  a  total  length  of  3".  (See  figure  III- 17).  The  air 
electrode  was  prepared  by  placing  ten  parallel  electrode/screen 
strips  on  a  GPS  -  Ion  Exchange  Membrane.  The  IEM  extended  out 
of  the  channel  for  reference  potential  measurements.  The  counter¬ 
electrode  was  prepared  in  a  similar  manner.  One  channel  end 
(exhaust  C^/^mixture)  was  closed.  This  may  affect  currents  at 
channel  end  because  of  local  turbulence.  In  the  case  of  open  exhaust 
channel,  -'nd-effects  were  observed,  corresponding  to  oxygen  back 
diffusion,  thus  to  higher  local  currents.  Flow  rates  for  in-conditions 
were  measured  with  calibrated  flow-meters  and  for  small  flow  rates 
with  a  soap-bubble  displacement  technique. 

The  assembled  experimental  set-up  is  represented  schematically  in 
figure  III-  18. 

4.  3.  4.  3  Experimental  Procedures: 

Under  galvanostatic  conditions  (for  currents  limited  to  700  ma, 
corresponding  to  current  densities  of  116  ma/cm  ),  and  for  steady-state 
operation,  local  cell  potentials  were  measured  for  various  flow  condi¬ 
tions.  In  order  to  show  effect  of  air  flow  and  to  account  for  convective 
processes,  similar  results  were  obtained  in  absence  of  forced  flow. 

The  presence  of  the  pre-channel  reduces  gas  diffusion  to  electrode  #1, 
allowing  for  small  local  currents,  whereas  if  reaction  channel  end  is 
open  to  the  atmosphere,  back-diffusion  of  oxygen  is  occurring,  even 
at  low  air  flow  rates,  thus  displaying  higher  local  current  densities 
in  position  10  as  compared  to  electrode  #1. 

4.  3.  4.  4  Experimental  Results: 
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Polarization  Characteristics 

Polarization  characteristics  for  individual  electrodes  in  the  reaction 
channel  are  represented  in  figure  III- 19  for  individual  cases  as 
related  to  electrode  position,  air  flow  rates  and  channel  geometry. 

On  the  basis  of  limiting  currents  observed  for  flow  conditions  over 
electrodes,  total  channel  limiting  current  should  be  about  800  mA, 
corresponding  to  16  cc7min.  of  air  flowing  through  the  channel. 

Current  distributions  are  reported  for  flow  rates  from  1  to  60  cc/min. , 
or  linear  velocities  of  0.08  to  1.0  cm/sec.  for  1/16"  channel  height. 

In  figure  III- 19  it  appears  that  appreciable  stoichiometric  reaction 
air  flow  rates  are  required  to  avoid  limiting  currents  smaller  than 
edge  limiting  currents  observed  on  electrode  #1  already  at  low  flow 
rates. 


Limiting  Currents: 

Figure  III- 20  represents  limiting  currents  (normalized  with  respect 
to  the  limiting  current  at  x  =  0)  as  a  function  of  position  and  air  inlet 
conditions.  These  results  were  obtained  from  polarization  characteris¬ 
tics  as  described  by  figure  III- 19  and  extrapolated  to  (iL^L'o)  =  1 
for  electrodes  near  air  inlet  channel. 


On  the  assumption  that  these  limiting  currents  can  be  derived  from 
gas  phase  diffusion  processes,  an  attempt  has  been  made  to  solve 
differential  equations  descriptive  of  the  process,  considering  defined 
boundary  conditions.  (For  details  of  the  derivation,  see  Appendix 
III-  3). 


On  the  basis  of  a  simplified  (approximately  correct  near  the  electrode 
surface)  linear  velocity  profile,  it  was  possible  to  derive  a  relation 
describing  limiting  currents  iL  as  a  function  of  position  (x  in  cm), 
gas  flow  rates  (Q  in  cm3  air/sec.),  channel  height  (yQ  in  cm),  gas 
bulk  concentration  (c0  in  moles/cm3),  channel  width  (W  =  const,  in  cm) 
and  diffusion  coefficient  (fi  in  cm  /sec.).  \y 

iLM  =  ZnFD^Co(^T)5 

If  Eqn.  Ill- 15  is  normalized  by  dividing  by  iL  (x  =  0)  -  130  ma/cm2, 

a  linear  relationship  should  be  obtained  for  the  log.  function,  i.  e. 


(Ill- 16) 
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Experimental  results  reported  in  figure  III-21  for  yQ  =  0.  16  cm,  air 
flow  rates  varying  from  1.0  to  6.9  cm^/sec,  and  to  electrode 
positions,  i.e.  #1,  2  and  5,  yield  a  linear  function  with  slope  equal 
to  2.  The  specific  extrapolated  conditions  (d/xy^0)  =  1,  yields 
log  A  =-1.0  from  figure  III-21  whereas  the  calculated  value  from 
eqn.  IU-  15  corresponds  to  log  A  =  0.  76. 

Presently,  experimental  results  suggest  the  application  of  a  semi- 
empirical  relationship  defined  as 

Description  of  the  polarization  curves 

As  in  Section  4.  3.  4.  2  (on  "Self-Breathing  Air  Electrodes",  regarding 
diffusion  polarization),  the  polarization  curves  can  be  completely 
described,  if  in  eqn.  IU-3, 


(-4)= 


(in -3) 


limiting  current  iL  (*)  can  be  analytically  derived  and  if  the  value 
of  4.  can  be  empirically  established  or  evaluated  from  its  physico¬ 
chemical  significance.  However,  in  this  case,  iL  W  is  not  only 
position-dependent,  but  also  dependent  on  air  flow  rates  and  gas 
properties.  In  principle,  all  experimental  results  for  various  flow 
conditions  and  channel  geometries  could  be  correlated  by  eqn.  III-3 
A  plot  of: 


(III -4) 


should  yield  linear  relationship  as  shown  in  figure  III-22,  with  slope 
^  =  -11  0.  Experimental  results  reported  in  figure  III- 22  are  related 
to  two  channel  geometries  (1/8"  and  1/16")  and  two  electrode  positions 
(#5  and  #8).  Data  scattering  is  not  unexpected,  since  diffusion  polari¬ 
zation  •>«  d  is  obtained  from  actual  observed  potentials  and  potentials 
observed  for  air  -  operated  cathode  in  absence  of  channel  effects, 
i.e.  very  high  air  flow  rates.  Voltage  scattering  can  also  occur 
because  of  initial  rest  potential  variation  and  dissimilitude  in  the 
various  electrode  segments  in  the  channel.  Significance  of  *  has  been 
discussed  in  Section  4.  3.  3.  3,  and  values  obtained  for  forced  flow 
conditions  are  comparable  to  *  =  10.  2  and  9.0  for  channel  heights  of 
1/16  and  1/8",  respectively,  in  the  case  of  self-breathing  electrodes. 
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In  view  of  the  previous  discussion,  the  following  semi-empirical 
lationship  can  be  used  to  describe  polarization  characteristics 


re 


in  channels 


i(*)  -  4.I03  iL(o)  1111 


where  i  =  Amp/cm  , 

i  (0)  =  limiting  current  density  (Amp/cm2)  in  absence  of  channel 

eVfects,  andd^n  =  0.28. 

In  addition  to  I-V  relationships,  approximate  equations  to  describe 
oxidant  concentration  distribution  can  be  used  in  relation  to  local 
current  densities: 

C0  7)FDc0  \nJ  Iffys) 


-18) 


III- 19) 


when  77*  =  0; 


•WALL 


(111-08) 


nFDCo  V  NJ  f(2/5) 


leading  to  the  limiting  current  expression  (Appendix  III-3),  when 


Cwall  =  °’  Ue’  i  =  iL 


4.  3.  4.  5  Conclusion: 


Experimental  results  and  analytical  result  interpretation  have  shown 
the  importance  of  channel  geometry  for  air-operated  current  collectors 
and  channels.  Appreciable  diffusion  polarization  and  poor  current 
density  distribution  can  only  be  avoided  for  high  air  flow  rates  i  long 
air  channels  are  to  be  used.  Designs  should  take  care  of  all  influen¬ 
tial  variables,  suggesting  large  channel  height  or  better  yet,  sma 
air  channel  length  to  obtain  good  pq2  distribution.  These  observations 
are  specially  true  for  operation  near  limiting  currents.  For  fractional 
parts  of  i.  ,  i.e.  0.  5  i^  flow  rates  *  2  times  stoichiometric  rates 
for  i.  are  sufficient  to  obtain  quasi-uniform  rate  distributions. 

i-i 

Air  electrode  can  operate  near  iL  for  very  high  flow  rates  for  short 
channels.  These  limiting  currents  can  then  be  increased  by  increasing 
electrode  surface  temperature  and/or  oxidant  total  pressure,  if  the 
partial  pressure  is  to  be  maintained  constant. 
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Over-all  Heat  and  Mass  Transport  processes  will,  of  course,  be 
affected  by  suggested  designs.  New  designs  are  discussed  in 
Task  IV,  where  experimental  results  regarding  modified  air 
current  collectors  are  presented. 

4.4  TASK  IV 

HEAT  AND  MASS  TRANSFER  INVESTIGATIONS  OF  THE 
AIR-BREATHING  I.  E.  M.  FUEL  CELL 

J.  P.  Dankese 


4.  4.  1 


Introduction: 

Substantially  all  of  the  necessary  data  for  a  heat  and  mass  transfer 
analysis  and  possible  correlation  are  available.  Since  the  heat 
transfer  data  reduction  and  analysis  are  not  completed,  the  results 
of  this  work  will  be  reported  in  the  Final  Report. 


Included  in  this  report  are  descriptions  and  photographs  of  the 
experimental  apparatus  and  fuel  cell  hardware;  the  performance 
characteristics  of  the  new  design  air  electrode  current  collector, 
the  results  of  a  mass  transfer  analysis  which  included  bulk  gas, 
pore,  and  gas-liquid  diffusion  mechanisms;  a  comparison  of  the 
water  transport  and  evaporation  rates  for  a  "GPS"  I.  E.  M.  fuel 
cell  membrane  and  pure  water;  an  empirical  equation  for  the  effect 
of  air  velocity  on  the  rate  of  water  transport  and  evaporation  for  a 
"GPS"  I.  E.  M.  fuel  cell  with  and  without  electrochemical  reaction; 
and. finally  possible  over-all  mass  transfer  correlating  relationships. 


4.  4.  2  Experimental  Apparatus  and  Procedure 
4.4.  2.1  A.  I.  E.  M.  Fuel  Cell  Hardware 


1.  Perforated  Metal  Screen 

Perforated  metal  screens  made  by  a  photo-etching  process  are  used. 
This  screen  may  be  made  from  any  metal  desired,  but  for  the  experi¬ 
mental  investigations  reported  herein,  the  materials  were  Carpenter 
20  Cb  stainless  steel  and  molybdenum.  The  practical  ranges  o  s 
thickness,  hole  diameter,  and  percentage  flow  area  for  fuel  ce  app  i- 
cations  appear  to  be  0.003  to  0.010  in.  and  5  to  50%,  respectively. 

The  range  of  the  number  of  holes  per  geometric  square  inch  of  screen 
area  is  then  about  160  to  18,  000.  The  screens  used  in  the  expenmenta 
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work  were  0.  005  in.  thick,  had  a  hole  diameter  of  0.  005  in. ,  and 
contained  from  2000  to  15,  000  holes  per  in2.  The  percent  flow 
areas  were  from  10%  to  35%. 

Representative  cross-sectional  profiles  of  perldfated  metal  sheets 
are  shown  in  Photograph  IV- 1.  The  top  profile  represents  the  35% 
flow  area  screen  with  0.  005  in.  holes  and  an  inlet  to  throat  diameter 
ratio  of  about  2.  5  to  1.  The  bottom  profile  represents  a  ratio  of 
3.5  to  1,  a  percent  flow  area  of  11%  for  a  0.  0075"  screen  thickness. 

The  holes  are  smooth,  without  burrs,  and  are  streamlined  with 
convergent-divergent,  convergent,  and  divergent  nozzle  type  shapes. 

2.  Perforated  Screen  Air  Electrode  Current  Collector 

A  close-up  photograph  of  the  new  design  perforated  metal  screen  air 
electrode  current  collector  is  shown  in  Photograph  IV-2.  This  side 
makes  contact  with  the  catalyst  by  means  of  the  perforated  metal 
flat  ribs  shown.  Air  enters  the  alternate  flat  rib  flow  passages 
perpendicular  to  the  long  dimension  of  the  collector,  passes  through 
the  rib  face  (45°angle),  across  the  catalyst  reaction  zone,  and 
reacted  air  plus  product  water  vapor  then  passes  through  the  opposite 
rib  face  and  out  the  alternate  ribs  (from  the  end  with  the  white  tabs). 

The  white  battery  separator  material  prevents  bridging  of  the  air 
flow  passages  with  product  water  and  aids  in  its  removal  from  the 
fuel  cell.  The  dark  silicone  material  is  used  to  form  the  desired  flow 
pattern.  A  typical  0.004  in.  insulated  Pt-Pt-10%  Rh  thermocouple 
is  shown  on  the  white  material.  The  collector  temperature  is 
measured  at  several  points  during  cell  operation. 

An  isometric  view  of  the  same  collector  showing  the  reaction  air  side 
up  and  the  cooling  air  side  down  is  Photograph  IV- 3.  The  alternating 
pattern  of  the  reaction  and  cooling  air  passages  is  clearly  shown. 

Two  collector  thermocouples  are  fixed  to  the  battery  separator  material 
with  thin  transparent  tape.  Another  view  is  shown  at  the  bottom  of 
Photograph  IV-4.  An  alternate  configuration  was  the  cooling  air  side 
ribs  perpendicular  to  the  reaction  air  side  ribs  instead  of  parallel 
and  permits  independent  air  flow  conditions  on  both  sides.  Thus  a 
cross-flow  heat  exchanger  is  incorporated  with  each  cell.  A  patent 
docket  has  been  written  for  the  heat  exchanger  use  and  copies  will 
be  submitted. 


-  82  - 


-83- 


A  12730 


A  12732 


lico  li'odc  current  Collet.  tor 


4.0 


FACTUAL  DATA  (Cor.t'd) 


3.  Ion  Exchange  Membrane  Fuel  Cell  ~  "GPS"  Type. 

A  typical  "GPS"  ion  exchange  membrane  fuel  cell  is  shown  in  the 
center  of  Photograph  IV-4.  Two  catalyst  screen  electrodes  one  on 
each  side  are  in  intimate  contact  with  the  membrane. 

4.  Perforated  Metal  Screen  Hydrogen  Electrode.  Current  Collector 

A  most  recent  version  of  the  hydrogen  electrode  current  collector 
(at  the  top  of  Photograph  IV-4)  incorporates  the  perforated  screen 
air  electrode  currert  collector  configuration  with  manifolds  to  distri¬ 
bute  the  hydrogen  more  uniform  during  operation  of  the  cell  and 
hydrogen  side  inert  gas  purge  periods.  The  battery  separator 
material  was  not  added  to  the  particular  collector  photographed. 

5.  Assembled  Heat  and  MasB  Transfer  Air-Breathing  I.  E.  M. 

Fuel  Cell 

Photograph  IV-5  presents  a  "GPS"  ion  exchange  membrane  fuel  cell 
in  assembled  air  breathing  form,  i.  e.  attached  to  the  hydrogen 
electrode  current  collector.  Thermocouples  for  measuring  the 
hydrogen  side  collector  and  catalyst  temperatures  at  points  opposite 
those  on  the  air  side  shown  face  up  are  Bealed  between  the  frame 
and  membrane.  Previous  knowledge  cf  the  temperature  profile 
across  the  catalyst  surface  indicated  that  the  temperature  at  the 
extreme  sides  of  the  catalyst  are  the  same;  hence,  only  one  thermo¬ 
couple  on  one  side  is  needed.  Tne  reaction  zone  catalyst  area  has 
been  marked  off  with  tape  to  the  air  electrode  geometric  area  so  as 
to  give  well  defined  catalyst  active  areas  for  all  analyses  and  accurate 
current  density  calculations. 

6.  Assembled  Single  Cell  I.  E,  M,  ard  Fixture 

The  assembled  air-breathing  I.  E.  M.  fuel  cell  is  placed  on  top  of  one 
current  and  backup  plate  and  then  the  sir  electrode  current  collector 
is  placed  with  ribs  down  on  the  air  side  catalyst  and  finally  a  second 
current  and  backup  plate  is  placed  cn  top  of  it.  Sufficient  mechanical 
pressure  is  a.pplied  by  means  of  small  threaded  tie-rod  screws  to 
give  minimum  cell  electrical  resistance.  The  inlet  air  end  is  shown  in 
the  photograph.  (Ref.  Photograph  IV- 6).  This  unit  is  then  placed  in  the 
white  insulated  pressure,  tight  box  shown  in  Photograph  IV-7  and  external 
load  connections  are  made  by  means  of  terminals  in  tne  cover  which 
are  attached  to  the  cell  current  plates. 
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4.  4.2.2  B.  Heat  and  Mass  Transfer  Investigations  With  Electro* 
chemical  Reaction 

The  experimental  apparatus  used  for  heat  and  mass  transfer  studies 
with  electrochemical  reaction  is  shown  in  Photographs  IV-7  and  1V-8. 
Air  at  70°  -  80°F  and  0%  relative  humidity  is  fed  from  a  compressed 
air  bottle  at  a  few  inches  of  water  pressure  to  a  calibrated  flowmeter 
and  then  to  the  white  insulated  fuel  cell  fixture.  For  higher  inlet  air 
humidities,  the  air  is  passed  from  the  flowmeter  through  a  gas 
saturator  (as  shown  in  Photograph  IV-7)  and  then  to  the  fixture. 

Higher  air  temperatures  are  obtained  by  controlling  the  heat  input 
to  a  short  glass  section  wrapped  with  electrical  heating  tape  by 
means  of  a  rheostat  and  a  thermocouple.  Higher  cell  temperatures 
are  obtained  for  a  given  current  density  in  a  similar  fashion.  Reacted 
air  plus  product  water  vapor  leaving  the  fixture  passes  through  a  14 
long  x  2"  diameter  silica  gel  bed  to  a  glass  control  trap  immersed  in 
an  acetone-dry  ice  bath  at  -70°C  located  in  a  Dewar  bottle.  Boil-off 
from  the  bath  due  to  heat  conduction  is  minimized  by  insulating  the 
top  of  the  bottle  with  cork  wrapped  with  aluminum  foil. 

The  water  balances  were  achieved  by  weighing  the  cell  fixture,  the 
silica  gel  column,  the  cold  trap,  the  hydrogen  bubbler,  and  in  the 
greater  than  0%  relative  humidity  cases,  the  air  saturator  before 
and  after  operation  of  the  fuel  cell  at  stable  current  densities  of  about 
42,  62,  and  83  amps/ft2  for  periods  of  up  to  eight  hours  at  an  air 
flowrate  sufficient  to  remove  100%  of  the  product  water.  This  air 
flowrate  was  determined  by  monitoring  the  cell  resistance  by  means  ol 
a  Kordesch-Marko  bridge.  The  air  flowrate  at  which  the  cell  resis¬ 
tance  and  performance  do  not  vary  with  time  is  the  100%  removal 
point.  The  relative  humidity  of  the  inlet  air  was  quantitatively 
determined  by  passing  the  bottled  air  through  the  absorption  column, 
by  the  loss  in  weight  of  the  air  saturator,  the  air  flowrate,  and  the 
temperature.  Humidity  sensors  were  also  used.  The  relative 
humidity  of  the  reacted  air  stream  was  determined  from  the  temper¬ 
ature  of  the  gas,  the  change  in  weight  of  the  fixture,  absorption 
column,  and  the  cold  trap,  and  the  air  flowrate  by  using  Psychrometric 
charts.  The  amount  of  product  water  produced  was  calculated  from 
the  current  and  the  length  of  the  run  in  hours.  The  over-all  water 
balances  were  95  to  99%  complete. 
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The  heat  transfer  temperature  data  were  determined  using  a  Honeywell- 
Rubicon  precision  millivolt  potentiometer  (range  0-80  mv)  and  0.004  in. 
Pt  -  Pt  -  10%  Rh  thermocouples  insulated  with  a  coating  li  a  fluoro¬ 
carbon  resin.  These  data  consisted  of  measurement  of  the  tempera¬ 
tures  across  the  length  and  width  of  the  air  side  catalyst,  the  middle 
of  a  channel  to  the  edge,  the  collector,  and  the  corresponding  tempera¬ 
tures  on  the  hydrogen  side.  The  inlet,  ambient,  rib  outlet,  exit  bulk, 
and  inlet  abosrption  column  air  temperatures  were  also  measured  as 
were  the  face  temperatures  of  the  fixture.  Data  were  recorded  every 
10  to  15  minutes. 

Steady- state  attainment  was  determined  by  n  zero  rate  of  change  of 
these  temperatures.  Usually  about  1/2  to  1-1/2  hours  were  required 
for  conditions  to  steady-out. 

The  Kordesch-Marke  and  D.  C.  voltages  were  read  on  a  Honeywell- 
Rubicon  potentiometer  for  cell  resistance  data  while  a  continuous 
record  of  D.  C.  voltage  was  obtained  with  a  Varian  recorder. 

The  D,  C.  voltage  could  be  maintained  within  -  10  to  25  m.v.  by 
manual  control  of  the  flowrate  but  attempts  at  automatically  controlling 
the  cell  performance  were  relatively  unsuccessful  resulting  in  a  pro¬ 
gressively  more  rapid  decrease  in  cell  voltage  for  some  yet  un¬ 
explained  reason.  In  the  interest  of  obtaining  heat  and  mass  transfer 
data  under  steady  conditions  manual  control  was  employed. 

4.  4.  2.  3  Rates  of  Water  Transport  Through  and  Evaporation  from  the  "GPS" 

I.  E.  M,  Fuel  Cell,  Membrane,  Catalyst-Electrode,  and  Pure  Water, 

A  comparison  of  the  relative  rates  of  water  transport  through  and  from 
a  "GPS"  I.  E.  M.  fuel  cell  (membrane  +  2  catalyst  electrodes)  and  a 
"GPS"  membrane  alone  was  obtained  by  first  taking  the  air-breathing 
cell  assembly  with  the  hydrogen  side  filled  with  water  and  measuring  the 
loss  in  weight  of  the  cell  assembly  sitting  on  the  pan  of  a  precision 
Mettler  balance  (-O.OOOlg)  under  constant  ambient  conditions  of 
70-72°F  and  56%  R.  H.  every  ten  minutes  until  a  steady  rate  was 
attained.  Then  the  same  experiment  was  performed  using  a  piece  of 
membrane  from  the  same  batch  number  taped  to  a  hydrogen  side 
collector  and  frame.  The  relative  humidity  was  determined  from 
the  wet  and  dry  bulb  temperature  readings  and  the  total  atmospheric 
pressure  using  Psychrometric  Tables.  The  temperature  of  the 
ambient  air  and  the  catalyst  and  membrane  surfaces  were  determined 
by  thermocouples  both  at  the  start  and  end  of  a  run.  By  filling  the 
hydrogen  gas  compartment  with  water  to  a  level  just  equal  to  the 
height  of  the  frame,  the  rate  of  evaporation  of  water  from  a  pure  water 
(distilled)  surface  under  the  same  experimental  conditions  was 
determined  for  comparison. 
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4.  4,  2.  4  The  Effect  of  Air  Velocity  on  the  Water  Transport  and  Evaporation 
Rates  for  a  "GPS"  I.  E.  M.  Fuel  Cell 

The  rates  of  water  transport  through  and  evaporation  from  the 
surface  of  a  "GPS"  I.  E.  M.  fuel  cell  as  a  function  of  air  flowrate 
were  determined  by  filling  the  hydrogen  gas  compartment  of  the 
air-breathing  cell  assembly  completely  with  water  (no  air  bubbles) 
and  then  tying-off  the  inlet  and  outlet  tubes.  The  cell  assembly 
was  weighed  and  then  put  into  the  insulated  white  fixture  and  air 
at  0%  R.  H.  and  72-80°F  was  passed  into  the  fixture  at  flowrates 
which  spanned  the  range  of  the  experimental  fuel  cell  operating 
conditions  and  then  into  the  silica  gel  absorption  column  and 
cold  trap.  The  inlet  and  outlet  air,  the  catalyst  surface,  the 
collector,  and  the  fixture  face  temperatures  were  recorded  every 
ten  minutes.  The  highest  flowrate  run  was  terminated  when  the 
temperatures  did  not  vary  with  time  (3.  27  hours)  and  all  other 
runs  at  lower  air  flowrates  were  of  this  duration.  The  cell  assembly, 
the  fixture,  the  absorption  column,  and  cold  trap  were  weighed  on  a 
precision  Seko  balance  (5000g  capacity  -  0.  0005g)  after  each  run. 

Four  air  flowrates  corresponding  to  about  3  to  30  times  the  stoichio¬ 
metric  air  requirement  for  a  cell  running  at  a  current  density  of 
81  amps/ft2  were  used. 

4.  4.  3  Results: 

4.  4.  3.  1  A  Comparison  of  the  Hydrogen-Force  Convected  Air  Performance  of 
"GPS"  Ion  Exchange  Membrane  Fuel  Cells  with  the  New  and  Old 
Design  Perforated  Screen  Air  Electrode  Current  Collector. 

Increasing  the  percentage  rib-catalyst  contact  area  from  about  1%  to 
about  5%  and  changing  the  angle  of  the  rib  face  from  80  to  45  have 
resulted  in  a  fuel  cell  fixture  which  has  an  electric  resistance  of 
about  0.  00  1a  compared  to  about  0.  020a  for  previous  designs  and; 
therefore,  a  considerably  higher  performance  capability.  A  compari¬ 
son  between  the  previous  maximum  hydrogen-force  convected  air 
performance  possible  with  "GPS"  ion  exchange  membrane  fuel  cells 
and  the  present  shows  that  at  the  100%  product  water  removal  condition 
the  former  "hot"  (120°F)  performance  is  about  equal  to^he  latter  cold 
(85°F)  performance,  i.e.  0.69  D.  C,  volt  at  50  amps/ft  and  0.^59 
D.  C.  volt  at  100  amps/ft2.  However  the  new  design  "hot"  (120  F) 
performance  is  Q.  75  and  0.  68  volt  D,  C.  at  current  densities  of  50 
and  100  amps/ft  .  Current  densities  and  D.  C.  voltages  in  this  range 
have  been  maintained  for  periods  up  to  eight  hours  on  a  few  different 
cells.  The  portions  of  the  polarization  curves  above  110  amps /ft 
are  shown  in  broken  line  form  because  the  points  are  short  time  points. 
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i.e.  15  minutes  at  200  A/ft2  to  1  minute  at  290  amps/ft  .  The 
limiting  current  density  for  hydrogen  and  air  at  100%  removal  of  ^ 
product  water  appears  to  be  290  amps/ft  compared  to  155  amps/ft 
for  the  previous  design. 

The  maximum  dead-ended  hydrogen-oxygen  performance  in  a  standard 
laboratory  fixture  is  0.  82-0.  87  and  0.  70  to  0.  80  volt  D.  C.  at  current 
densities  of  50  and  100  amps/ft2,  respectively.  (Ref.  Figure  IV- 1). 


4.  4.  3.  2  Heat  and  Mass  Transfer  Analysis 

Although  the  heat  transfer  data  has  been  taken,  the  reduction  of  data 
and  analysis  was  incomplete  at  the  time  of  the  writing  of  this  report 
and  will,  therefore,  be  included  in  the  Final  Report. 

The  details  of  a  mass  transfer  analysis  of  the  possible  rate  limiting 
mechanisms  in  the  diffusion  of  oxygen  from  an  air  stream  to  the 
surface  of  the  oxygen  electrode  of  a  "CPS’'  ion  exchange  membrane 
fuel  cell  appear  in  Appendix  IV.  The  main  results  of  the  analysis 

are  as  follows: 

1.  CounterdiffuBion  of  Oxygen  and  Water  Vapor  in  the  Presence  of 
Inert  Nitrogen  Gas  -  Gas  Film  Controlling.  (45,  4fe,  47,  48,  56) 

The  rigorous  treatment  of  the  simultaneous  diffusion  of  two  gases 
in  the  presence  of  a  stagnant  film  of  inert  gas  leads  to  complex 
equations  even  for  this  relatively  simple  case.  However,  useful 
approximations  are  obtained  by  simple  relationships  which  are  based 
upon  the  assumption  that,  in  a  complex  system  of  diffusing  gases, 
the  diffusional  gradient  established  for  any  component  A  is  equal  to 
the  sum  of  the  gradients  which  would  result  from  the  separate 
diffusion  of  A  with  each  of  the  other  components  in  separate  binary 
systems  in  which  the  concentration  and  rate  are  the  same  as  in  the 
complex  system.  Thus  a  complicated  equation 
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for  four  diffusing  gases  and  an  inert  gas  I  reduces  to  the  form 
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which  can  be  applied  to  any  diffusing  component. 

a.  Based  on  Bulk  Gas  Diffusion  Coefficients  (45,  46,  47,  48,  56) 

On  the  basis  of  bulk  diffusion  coefficients  alone,  the  theoretical  limiting 
current  density  for  the  new  design  perforated  metal  air  electrode 
current  collector  operating  on  air  and  hydrogen  with  100%  product 
water  removal  is  13,  300  amps/ft2,  based  on  the  maximum  bulk  rate 
of  diffusion  of  oxygen.  The  current  density  at  which  the  maximum 
diffusion  path  for  water  vapor  becomes  equal  to  one  half  the  air  gap 
is  2,  430  amps/ft2. 

The  effect  of  the  counterdiffusion  of  water  vapor  on  the  rate  of  diffusion 
of  oxygen  is  negligible,  i.e.,  about  1.5%.  However,  the  rate  of  water 
vapor  diffusion  is  reduced  12.  1%  by  the  counterdiffusion  of  oxygen. 

b.  Based  on  Pore  Gas  Diffusion  Coefficients  (50,51,  55) 

The  problem  with  calculating  rates  of  gas  diffusion  in  pores  is 
determining  whether  the  flow  is  in  the  laminar,  intermediate,  or 
Knudsen  regimes  since  the  diffusion  coefficients  depend  on  the  mode 
of  flow.  Thus  for  laminar  flow  the  bulk  diffusion  coefficient  for  gases 

(  0.  707  \  l  sv~  \ 

Db  =  V  3  1L  7  CT  /  <IV-15> 

has  the  range  for  most  gas  mixture  at  atmospheric  pressure  and  20°C 
of  from  0.  1  to  1.0  cm2  .  It  is  dependent  on  the  total  pressure, 

the  total  concentration  of  molecules  and  the  molecular  size  of  the 
molecules  in  the  mixture.  For  the  Knudsen  flow  regime  the  gas  pore 
diffusion  coefficient 

2r  AT" 

Dk  =  - - -  (IV- 18) 

has  the  range  of  about  0.01  to  1.0  for  0.  00  \JL  and  lyU  pore 

radii,  respectively. 
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Therefore,  an  expression  which  shows  both  the  transition  between 
Knudsen  and  laminar  flow  and  reduces  to  these  two  special  cases 
when  r  X  and  r  ^  X  is  needed  to  calculate  the  over-all  pore 
diffusion  coefficient.  Equation  (IV- 19)  in  the  form 


D 


(0.  A. ) 


-2r  V/  3Dj}j 


(IV- 19) 


does  just  this.  The  effect  of  average  pore  radius  on  the  over-all 
pore  diffusion  coefficient  for  the  standard  G.  E.  ion  exchange 
membrane  fuel  cell  catalyst- screen  electrode  is  shown  in  Figure  IV-2. 
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for  the  catalyst  electrode. 

Based  on  the  average  oxygen  pore  diffusion  coefficient  alone,  the 
theoretical  limiting  current  density  is  9,  600  amps/ft  for  the  new 
design  perforated  metal  screen  air  electrode  at  100%  removal  of 
product  water.  Theoretically,  the  current  density  at  which  water 
vapor  diffusion  from  the  bottom  of  the  pores  to  the  air  stream  becomes 
the  limiting  factor  in  the  rate  of  removal  of  product  water  from  the 
pores  is  1,  450  amps/ft2.  (Ref.  Figure  IV-2). 

2.  Counterdiffusion  of  Oxygen  and  Water  Vapor  in  the  Presence  of_ 
Inert  Nitrogen  Gas  -  Liquid  Film  Controlling  (45,  46,  47). 


Using  Fick's  law  in  the  form 
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for  the  rate  of  diffusion  of  oxygen  through  water,  the  limiting  current 
density  for  a  "GPS"  I.  E.  M.  fuel  cell  operating  on  hydrogen  and  force 
convected  air  at  100%  product  water  removal  appears  to  be  from 
70  amps/ft2  for  a  non  wet-proofed  electrode  to  a  maximum  of 
500  amps/ft2  for  a  wet-proofed  electrode.  Poorly  or  non-uniformly 
wet-proofed  electrode  limiting  current  densities  should  fall  in 
between  these  two  limits. 

At  a  current  density  of  81  amps/ft2  the  reaction  zone  appears  to 
occupy  about  0.031%  of  the  total  pore  geometric  surface  area  and 
from  6.  4%  to  0.  006%  of  the  pore  length  for  the  lb/*-  and  0.  08ytc- 
average  pore  sizes,  respectively.  The  thickness  of  the  liquid  film 
is  about  8.  7  x  10'^cm. 

4.4.  3.  3  Rates  of  Water  Transport  Through  and  Evaporation  Rates  from  the 
"GPS"  I.  E.  M.  Fuel  Cell,  Membrane,  Catalyst-Electrode,  and 
Pure  Watfer.  (53) 

The  results  of  the  measurement  of  the  steady-state  rates  of  water 
evaporation  from  the  surfaces  of  a  "GPS"  I.  E.  M.  fuel  cell,  mem¬ 
brane,  and  pure  water  under  the  same  ambient  conditions  of 
70°-72°F  air  at  56%  relative  humidity  are  shown  in  Table  IV- 1. 

The  ambient  air  velocity  varied  randomly  from  0.25  to  1.0  ft/sec. 

TABLE  IV- 1 


Evaporation 

(gH20/Hr) 

Rate 

Geom.  Area 
(cm2) 

kH20  in  ,  1 

g-moles  H20/cm^-sec 

1 

[kH2C>)  w 

ph2o 

(kH2°)  i 

v.p.  HzO 

Water 

346 

104.6 

1. 98x10*7 

1 

1 

Membrane 

0.  785 

74.  2 

1. 63x10 ”7 

0.822 

0.86 

Fuel  Cell 

0.  808 

74.  6 

1. 672xl0'7 

0.845 

0.  962 

From  the  data  the  evaporation  rate  of  water  from  a  fuel  cell  (CRM  7) 
is  84.  5%  of  that  observed  with  pure  water  and  the  effective  water 
vapor  partial  pressure  at  the  surface  is  96.  2%  of  the  vapor  pressure 
of  water  at  the  surface  temperature.  For  the  membrane  (CRM  7;  I.  E.  C. 
=  0.  64;  %  water  dry  basis  =  18)  alone  the  water  evaporation  rate  is  only 
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82.  2%  of  the  pure  water  rate  and  the  partial  pressure  of  water 
vapor  at  the  surface  is  86%  of  the  vapor  pressure  of  water  at 
the  surface  temperature.  For  the  membrane  and  fuel  cell  the 
evaporation  rates  are  equal  to  the  water  transport  rates  through 
the  particular  medium  also. 

Since  the  rate  of  water  transport  and  evaporation  rate  changes  only 
slightly  from  the  fuel  cell  to  the  membrane  case,  the  transport  and 
evaporation  properties  of  the  catalyst-electrode  are  assumed  to  be 
the  same  as  for  pure  water. 

4  4  3  4  The  Effect  of  Air  Velocity  on  the  Water  Transport  and  Evaporation 
Rates  for  a  "GPS"  I.  E.  M.  Fuel  Cell. (49,  53) 

The  effect  of  air  velocity  on  the  steady- state  water  transport  and 
evaporation  rates  for  a  "GPS"  L  E.  M.  fuel  cell  in  the  absence  of 
electrochemical  reaction  (No  E.  C.  R.  )  is  shown  in  Figure  IV- 3. 

The  range  of  inlet  air,  outlet  air,  and  catalyst  temperatures  were 
70°  to  80°F,  68°  -  73°F,  and  66°-  73°F  respectively.  The  inlet 
air  relative  humidity  was  0%  and  the  average  relative  humidity  in 
the  cells  at  steady- state  conditions  was  from  28%  to  44%for  the 
highest  to  lowest  flowrates  which  represented  from  3  to  30  times 
the  stoichiometric  air  requirement  for  a  cell  operating  at  a  current 
density  of  81  amps/ft2.  The  data  were  obtained  with  the  fuel  cell 
in  the  insulated  white  fixture  using  the  new  design  perforated  metal 
screen  air  electrode  current  cpllector.  The  data  were  corrected 
for  the  variable  relative  humidity  conditions  by  using  the  data 
from  section  4.  4.  3.  3  and  dividing  the  steady- state  rates  by  the 
water  vapor  partial  pressure  gradient  from  the  catalyst  to  the  air 
stream  to  obtain  the  water  transport  and  evaporation  coefficients 
kH  o-  The  data  are  rpPresented  by  the  empirical  equation 

-6  0.76 

kH2o  =  2-  x  10  V 

g-moles 

where  ku  q  =  - 7 - 

n2^  cm  -sec-atm 

(cm  \ 

-  J 

S  6  C  J 

The  data  point  used  in  the  mass  transfer  analysis  (Ref.  Appendix  IV) 
using  a  different  "GPS"  I.  E.  M.  fuel  oell  in  the  same  hardware  but  at 
a  current  density  of  81  amps/jt2  (E.  C.  R. )  yields  a  steady-state 
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evaporation  rate  at  the  100%  product  water  removal  point  which  is 
very  close  to  the  data  taken  with  no  E.  C.  R.  Correction  for  the  effect 
of  the  counterdiffusion  of  oxygen  brings  the  data  more  closely  to  the 
no  E.  C.  R.  data.  The  average  air  temperature  was  90°F,  the  average 
percent  relative  humidity  was  35%,  and  the  catalyst  temperature  was 
115°F  for  the  E.  C.  R.  point. 

4.  4.  4  Discussion  of  Results 

4.  4.  4.  1  Performance 

The  main  factors  contributing  to  the  markedly  improved  performance 
observed  with  the  new  design  perforated  metal  screen  air  electrode 
current  collector  operating  on  hydrogen  and  force-convected  air  at 
70°-80°F  and  0%  R.  H.  are  believed  to  be  the  following:  much  more 
uniform  removal  of  product  water  with  the  45°  rib  face  angle  as 
».  indicated  by  a  maximum  variation  of  catalyst  temperature  of  1  to  2°F 

at  a  current  density  of  62  amps/ft^  and  5°F  at  8  1  amps/ft2  compared 
to  10°  to  20°F  at  50  amps/ft2  for  the  old  design;  lower  cell  resistance; 
higher  air  side  heat  transfer  coefficients;  and  higher  current  carrying 
capability.  With  the  present  low  rib-catalyst  contact  area,  the  mem¬ 
brane  resistance  is  the  total  resistance.  Further  study  above  110  amps/ 
ft*"  at  a  constant  temperature  of  about  120°F  and  100%  water  removal 
is  warranted  to  prove  out  high  current  density  capability. 

4.  4.  4.  2  Mass  Transfer  Analysis: 

The  results  of  the  mass  transfer  analysis  which  covered  both  bulk  and 
pore  gas  diffusion  accounting  for  the  counterdiffusion  of  oxygen  and 
water  vapor  and  gas-liquid  diffusion  indicate  that  the  maximum  current 
density  possible  with  the  present  perforated  air  electrode  current 
collector  operating  with  100%  product  water  removal  is  limited  by  the 
rate  of  diffusion  of  oxygen  through  a  thin  lilm  of  water  at  the  catalyst- 
electrolyte  interface  or  reaction  zone  which  covers  about  0.03%  of  the 
total  pore  geometric  surface  area.  The  estimated  range  of  the  limiting 
current  density  is  70  to  500  amps/ft2.  The  highest  experimentally 
observed  limiting  current  density  is  290  amps/ft  2  with  a  catalyst 
containing  about  50%  by  volume  wet-proofing  agent.  Air  flowrate 
must  be  of  sufficient  magnitude  to  remove  100%  of  the  product  water 
from  the  oxygen  electrode  and  the  air  reaction  passages.  Water  vapor 
diffusion  from  the  pores  does  not  become  limiting  until  a  current  density 
of  1,  450  amps/ft2  is  reached,  a  value  which  is  impossible  to  reach  with 
air.  However,  the  maximum  pure  oxygen  performance  at  100%  removal 
of  product  water  by  recirculation  of  the  stream  and  absorption,  conden- 
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sation,  or  mechanical  separation  of  product  water  is  estimated  to  be 
from  330  to  2,  380  amps /ft2  based  on  liquid  film  control.  Since 
water  vapor  diffusion  from  the  pores  is  limiting  above  1,450  amp/ft  , 
this  current  density  would  appear  to  be  the  maximum  value  for  pure 
oxygen  under  ideal  conditions  at  temperatures  up  to  130°F. 

On  this  basis,  the  present  perforated  metal  screen  air  electrode 
current  collector  has  high  current  density  capability  and  is  suitable 
for  large  size  cells.  The  mode  of  removal  of  product  water  is 
adaptable  for  pure  oxygen  systems  also. 

The  self -diffusion  type  air  electrode  is  limited  by  bulk  gas  diffusion 
coefficients  to  a  maximum  round  shape  of  7.0  and  3.5  in.  in  diameter 
at  current  densities  of  81  and  162  amps/ft  ,  respectively.  The  corres¬ 
ponding  air  gap  widths  are  ^  0.007  in.  ^  0.632  in.  and  ^  0.024  in.  ^ 

0.  316  in.  Therefore,  this  type  device  is  limiting  to  small  sizes  and 
low  current  densities. 

4.  4.  4.  3  Comparison  of  Water  Transport  and  Evaporation  Coefficients 

This  data  indicates  that  the  "GPS"  type  membrane  fuel  cell  transports 
and  evaporates  water  into  an  air  stream  at  about  85%  of  the  pure  water 
rate.  The  suppression  of  water  vapor  pressure  at  the  surface  is  about 
0.96  times  the  pure  water  value.  Data  in  the  literature  (52)  indicate 
that  as  little  as  1%  porosity  in  a  plastic  film  which  does  not  contain 
water  will  cause  an  evaporation  rate  equal  to  90%  of  the  water  rate. 

Since  the  I.  E.  M.  contain  up  to  20%  water  the  experimental  data  are 
entirely  reasonable.  The  catalyst  layer  offers  essentially  no  resistance 
to  the  transport  of  water  and  acts  like  a  pure  water  surface. 

These  results  are  in  agreement  with  the  fact  that  dry  I.  E.  M.  cells 
(0.  250-*)  have  been  wet  with  product  water  by  operating  at  a  low  air 
flowrate  (wetting  condition)  until  the  cell  resistance  reduces  to 
normal  values  (0.015*).  By  flowing  air  at  a  high  flowrate,  wet  cells 
can  be  dried  very  quickly  as  well  known. 

4,  4.  4.  4  Air  Velocity  Effects  on  Water  Transport  and  Evaporation  Rate  (46) 

The  empirical  equation  ^  o  76  represents 

kn20  =  2.95  x  10  V  ‘ 

the  water  transport  rate  through  and  evaporation  from  a  ''GPS1'  1.  E.  M. 
fuel  cell  in  the  absence  of  and  very  likely  with  electrochemical  reaction 
is  in  agreement  with  kH_0  =  ,c  V0,  5  for  laminar  air  flow  parallel  to 
‘many  surfaces,  kj^O  =  C  V°’  65  for  laminar  flow  perpendicular  to 


-  104  - 


4.0 


FACTUAL  DATA  (Cont'd) 


various  types  of  surfaces,  and  for  laminar  air  flow  past  a  surface 
facing  downstream  at  an  angle  of  40°  kH20  =  C  V  *  which  is 
the  maximum  value.  With  the  45°  air  impingement^i^le  a  value  of 
V°*  might  be  expected  on  the  inlet  air  side  and  V  '  on  the  exit 
side  giving  an  average  value  of  V0-70.  The  slightly  higher  value  is 
probably  due  to  the  effect  of  the  perforation  close  to  catalyst  surface. 
The  close  agreement  between  the  no  E.  C,  R,  and  E.  C.  R,  data  further 
substantiate  the  observation  that  the  catalyst  layer  offers  no  addi¬ 
tional  resistance  to  water  transporter  evaporation  as  was  first  shown 
under  C. 


In  view  of  the  nature  of  the  empirical  equation,  the  conclusions 
presented  under  section  4.  4.  4.  2,  and  previous  mass  transfer 
analysis  (54)  of  binary  gas  systems  on  phenolic  type  I.  E.  M.  fuel 
cells  at  lower  current  densities,  there  exists  a  distinct  possibility 
that  the  over-all  mass  transfer  mechanism  for  the  oxygen  electrode 
may  be  represented  by  a  modified  mass  transfer  factor  of  the  iorm. 

^4=  C  Re'0,8JlO2_1'0 


where  C  is  a  const.,  Re  is  the  Reynold's  number,  and  q.,  is  the  mole 
fraction  of  oxygen  gas  in  the  liquid  film  at  the  catalyst-electrolyte 
interface.  Using  the  empirical  equation,  an  expression  of  the  form 


C  k 


-1.  315 

h2o 


-b 

2 


may  apply  for  force  convected  air  systems  or  any  recirculating  gas 
system  which  removes  100%  of  the  product  water  from  the  oxygen 
electrode  for  all  practical  current  densities.  If_  a  pure  oxygen 
performance  of  1,  450  amps/ft2  could  be  obtained,  the  rate  of  water 
evaporation  would  no  longer  depend  on  the  gas  recirculation  rate. 

The  possibility  of  the  existence  of  a  mass  transfer  correlation  of 
this  form  will  be  explored  and  reported  in  the  Final  Report  as  will 
a  similar  heat  transfer  correlation  and  an  over-all  correlation  of 
heat  and  mass  transfer  mechanisms,  if  possible  in  the  time  remaining 
under  the  contract. 

On  the  basis  of  the  water  transport  and  evaporation  coefficients 
obtained  for  "GPS"  membranes,  they  may  be  used  in  a  system  to 
conserve  product  water  and  humidify  reaction  feed  air.  An  invention  disc 
losure  will  be  written  along  these  lines  in  the  immediate  future. 
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4.  4.  5  Specific  Conclusions: 

1.  The  new  design  perforated  metal  screen  air  electrode  current 
collector  has  permitted  substantial  increases  in  "GPS"-type  ion 
exchange  membrane  fuel  cell  performance  on  force-convected  air 
at  70* -100* F  and  hydrogen,  i.e.  0.  63-0.  68  volt  D.  C.  at  a  current 
density  of  81  amps/ft  and  a  limiting  current  density  which  appears 
to  be  290  amp s /ft1  mainly  due  to  more  uniform  evaporation  of  pro¬ 
duct  water  and  decreased  fixture  resistance. 


2.  The  results  of  a  mass  transfer  analysis  which  included  bulk  gas, 
pore,  and  gas-liquid  diffusion  mechanisms  indicate  that  the  limiting 
current  density  range  for  force-convected  air  at  100%  removal  of 
product  water  is  from  70  to  500  amps/ft2  depending  on  whether  the 
catalyst  is  non,  poorly,  or  well  wet-proofed.  The  maximum  current 
densities  appear  to  be  limited  by  the  rate  of  diffusion  of  oxygen  through 
a  thin  aqueous  film  (which  covers  about  0.03%  of  the  total  geometric 
pore  surface  area)  at  the  catalyst-electrolyte  interface.  The  corres¬ 
ponding  limiting  current  density  range  for  pure  oxygen  appears  to  be 
330  to  2400  amps/ft2  at  100%  removal  of  product  water  but  may  be 
limited  further  to  1,450  amps/ft2  by  water  vapor  diffusion  from 

the  pores.  On  this  basis,  the  new  design  perforated  screen  air 
electrode  current  collector  is  suitable  for  large  size  cells  and  oper¬ 
ation  at  high  current  densities. 

3.  A  comparison  of  the  rates  of  water  transport  and  evaporation  for 
"GPS"  I.  E.  M.  fuel  cell,  membrane,  and  pure  water  indicate  that  the 
fuel  cell  and  membrane  rates  are  about  85%  of  the  pure  water  rate 
and  that  the  catalyst  layer  rate  is  substantially  100%  of  pure  water. 

4.  The  effect  of  air  velocity  on  the  rate  of  evaporation  of  water  from 
a  "GPS"  I.  E.  M.  fuel  cell  with  or  without  electrochemical  reaction 
appears  to  be  represented  by  an  empirical  equation  of  the  form 

*H20  =  2-95x  lO'6  V0*76 


Based  on  past  and  present  work,  there  exists  the  possibility  that  the 
over-all  mass  transfer  mechanism  at  the  oxygen  electrode  may  be 
represented  by  a  modified  mass  transfer  factor  of  the  form. 


or 
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4.  4.  A  Addition  to  TASK  IV 

AN  OPTIMIZING  CONTROL  DEVICE  FOR  USE  IN  LABORATORY 
EVALUATION  OF  HYDROGEN -AIR  FUEL  CELL  OPERATION. 

T.  Hovious 

This  Section  represents  an  addition  to  contract  requirements  as  an  off¬ 
spring  of  investigations  on  water  transport  in  membranes,  internal  cell 
resistance  and  the  response  to  temperature  and  air  flow  rate. 

4.  4.  A.  1  General 

This  control  device  is  designed  to  adjust  air  flow  rate  to  an  operating 
hydrogen-air  fuel  cell  in  such  a  manner  as  to  maintain  near  optimum 
DC  output  voltage  at  any  given  current  level. 

4.  4.  A.  2  Control  Mode 

A  typical  fuel  cell  DC  voltage  versus  total  cell  water  content  is 
illustrated  by  Figure  4A-1.  Cell  voltage  for  a  given  current  output 
will  be  less  than  the  maximum  attainable  if  the  fuel  cell  is  either  too 
wet  or  too  dry.  The  solid  line  curve  of  Figure  4A-2  illustrates 
qualitatively  the  air  flow  rate  required  at  any  given  output  current  so 
that  the  water  removal  rate  just  equals  the  water  production  rate. 

An  air  flow  rate  greater  than  that  given  by  the  solid  line  curve  will 
decrease  the  cell  moisture  content,  whereas  a  smaller  air  flow 
rate  will  allow  the  cell  moisture  content  to  increase.  The  effects 
of  variation  in  quantities  such  as  cell  vapor  pressure  suppression, 
cell  operating  temperature,  and  inlet  air  relative  humidity  are 
neglected  in  the  foregoing  description. 

To  illustrate  the  mode  of  control  to  be  used,  consider  the  following 
assumption.  A  fixed  output  current  of  value  I  (Figure  4A-2)  is 
established.  The  control  would  be  adjusted  so  that  either  air  flow 
rate  A  or  air  flow  rate  B  (Figure  4A-2)  could  exist  for  the  current 
of  value  I.  Let  it  also  be  assumed  that  at  start-up  the  cell  moisture 
conditions  correspond  to  that  of  point  Wi  (Figure  4A-1)  and  that 
the  control  supplies  an  air  flow  rate  equal  to  point  A  (Figure  4A-2). 
Thus,  the  cell  water  content  is  too  large,  but  a  drying  condition  has 
been  applied  which  causes  the  total  cell  water  content  to  decrease 
with  time.  The  control  senses  this  positive  rate  of  change  of  output 
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voltage  and  maintains  the  drying  air  flow  rate  A  (Figure  4A-2). 

Eventually  the  cell  moisture  content  decreases  to  a  value  Wc 
(Figure  1}  and  the  output  voltage  reaches  its  maximum.  As  the 
cell  further  dries,  the  output  voltage  decreases  until  it  has 
reached  an  amount  equal  to  AV  (Figure  4A-1)  less  than  the  peak 
voltage  value.  This  corresponds  to  a  water  content  W2  (Figure  4A-1). 

At  this  point,  the  control  changes  the  air  flow  rate  to  the  wetting 
condition  (point  B,  Figure  4A-2).  The  cell  water  content  increases 
until  a  water  content  (Figure  4A-1)  is  reached. 

At  this  point,  the  control  flips  back  to  an  air  flow  rate  which  allows 
cell  drying.  Thus,  a  continuous  cycling  between  cell  water  content 
conditions  W2  and  W3  (Figure  4A-1)  is  maintained.  The  control 
will  be  able  to  maintain  AV  as  small  as  0.01  volts  if  desired. 

The  control  also  incorporates  a  timer  which  operates  whenever  the 
rate  of  change  of  voltage  is  negative.  If  the  rate  of  change  of  voltage 
does  not  become  positive  within  a  pre-selected  period  of  time,  the 
control  will  automatically  select  the  opposite  air  flow  rate  condition. 

This  ensures  attainment  of  the  proper  operating  mode  when  start-up 
or  transient  loading  upsets  the  normal  control  mode. 

4.  4.  A.  3  System  Description 

The  system  to  be  used  is  shown  in  the  block  diagram  of  Figure  4A-3. 

In  using  this  system,  the  desired  fuel  cell  output  current  level  is  first 
determined.  The  drying  and  wetting  air  flow  rates  (points  A  and  B  of 
Figure  4A-2)  corresponding  to  the  selected  current,  I,  are  then 
established.  The  setting  air  flow  rate  (point  B  of  Figure  4A-2)  is 
then  set  by  adjustment  of  valve  1  (Figure  4A-3)  with  the  solenoid 
operated  valve  closed.  The  solenoid  operated  valve  is  then  opened 
and  valve  2  (Figure  4A-3)  is  set  so  that  the  drying  air  flow  rate 
(point  A  of  Figure  4A-2)  is  established.  The  electronic  control  will 
then  operate  the  solenoid  valve  so  as  to  alternate  between  drying 
and  wetting  air  flow  rates  and  maintain  the  cell  water  content  such  as 
to  yield  maximum  output  DC  voltage. 

This  system  is  simple,  and  provides  for  optimizing  only  at  each  selected 
fuel  cell  output  current  condition.  It  does  not  incorporate  automatic 
adjustment  of  the  air  flow  rates  as  the  current  is  varied  over  a  wide 
range.  This  air  flow  rate  adjustment  must  be  manually  done  whenever 
a  major  change  in  output  current  is  made.  Automatic  compensation 
for  varying  electrical  loads  can  be  incorporated  in  the  control.  The 
major  complication  would  be  that  of  providing  a  means  of  modulating 
the  air  flow  proportional  to  an  electrical  signal. 
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Operation  and  Circuitry  of  the  Electronic  Control 
Operational  Control  Equations: 

The  operational  equations  for  amplified  A  are: 
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where  S  is  the  Laplace  one-sided  operator 


Similarly,  where  the  region  of  operation  is  outside  of  breakdown 
range  of  D  j  and  D2,  we  conceive  the  transfer  Vg/Vy^  .  Diode 
Dj  and  D2  are  shown  in  Figure  4A-4  which  illustrates  the  essentials 
of  the  control  unit. 
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Consequently,  the  overall  control  equation  from  Vg  to  Vj  is: 
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The  transfer  function  Vg/Vj  determines  the  signal  at  the  logic  selection 
diodes  D  3  and  D4,  The  gain  to  DC  voltage  change  is  5000  and  this 
determines  the  control  band  sensitivity.  A  change  of  10  millivolts 
causes  Vg  to  change  by  50  volts  which  is  more  than  sufficient  to  pull 
in  relay  CRj  or  CR  2  depending  on  the  sign  of  the  change  in  the  fuel  cell 
voltage.  If  the  fuel  cell  voltage  increases  by  10  millivolts,  the  relay 
CR2  is  selected  by  diodes  D4  and  pulls  in.  Pulling  in  CR2  causes  no 
change  in  the  control  sequence.  When  the  fuel  cell  voltage  drops  by 
10  millivolts,  diode  D3selects  CR  1  relay  and  the  control  sequence  begins. 
Pulling  in  CRj  causes  the  bistable  relay  LR  (Figure  4A-5)  to  change 
in  state  and  begins  to  supply  signal  power  to  the  time  delay  circuit 
composed  of  the  time  delay  relay  Tq,  the  logic  diode  D9  and  storage 
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element  C If  the  change  in  fuel  cell  voltage  reverses  sufficiently 
(i.e.  10  millivolts  or  greater)  to  drop  out  CR^,  then  the  control 
returns  to  its  primary  state.  If  the  fuel  cell  voltage  continues  to 
drop  off  for  a  period  greater  than  the  time  out  period  of  the  time 
delay,  then  the  time  delay  relay  pulls  in,  pulsing  the  bistable  relay 
LR  which  changes  the  state  of  its  contacts  which  changes  the  air 
flow  rate  through  activation  of  the  previously  mentioned  solenoid 
control  valve.  This  time  delayed  check  is  a  safety  device  to  insure 
proper  mode  of  operation  under  a  non  steady-state  condition.  This 
control  is  essentially  the  same  as  a  unit  designed  and  tested  to 
control  a  hydrogen-oxygen  fuel  cell  by  varying  mixed  water  coolant 
to  the  cell,  thus  controlling  evaporation  of  water  from  the  cell. 

However,  the  air-breathing  cell  in  the  operation  described  tends  to 
have  lower  rates  of  change  of  voltage  with  change  in  control  mode 
than  does  the  hydrogen-oxygen  cell  described;  thus  more  care  was 
necessary  in  evaluation  of  integrater  drift  problems  in  the  electronic 
circuitry. 

4.  5  TASK  V: 

DESIGN  PARAMETERS  AFFECTING  AIR-OPERATED  IEM  FUEL  CELLS 

.  D.  W.  Craft 

4.  5.  1  Introduction: 

The  optimum  operation  of  an  IEM  fuel  cell  system  on  air  requires  that 
the  following  conditions  be  approached: 

1.  Removal  of  the.  product  water  from  the  IEM  at  the  rate  of  generation. 

2.  Retention  of  an  optimum  water  content  within  the  IEM  to  attain 
maximum  voltage  output. 

3.  Minimization  of  temperature,  partial  pressure,  water  content  and 
voltage  gradients  within  the  IEM. 

An  experimental  investigation  has  been  initiated  to  determine  how 
successful  the  approach  of  using  a  recirculating  air  flow  within  a  fuel 
'cell  system  will  be  in  attaining  the  above  requirements. 

4.5.2  The  Recirculating  Air  Flow  Fuel  Cell  System 

The  fuel  cell  system  under  consideration  is  shown  in  Figure  V-l.  The 
gases  within  the  system  which  are  on  the  air  side  of  the  ion  exchange 
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membrane  (IEM)  are  continuously  recirculated  by  Fan  R  in  order 
to  minimize  temperature  and  partial  pressure  gradients  and  to 
eliminate  the  possibility  of  stagnating  the  flow  in  some  region 
of  the  fuel  cell  stack. 

The  gases  within  the  recirculation  system  will  be  referred  to  as 
air  in  this  discussion  since  the  constituents  are  the  same  as  air  but 
with  a  different  fractional  composition  due  to  the  fuel  cell  reactions. 

Fan  R  established  static  pressures  within  the  recirculating  system 
which  are  above  and  below  the  communicating  atmospheric  pressure, 
thus  allowing  the  intake  of  air  at  station  6  and  the  discharge  of  air  at 
stations  7  or  8.  The  discharge  of  air  at  station  7  is  preferable  from 
the  standpoint  that  the  oxygen  partial  pressure  is  a  minimum  at  this 
station  and  therefore  represents  a  minimum  flow  of  oxygen  from  the 
system  for  a  given  concurrent  flow  rate  of  water  vapor.  Discharge 
from  station  7  allows  the  establishment  of  a  higher  average  oxygen 
partial  pressure  within  the  system.  The  use  of  station  7  as  a  discharge 
location  has  additional  advantages  from  the  standpoint  of  analysis  in¬ 
asmuch  as  (1)  the  recirculation  system  can  be  directly  compared  with 
a  through-flow  system  (no  recirculation  or  as  though  there  was  a 
blockage  at  station  3),  and  (2)  the  bleed  conditions  are  unaffected  by 
the  rate  of  recirculation.  The  main  reason  the  recirculation  system 
can  minimize  tne  temperature  gradients  within  the  fuel  cell  stack  is 
that  the  mode  of  heat  rejection  is  by  the  evaporation  and  condensation 
of  water  that  is  retained  within  the  system.  The  relatively  high  rate 
of  mass  flow  in  the  recirculation  circuit  (i.  e.  100  times  the  stoichio¬ 
metric  requirement)  allows  the  necessary  water  for  heat  rejection  to 
be  conveyed  in  the  recirculation  loop  as  water  vapor.  The  air  leaving 
the  fuel  cell  stack  is  essentially  saturated  and  after  additional  air  is 
admitted,  mixed,  and  heated  by  Fan  R,  the  air  is  still  essentially 
saturated  as  it  enters  the  condenser.  As  the  recirculation  air  passes 
through  the  condenser,  the  external  cooling  air  reduces  the  temper¬ 
ature  of  the  gases  and  condensation  takes  place  accompanied  by  a 
small  temperature  drop  in  the  recirculating  air.  The  temperature 
change  is  small  due  to  the  latent  heat  of  vaporization  of  water. 

The  water  condensed  on  the  condenser  surfaces  is  absorbed  by  a 
porous  wicking  material  and  the  water  is  transported  to  the  evapora¬ 
tive  surfaces  of  the  fuel  cell  and  by  capillary  conduction  through  the 
same  wicking  material.  The  evaporative  surfaces  of  the  fuel  cell 
assembly,  as  referred  to  herein,  are  the  metal  current  collector 
plates  which  form  the  hydrogen  chamber  in  each  fuel  cell  assembly 
as  shown  in  Figure  V-2.  The  air  side  of  the  IEM  is  also  an  evapora¬ 
tive  surface,  but  this  surface  will  not  be  referred  to  as  such  herein. 
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inasmuch  as  it  is  planned  that  in  the  system  under  consideration,  the 
rate  of  evaporation  from  the  IEM  will  be  equal  to  the  rate  of  product 
water  generation  within  the  air  side  electrode. 

4.  5.  3  Product  Water  Management 


Early  air-operated  IEM  fuel  cells  required  that  water  be  periodi¬ 
cally  added  to  the  system  to  offset  the  excessive  drying  of  the 
membranes.  The  excessive  drying  was  caused  by  a  high  through- 
flow  rate  of  air  in  and  out  of  the  system.  The  high  through-flow 
rate  was  required  for  a  maximum  oxygen  partial  pressure  within 
the  system. 

The  recirculation  system  shown  schematically  in  Figure  V- 1  avoids 
the  necessity  for  periodic  water  addition  by  reducing  the  through- 
flow  rate  of  air  until  the  amount  of  water  leaving  the  system  is 
equal  to  the  rate  of  product  water  generation.  Appendix  V,  Sections 
1-6,  gives  an  analysis  whereby  it  is  shown  that  the  flow  rate  in  and 
out  of  the  system  (through-flow  rate)  is  essentially  dependent  upon 
the  fuel  cell  stack  operating  temperature  and  the  conditions  of  the 
entering  air. 

Figure  V-3  shows  the  variation  of  the  ratio  of  the  through-flow  to  the 
stoichiometric  air  flow  with  various  operating  conditions  when  the 
product  water  is  removed  from  the  system  at  the  same  rate  as  it  is 
generated.  The  data  in  Figure  V-3  was  calculated  from  Equation  V-22. 
Figure  V-3  shows  that  if  the  dry  air  entered  the  system,  the  air  flow 
requirement  would  only  be  dependent  on  altitude  (operating  pressure) 
and  the  stack  operating  temperature. 

The  relative  humidity  HR^of  the  entering  air  also  has  considerable 
influence  on  the  air  flow  requirement.  Figure  V-3  is,  however, 
plotted  on  reciprocal  coordinates  showing  an  infinite  air  flow  rate 
when  the  stack  temperature  is  equal  to  the  inlet  temperature  and 
HRl  =  100  percent.  The.  aforementioned  condition  is  unrealistic 

since  the  stack  temperature  would  be  above  the  inlet  temperature. 

Figure  V-4  shows  the  variation  of  the  oxygen  partial  pressure  at  the 
exit  of  the  fuel  cell  stack  for  various  air  inlet  and  stack  temperature 
conditions.  The  data  in  Figure  V-4  was  calculated  from  Equation  V-9. 
Although  Equation  V-9  is  relatively  complex,  it  is  shown  in  Section  5 
of  Appendix  5  that  a  linearized  version  of  Equation  V-9  (cf  Equation  V-12 
and  V-13)  is  quite  accurate. 
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The  linearized  versions  of  Equation  V-9  suggest  the  possibility  of 
providing  a  control  which  would  (1)  compute  the  oxygen  partial 
pressure  required  for  water  balance,  (2)  sense  the  true  oxygen 
partial  pressure,  and  (3)  adjust  the  through-flow  rate  to  maintain 
the  proper  water  balance.  The  major  difficulty  of  applying 
Equation  V-13,  for  example,  to  a  control  scheme  is  the  complexity 
of  accurately  measuring  the  vapor  pressure  of  the  water  in  the 
entering  air. 

A  more  practical  approach  which  would  avoid  the  necessity  of 
measuring  the  vapor  pressure  of  water  of  the  entering  air  would 
be  to  alternately  produce  a  wetting  and  drying  condition  within 
the  system  by  controlling  the  air  flow  rate.  Whether  the  air  flow 
rate  should  be  increased  or  decreased  could  be  determined  by  a 
measurement  of  the  fuel  cell  Btack  internal  resistance  as  a  criterion 
for  optimum  IEM  water  content. 

4.  5.  4  Operation  oi  the  IEM  Fuel  Cell  in  a  Uniform  Environment 

The  main  benefit  of  the  recirculation  of  air  within  the  fuel  cell 
system  under  consideration  is  to  minimize  the  temperature, 
partial  pressure,  water  content,  and  voltage  gradients  within  the 
IEM.  The  example  given  for  the  heat  balance  across  the  fuel  cell 
stack  in  Section  10  of  the  Appendix  A  shows  the  following  conditions 
entering  and  leaving  the  stack  for  a  recirculation  rate  f  c  =  1°°  tlme 
the  stoichiometric  requirement: 

(P02)  1  = 

(f>h2o)  {  =  1.  6023  psia 
tl  =  118. 05F 

In  a  fuel  cell  system  with  the  same  through-flow  rate  <rc  =  3.42  as 
shown  by  Figure  V-3,  the  oxygen  partial  pressure  at  the  stack 

entrance  would  be  (p0z)  i  =  3-016  Psia  as  shown  in  FiSure  v’4- 
Also,  if  all  the  cooling  were  done  in  the  system  without  recircula¬ 
tion  with  a  through-flow  rate  <rc  =  100  and  no  evaporation,  the 
temperature  rise  would  be  approximately  13  deg.  F  rather  than 
2  deg.  F  as  shown  in  the  above  example. 


M 

(P»2^ 


2  =  2.  058  psia 


=  1.  6924  psia 


t2  =  120  F 
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The  greatest  difficulty  in  any  air  operated  fuel  cell  system  is  to 
avoid  a  non-uniform  evaporation  rate  from  the  IEM  since  the 
rate  of  water  redistribution  within  the  IEM  is  slow  and  local 
drying  or  wetting  is  to  be  avoided.  Even  though  the  above  values 
of  the  partial  pressure  of  water  vapor  seem  to  differ  slightly,  the 
driving  force  to  evaporate  water  from  the  IEM  is  higher  at  the 
entrance  of  the  fuel  cell  stack  than  it  is  at  the  exit.  The  differ¬ 
ence  is  to  be  offset  by  using  a  mask  which  has  graduated  perfor¬ 
ations  spaced  a  distance  away  from  the  IEM  and  is  to  be  used  for 
controlling  the  evaporation  rate,  thus  attaining  uniformity.  The 
use  of  a  fixed  mask  only  permits  uniformity  of  a  ,-apn  ration  from 
the  IEM  at  one  current  density  in  a  given  system.  Figure  V-5 
shows  a  hypothetical  polarization  curve.  The  voltage  differences 
at  a  given  current  density  are  proportional  to  the  power  output 
and  the  heat  rejection.  Figure  V-5  shows  that  heat  rejection  due 
to  product  water  evaporation  is  proportional  to  a  voltage  &E 
0.228  volts  regardless  of  current  density,  whereas  the  remainder 
of  the  heat  rejection  varies  with  current  density.  It  is  expected 
that  the  changes  in  the  required  evaporation  rates  will  be  offset 
by  the  vapor  suppression  characteristics  of  the  IEM  and  wicking 
materials. 

4.  5.  5  Experimental  Equipment  for  Evaluating  the  Fuel  System  with 
Air  Recirculation 

Drawings  1076524-536  Sh.l  and  Sh.6  show  the  arrangement  of  the 
principal  parts  of  the  experimental  system  for  demonstrating  the 
concept  of  using  recirculating  air  flow. 

The  equipment  consists  of  closed  loop  air  flow  system  with  the  air 
being  recirculated  by  an  internal  centrifugal  fan.  The  air  flow  rate 
can  be  regulated  by  adjusting  the  voltage  input  to  the  6  volt  d.c. 
fan  motor.  The  fan  can  be  located  at  an  alternate  position  in  the 
system  in  order  to  reverse  the  flow  direction  of  the  air  through 
the  fuel  cell  stack  and  condenser. 

The  fuel  cell  stack  and  the  condenser  are  positioned  in  a  7  in.  I.  D. 
acrylic  tube.  The  recirculation  flow  path  is  through  a  1.5  in.  I.  D. 
acrylic  tube  which  is  used  to  measure  the  recirculating  mass  flow. 

The  fuel  stack  and  ccndenser  units  are  suspended  on  circular 
diaphragms  which  seal  against  the  7  in.  I.  D.  acrylic  tube  where 
necessary  to  prevent  the  recirculating  flow  from  by-passing  either 
unit.  Only  one  sealed  diaphragm  was  needed  for  the  fuel  cell  stack, 
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but  two  sealed  diaphragms  were  necessary  at  the  condenser  to 
prevent  water  vapor  condensation  on  the  coolant  lines  and  the 
exterior  of  the  condenser. 

Considerable  care  was  taken  in  the  design  to  insure  no  uncontrolled 
air  leakage  from  the  system  which  would  make  a  heat  and  mass 
balance  difficult.  Due  to  the  U>w  power  output  level  of  the  fuel  cell 
stacks  of  approximately  9  watts  at  50  amps/sq.  ft.,  the  closed  loop 
system  will  be  enclosed  in  a  temperature  controlled  case  so  that 
heat  leakage  will  not  prevent  a  reasonably  accurate  system  heat 
balance. 

4.  5.  5.  1  Fuel  Cell  Stack  Design 

The  experimental  fuel  cell  stack  consists  of  four  fuel  cells  each 
with  3  in.  x  3  in.  catalyst  areas.  Drawing  1076524-536  Sh.6 
shows  the  arrangement  of  the  cells  in  the  stack.  The  cells  are 
positioned  horizontally  to  facilitate  the  water  transport  from  the 
condenser  through  the  wicking. 

In  order  to  cool  four  cells,  five  air  flow  passages  are  required. 

The  uppermost  flow  passage  provides  air  flow  only  for  cooling  the 
hydrogen  collector  of  the  uppermost  cell.  The  lowest  air  flow 
passage  provides  air  flow  only  for  evaporating  the  product  water 
from  the  lowest  cell.  A  modified  mask  therefore  has  to  be  used 
for  the  lowest  cell.  The  two  innermost  cells  will  be  used  for  close 
comparison  of  performance  since  the  heat  transfer  conditions  for 
these  cells  should  be  reasonably  comparable. 

4.  5.5.  2  Fuel  Cell  Design,  Individual  Assemblies 

The  design  of  individual  fuel  cell  assemblies  is  best  illustrated  in 
Figure  V-2  which  is  actually  distorted  dimensionally  for  pictorial 
purposes.  Drawing  1076524-536  Sh.6  shows  the  details  in  true 
proportion. 

The  hydrogen  is  supplied  through  appropriate  porting  through  an 
epoxy  glass  laminate  frame.  A  closed  hydrogen  chamber  is  formed 
between  a  "hat-shaped"  titanium-palladium  collector  plate  and  the 
IEM.  A  standard  catalyst  and  3  Ta  10-3/0  expanded  metal  tantalum 
screens  are  pressed  into  each  side  of  the  IEM  to  form  the  electrode 
surfaces.  The  IEM  is  separated  from  the  hydrogen  collector  by 
means  of  two  layers  of  screens  to  allow  free  hydrogen  flow  to  all 
areas  of  the  anode  and  provide  a  minimum  distance  between  the  IEM 
and  the  cooled  hydrogen  collector  plate. 
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The  air  side  of  the  IEM  will  be  maintained  at  a  pressure  of  approxi¬ 
mately  20  in.  water  above  the  hydrogen  pressure  in  order  to  make 
satisfactory  electrical  contact  between  the  hydrogen  electrode  and 
the  double  layer  of  screens. 

For  a  simplified  heat  transfer  analysis,  the  principal  heat  generation 
in  the  cell  may  be  assumed  to  be  at  the  air  electrode.  Approximately 
one-third  of  the  heat  rejection  is  by  means  of  the  evaporation  of  the 
product  water  from  the  air  electrode,  as  shown  in  Figure  V-5.  The 
balance  of  the  heat  rejection  in  the  design  illustrated  in  Figure  V-? 
is  by  conduction  through  the  IEM,  the  hydrogen  gas  gap  (approximately 
0.007  in.),  and  the  hydrogen  collector  to  the  wicking  where  additional 
water  us  evaporated  for  cooling.  At  a  current  density  of  50  amps/ 
ft2,  the  temperature  drop  through  the  IEM  is  approximately  0.5  deg.  F 
and  the  drop  through  the  hydrogen  gas  gap  is  approximately  0.05  deg.  F. 

4.  5.  5.  3  Condenser  Design 

A  water-cooled  condenser  is  used  for  simplicity  in  the  experimental 
equipment  rather  than  a  more  complex  gas  to  gas  heat  exchanger. 

Drawing  1076524-535  Sh.6  shows  the  outline  of  one  of  the  coolant 
passages  located  in  each  of  four  condenser  sections.  The  coolant 
passages  are  machined  in  epoxy  glass  laminate  stock  and  titanium 
covers  are  bonded  on  to  seal  the  passages  against  leakage.  The 
condenser  coolant  passages  are  connected  in  series  to  insure  against 
flow  stratification.  The  individual  condenser  sections  are  inter¬ 
connected  by  means  of  "o-ring"  seals  for  easy  disassembly  and 
replacement  of  the  wicking.  The  water  coolant  flow  will  only  have 
a  temperature  rise  of  0.  1  deg.  F  with  a  flow  rate  of  1  gpm  and  a 
fuel  cell  stack  current  density  of  60  amps/ sq.  ft. 

4.  5.  5.  4  Water  Transport  Wicking 

Woven  dacron  wicking  approximately  0.015  in.  thick  is  used  to 
transport  the  water  from  the  condenser  to  the  fuel  cell  stack  as 
shown  in  Drawing  10765240536  Sh.  1  and  Sh.6.  A  continuous  strip 
of  wicking  is  used  to  transport  the.  water  from  the  metal  surface 
of  each  condenser  passage  to  the  corresponding  hydrogen  collector 
plate  of  each  fuel  cell  passage. 


-  127  - 


4.0 


FACTUAL  DATA  (Cont'd) 


Additional  wicking  is  used  to  line  the  epoxy  glass  laminate  surfaces 
of  the  condenser  flow  passages  to  facilitate  the  removal  of  condensed 
water.  The  major  condensing  and  evaporative  surfaces  are  arranged 
on  the  lower  side  of  the  passages  to  minimize  the  effects  of  gravity 
on  the  water  transport  through  the  wicking.  The  "hat-shaped"  design 
of  the  hydrogen  collector,  however,  utilizes  gravity  to  improve  the 
water  transport  through  the  wicking. 

4.  5.  5.  5  IEM  Diffusion  Masks 

As  described  in  Section  4.5.4,  masks  with  graduated  perforations  are 
used  to  improve  the  uniformity  of  evaporation  of  product  water  from 
the  entrance  to  exit  of  the  fuel  cell  stack.  Drawing  1076524-536 
Sh.6  shows  how  the  masks  are  ganged  together  by  supporting  screws. 

The  masks  can  be  adjusted  uniformly  to  change  the  relative  spacing 
of  the  masks  between  the  evaporative  surface  and  IEM  in  order  to 
vary  the  relative  evaporation  rates  from  the  opposing  surfaces. 

The  masks  are  made  of  epoxy  glass  laminate  to  eliminate  the  danger 
of  leakage  currents.  Lip  seals  are  used  on  the  end  of  each  mask  to 
minimize  flow  in  the  vicinity  of  the  IEM. 

4.  5.  5.6  Fuel  Cell  Stack  Current  and  Voltage  Connections 

The  current  collection  on  the  air  side  is  through  the  tantalum 
electrode  screens.  The  estimated  voltage  drop  from  the  center 
of  each  cell  to  the  edge  is  6  millivolts  at  50  amps/ft  ,  The  voitage 
drop  would  produce  a  total  variation  of  current  density  of  approxi¬ 
mately  5  percent  in  each  cell.  It  may  be  possible  to  later  use 
molybdenum  screens  to  further  reduce  the  current  density  gradients. 
This  type  of  current  collection  is  used  since  the  cells  are  small  and 
maximum  freedom  of  space  was  desired  in  the  air  flow  passages  for 
mask  change  and  adjustment.  The  current  collection  on  the  hydrogen 
side  is  through  the  double  layer  of  tantalum  screens  to  the  "hat- shaped" 
hydrogen  current  collector.  The  hydrogen  current  collectors  are 
extended  on  each  side  of  the  cell  assemblies  and  wrapped  around  the 
cell  assembly  spacers  to  make  contact  with  the  screens  of  the 
adjacent  cells. 

Potential  connections  will  be  made  to  each  cell  in  the  stack.  Cell 
potential  and  thermocouple  potential  connections  will  be  made  through 
a  55  pin,  Bendix  hermetic  connector.  The  current  connections  will 
be  made  through  potted  10-32  threaded  posts  in  the  end  plate. 
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4.  5.  5.  7  Mass  Flow  Measurement 

Two  alternate  methods  of  measuring  the  recirculation  mass  flow 
rate  will  be  available. 

One  method  will  be  by  means  of  measuring  the  pressure  drop  for  a 
12  in.  distance  in  a  0.  375  in.  I.  D.  tube.  A  flow  rate  of  1. 5  cfm 
will  produce  a  pressure  drop  of  approximately  0.  1  in.  water  in 
this  flow  meter  assuming  the  flow  is  essentially  laminar. 

The  alternate  device  illustrated  in  Drawing  1076524-536  Sh.6  is  a 
Thomas  type  meter  in  which  the  temperature  rise  of  a  heated  stream 
is  used  to  determine  the  flow  rate  thermodynamically.  The  Thomas 
meter  will  be  constructed  of  three  planes  of  resistance  wire  grids. 

The  center  grid  is  to  be  heated  electrically  and  the  additional  grids 
are  used  as  resistance  thermometers. 

The  recirculation  mass  flow  meters  are  to  be  calibrated  prior  to 
use  by  means  of  calibrated  rotameters.  Rotameters  could  not  be 
used  satisfactorily  for  measuring  the  recirculating  flow  rate  due  to 
their  relatively  high  pressure  drop. 

The  flow  rate  of  air  into  the  system  is  to  be  measured  by  a  Fischer 
and  Porter  Tri-Flat  flowmeter  with  a  0.062  in.  dia.  spherical  float. 

4.  5.  6  Conclusions 

The  analyses  of  the  air  operated  IEM  fuel  cell  system  with  recirculation 
indicate  system  feasibility  and  that  a  high  degree  of  uniformity  of 
operating  conditions  can  be  achieved  in  this  system.  More  definite 
conclusions  can  be  made  after  the  experimental  results  are  obtained 
in  the  next  phase  of  the  work. 
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APPENDIX  1-1 

MICROGRAPHS  FROM  ADHESION  TEST 


This  section  comprises  additional  micrographs  of  "scratches" 
made  with  the  adhesion-te sting  device  of  Photograph  I  -  1,  and 
are  described  in  detail  in  Section  4.  1.  3.  3. 
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APPENDIX  1-2 


THE  SHARPE  AND  SCHONHORN  ADHESION  THEORY 


According  to  Sharpe  and  Schonhorn  (21),  there  exists  a  one-to-one 
corr espondence  between  thermodynamic  spreading  and  strong 
adhesion.  The  former  concept  is  usually  expressed  as  requiring 
the  spreading  coefficient,  Ssp  to  be  positive  or  zero  (22).  Since: 


sp 


Y„  -  I Y 


L2V2 


V2 


(1-12) 


the  liquid,  L will  only  spread  on  the  solid,  Sj»  if: 


> 

Yc  =  Y,  y 
SI  L2V2 


+  Y 


S1L2 


(1-13) 


i.e.  if  the  sum  of  the  surface  free  energies  of  the  liquid,  contact 

with  its  saturated  vapor),  and  that  of  the  solid/liquid  interface  is 
less  than  that  of  the  solid. 


This  can  be  shown,  from  the  Young-Dupre  equation  (22),  to  be 
thermodynamically  equivalent  to  having  the  work  of  adhesion, 

W adh»  equalling  or  exceeding  the  work  of  cohesion  of  the  spreading 
substance,  i.e.: 


W 


adh 


=  W 


coh 


(1-14) 


so  that  strong  adhesion  can  only  occur  when  the  contact  angle  is  zero 
and  the  liquid  can  spread  on  (or  "wet”)  the  solid  (21).  Since  neither 
Ygj  ,  nor  vs  L2  can  bc  directly  determined  in  any  satisfactory 
manner,  the  empirical,  but  operationally-uscful  concept  of  c ritical- 
surface  tensions  (19)  was  adapted  by  these  authors  in  the  form: 


Ms,) 


>  Y 


L2V2 


(1-15) 


along  with  the  further  assumption  of  approximate  equality  between 
y  and  YLy  for  the  solid  and  fluid  forms,  respectively,  of  the  same 
material. 
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LIMITING  CURRENTS  FOR  VARIOUS  DIFFUSIONAL  PROCESSES 

In  order  to  restrict  possibilities  of  rate-controlled  gas  diffusional 
processes,  various  differential  equations  have  been  solved  and 
limiting  current  calculated  for  defined  systems  and  for  assumed 
boundary  conditions. 

1.  Cylindrical  Diffusion:  Linear  diffusion  of  02  in  the  O2/N2/H2O 
system  for  electrodes  exposed  to  semi-infinite  space,  (Figure  III-  IB  -  1) 
does  not  adequately  represent  the  process,  since  radial  diffusion  can 
also  occur.  For  the  following  boundary  conditions: 


c  (r,  ®  )  =  C° 

C  (r,  o  )  =0 

C  (b,  2)  =  C° 

the  solution  of  Laplace's 
eqn.  : 


dfc 

dr2 


_1_  d£  +  d2c  _  q 


dr  dr 


C  =  c”  1  -  2 


becomes 

-2-  r/fc) 
o  m 


■  ■  '■  \  J .  (k  ) 

m- 1  m  1  m 


(III  -  20) 


where  C  =  partial  pre  ssure,  i.  e.  p^ 
C°  =  bulk  partial  pressure, 

X.  =  roots  of  J  (\  ) 

m  o'  m' 


2  . 


The  local  current  density  becomes: 

.0  S —  J  (x  r/b) 


j(r)  =  nFD 


dc 

dz 
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z  =  0 
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m  =  l 
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( III -21) 


and  the  total  current: 

,b  c  V”  1 

I  =  Z-rr  I  rj(r)dr  =  4imFDC  Z _  £ 

m=  1  m 
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The  previous  result  yield  current  dependence  on  electrode  diameter 
rather  than  surface.  In  fact,  the  boundary  conditions  suggest  that  a 
large  contribution  to  electrode  current  is  due  to  radial  diffusion.  The 
thickness  of  the  diffusion  layer  as  determined  from: 


I  =  nFDA 


=  nFDA 

z  =0 


(111-23) 


becomes: 


b/4 


(III -24) 


The  total  current  I  =  4000  mA  for  an  electrode  surface  of  0.  7  cm2 
(D  =  0.  25  cm2/sec.  ,  n  =  4,  C°  =  0.  9  10*5)  and  for  air,  whereas  the 


measured  value  is  100  mA. 


2,  Cylindrical  Diffusion  In  A  Pore:  Analogous  to  case  1,  would  be 
the  cylindrical  diffusion  in  a  pore.  However,  the  B.  C.  now  becomes: 


where  the  condition  Co  =  O  becomes  a 


C  (r,  h)  =  C° 

C  (r.o)  =  C° 

C  (b,  z)  =  O 

and  the  following  concentration 

distribution  can  be  obtained: 

C°-C°  cosh(\  h/b) 

.  o _  m 

~  +  sinh  (\  m  h/b) 

(III  -  2 r ) 

special  case  of  eqn.  Ill  -  2  5 . 


nh(\  m  z/b)' 
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The  local  current  density  becomes: 


i(z) 


r=b 


2nFD 
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H  cosh(Xmz/b)  + 

m=  1 


sinh(\  m  z/b) 

o 

and  the  total  pore  current  becomes,  for  the  case  Co 


C°-  C°  cosh(\  m  h/b) 
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Y“  cosh(\  h/b)?l 
)  m 

I  (pore)  =  4rnFD  bC  —  Sh(Km  h/b) 

m=  1  rn  m 


(III  -  2  T) 


where  represents  an  effective  diffusion  coefficient.  Thus  the  pore 
current  will  depend  on  pore  diameter  and  on  a  function  of  h/b.  If  the 
dimensionless  group  (*-mhA>)  is  large,  eqr..  III-27  can  be  reduced  to 
eqn.  111-22  derived  abo'  e. 

3 .  Diffusion  trrougn  a  liquid  layer  including  particles. 

The  rate  for  cathodic  oxygen  reduction,  according  to: 

Oz  -  2  H20  *  4e  - ►  4  0(4. 


<  n  be  expressed  for  large  ic  v-lues  as 

,  .  iHl  E 

:  =  nFACk  (rl2 0)“  e  RT 


(III  -  28) 


lr.  fa  t  in  acidic  media  for  pressures  up  to  70  atm.  ,  Krasilschchiko  v  (35) 
and  Wir.kelnv. n  (36)  have  shown  ‘-a*: 

.  Ailfil  r 
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The  two  (  qr. .  h.omf  idt  r.ti  il  il  k'kj  =  qFA  (H2O,  and  if 

pn  ^  1  .Q  =  0  5  as  determined  by  these  authors.  Eqn.  111-28  can  be 

rewritten,  if  related  to  equilibrium  conditions,  as. 

anF 


=  nFA:'  /C'Ce 

V  c'  c, 


RT 


(III- 30) 


w*’ere  C  =  actual  02  tor  . 

Co  =  equilibrium  07  '  on:  . 


for  a  or.e -dimensional  analysis 


(III  -  3 1) 

De  -  effective  diffusion  coefficient 
i  =  current/ area 

k  -volume/area  =  cm,  represents 
a  property  of  the  medium,  as 
already  considered  by 
C.  Wagner  (40) 
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an  F 

If:  2  1  1  *e  RT 

B  De  K  C0  6 

eqn.  Ill- 31  becomes: 

2 

1-4  -  B2C  =  0 

dx*" 


(III -32) 


(III -33) 


Solution  of  eqn.  Ill  -  33  and  subsequent  differentiation  to  determine 
rates,  yields: 

(III -34 

I  =  nFDe  SCQ  B  tanh  (Bd) 
where  S  represents  the  transfer  area. 


If  Bd  »  0, 


I>S  =  nF  (D  C  i  /K) 


(III- 35) 


suggesting  a  current-dependent  e  on  De  1/Zand  on  the  electrode 
property,  i.  e.  thickness  of  the  film  or  structure  of  the  electrode 
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APPENDIX  III- 2 

DETERMINATION  OF  CONCENTRATION  DISTRIBUTION  IN  CHANNELS 


^Ci.  =0 

lCi-1 


(III- 36) 
(III- 37 ) 


Ci  represents  partial  pressures  in  the  three  component  systems. 
Boundary  conditions  are: 


yG 


Ci  (o,*)  -  c\  j  e.c.i 
CL  (u,y)-CL J 


cQi  (*,  o)  -  0 
Cfj]  (X,  o)5''^t)(6,C.  2 
Cw  (*,<>)'•  Ht) 


•*—  d.  C( 3 


oL  Ca/, 


2b 


Electrode 


7  7 


^ 


n 


rt° 
.’•-fV 


B.C.3 

6.C.4 


The  concentration  of  component  (i)  is  related  to  the  partial  pressure  by. 

£■  -  — t-  •  —  (III- 38) 

V  RT  1 


Solution  of  equation  III  -  36  are,  for  the  different 


v  c 


I  -  1  4/7  SlO  c*„  X  (cesho<„ £  -  than  ^ 


components: 
-*>)]  ,m-39) 


C\i  •Cl>  *  I 
Z  rt.i 

/$(*)••(  (r)  c.osh  ttn 

\ 


cosh 

+  ^)x  S*»«„ 


(III  -  40 ) 


-  X  X)/)  5//)  ^  x 


n</ 


(III-  41) 
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with 


o+n 


njl 

^ n  ~  zb 


Local  current  densities  can  now  be  obtained  from: 


jd)r  -nFD 


1 


rl 


and  are,  (for  O2  diffusion) 

FDC'  I  <*„  S<o<xn  X  rh<*n  £ 


(III-42) 


(III -43) 


and  for  H^O  diffusion 


j.  (t  •-x'PD  ^-<*n  An  °*/)  X  s/> 


P  ’ 


(III  -  44) 


Average  channel  current  densities  are: 

<r  f (•On-iT  M  ,  Xo  \ 

brr  £,  ~ —  &  Tb) 


(III- 45) 


and 


n'FP 
H,0  -  t>ir  * 


nil 


71 


}(T)-f(r)c,ihoin 


p 


n  3t 


(III  -  46) 
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APPENDIX  III -3 

DERIVATION  OF  EQUATIONS  FOR  FORCED  FLOW  CONDITIONS 
IN  r.HANNELS  WITH  RECTANGULAR  CROSS-SECTION.  DETER- 
MINATION  OF  LIMITING  CURRENT  DENSITIES. 


From  a  general  continuity  equation  for  transport  processes. 

VC/C  +  W  <±-  -!>>(&+&) 

V^dy  -  +dytJ 


(III  -  47) 


based  on  channel  geometry  described  in  Figure  III-23  and  taking  into 
account  simplifying  assumptions,  i.  e. 


i/  dc  _  dzC  _  n 
vy  dy  -  d&~° 


(III  -  48 ) 


a  general  equation  can  be  obtained  to  describe  diffusional  processes  in 
the  channel: 

vxf  .  D  djc 

xc/x  dy 2  (hi  -  49) 

where  is  only  y-dependent,  and  has  to  satisfy  the  B.  C.  (Figure  III-23) 

Vy=0=  -0  ( 


Yy. 


=  V. 


-/:y0/z  ~  Vmax 

D  represents  the  diffusion  coeffi  ier.t  of  the  binary  system  N2/02,  the 
concentration  c  is  related  to  the  partial  pressure  of  oxygen  by  means  of 


p  =  C  RT 


(III  -  50) 


Analytical  solutions  of  eqn.  III-49  are  possible,  in  principle,  providing 
Vx  is  a  simple  function  of  y.  Different  assumptions  regarding  such 
velocity  profile  under  laminar  flow  conditions  will  be  considered. 


1.  Parabolic  Velocity  Profile  represented  by: 

v=ay(y0-y) 

and  satisfying  B.  C.  1,  yields: 


V 


MAX 


adZ/4 


(III- 5 1) 


(III- 52) 
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Conservation  of  volumetric  flow  rates  (nearly  correct  for  already  low 
flow  rates): 


Q  ~  Z'  V  AVE  ds 


allows  to  calculate 

VavT  =  3 

yielding  the  modified  eqn. : 


6  Vm  w(w-y)dc 

Dyl  y dy 2 


(III- 53) 


(III -54) 


(III- 55) 


!  58 


The  change  of  variable  y/yo  =  Y/  and  adoption  of  the  dimensionlei 
group: 

DVaVE  4  &H  ^AVE  (III- 56) 


Re  Sc 


£ 


£> 


where  the  hydraulic  radius 


RH-*2b. 

W*y0 


Vy  =  channel  width, 

V  =  channel  height 
'JQ 


(III  -  57) 


yields  simplified  differential  equation  III-59. 

If  V  «  V/  eqn.  Ill  -  5  5  becomes: 

/,  ydc  _  /  dk  mi.58) 

y/'  d,)d*  3y0(ReSc)  dy* 

Now  with  the  change  of  variable:  X,  =  ^/^/Qt  and  defining  3  ReSc  =  N; 

..  //-w  )dc  =  l  dfC  (HI - 59) 

v  Wdx,  n  dyf 

with  the  following  B.  C. :  ,  ^  ^ 

c  =0  at  y,  - 0  anc/  X;  >0  i  r  r 
C  =c°atXi<  0  and  onyy,  j  b‘C\ 

Eqn.  Ill  -  59  does  not  yield  an  analytical  solution.  However,  numerical 
solutions  by  iteration  techniques  would  be  possible.  No  attempt  has  been 
made  to  apply  such  techniques. 
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Simpler  differential  equations  can  be  derived  if  non-parabolic 
velocity  profiles  are  assumed. 

2.  Angular  Velocity  Profile 

Close  approximation  to  parabolic  profiles  are  obtained  from: 


v  ■  Vmax  0 -y,2) 


(III  -  60) 


(see  figure  III  -  2 3) 

where:  VAVf 


V 


MAX 


(III  -  61) 


For  X,  =  X/y  and  ^  (ReSc)  =  N,  an  eqn.  similar  to  eqn.  HI  -  58 

1  Jo  /6 


can  be  derived: 


A  2\<k  =1 
V  y>  '  d'K,  N  dyf 


(III  -  6 1) 


with  the  following  B.  C. 


C  «0  aty,  -  +1  and  X ,  >*0l  g  ^  ^ 

C  ~Ccat  x,<0  and  any  yj 

By  rearrangement  of  variables,  the  dimensionless  group  N  can  be 
eliminated  to  yield  ? 

o-y,2)dA  -  $ 

This  differential  eqn.  allows  to  separate  variables  to  yield: 


X'_  Y 

x  (Ffjr 

requiring  the  solution  of: 


=  -  A 


y'Va(/-^)y-o 


(III  -  6  3) 


(III -64) 


with:  Y (/)  =  Y O)~0,  which  could  be  appreciably  simplified  if 
y /+/)=  Y"(-l )  “0  8Uggesting  a  Plane  of  symetry  at^  =0.  Practically, 
this  would  suggest  to  determine  limiting  current  densities  for  two 
parallel  flat  electrodes.  If  now, 
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it  is  possible  to  show  that  equation  III -66: 


is  a  Whittaker-type  differential  equation  (41): 

dli.fL+'H  +  lL 

■  ■  ’+'~4  +  4tjj  /(,  (f)2 


where 


d 

(p  -  -Tat 


(III- 65) 

)z=o 

(III- 66) 

)Zr=0 

(III -6  7) 

Solutions  of  this  eqn.  III-67  are  available.  Additional  work  will  be 
required  to  show  applicability  to  reported  experimental  results. 

3.  Linear  Velocity  Profile 

Assumption  of  linear  profiles  yield,  in  essence,  simpler  differential 
equations  and  as  a  result  have  been  the  subject  of  investigations  on 
heat  and  mass  transfer  processes.  If,  as  represented  in  figure  III-23, 

>,  0  (III - 68) 


V=  Ky  for  y 


the  following  eqn.  can  be  derived: 

dc=  D  j_<A 

d*  Ky2o  y;  rfyf 


(III-  69) 


with: 

C  -  0  at  y,  =  0  and  any  X 
C *C °at  x  ^-0  and 
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Development  of  eqn.  111-69  depends  on  the  velocity  gradient  K.  For 
laminar  flow:  „  $6 

JU.  (111-70) 


$k-o  '  K 


where: 


H¥) 


(III-71 ) 


d  =  film  thickness;  Q  =  flow  rate,  yU.  =  kinematic  viscosity, 
^  -  acceleration  in  the  x-direction 

For  eqn.  Ill  -  60  f  o 

D_  =z6i  J_ 

Kf  3  ReSc  fa 

and,'^^;  °D 

/<y2  /?e5c 

I/O 

For  small  flow  rates:  OC.  —  x2 


(III  -  72) 
(111-73) 


and  eqn.  Ill  -  69  becomes: 


dc  _  ]_  J_  d_c 
dx  N  y,dyf 


(III- 74) 


For  the  change  of  variables:  C  mC  ^Qeqn.  HI  -  74  becomes: 

dc*=  i  J_  d2C*  (III-75) 

dx  N  y,  c/y(2 

with  the  B.  C.  : 

C*=  0  /or  x ^ 0  Gmdy,>'0 
C*  =  Co  for  X  >  0  tfnc/  y,  *  0 

This  type  of  equation  has  been  studied  by  Kramers  &  Kreyger  (42). 

The  rate  of  transport  at  the  interface  can  be  shown  to  be: 

X\  _  /->  \  /~i 


dc_)  _  ~_3 


Co 


3x 


(III-76) 


dy,  hpvQ/s) 
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In  terms  of  original  variables,  eqn.  111-76  becomes: 

fdc\  _  ~_3_  Cfl. /_3X  \'/3 
[dyly.o  r(l/3)  yo  \2ReSc/ 

The  position-dependent  current  can  be  derived  as: 


(111-77) 


j  =-rj  rnf — )  =3gEQ-  Qs_(  ) 
^  r)hU\dy/y=o  r(f/3)  y°  \ZReSc) 


-hs 


(111-78) 


In  terms  of  the  experimental  variables  for  operation  on  air,  at  constant 

T),D,Co;  2/  /  /Q  \^5 

2 L(x)=2r*FD  3W  3C0{^j) 


(111-79) 


Eqn.  Ill  -  7  9  does  not  describe  limiting  currents  at  X  =Q  If  the  function 
is  normalized  by  taking  ij^  (o)  =  130  mA/cm*'  (as  determined  from 
limiting  current  density  measurements  for  flow  and  non-flow  conditions 
of  air)  linear  relationship  should  be  obtained  between  ^ 

toy  [  1 1  (y)/iL  (o)]  and  (Q/x y/)  5 

Constant  flux  at  the  Electrode 

In  the  case  of  galvanostatic  techniques,  constant  and  uniform  current 
distribution  are  obtained,  corresponding  to  (oC/(jy)y~0  =  Const- 
Under  such  conditions,  equation  111-75  can  be  rewritten: 

dzC*_  j  d 
dydx  N  dy, 

and  if  (dcy/dy )  =  y/ 

djp=  l 
dx  N 


y> 


1  f  L  4l 

U  dy,  * 


(III -80) 


(III -8 1) 
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Eqn.  Ill -81  has  been  studied  in  the  case  of  heat  transfer  for  non-Newtonian 
liquids  for  laminar  flow  (43).  Now,  the  B.  C.  becomes: 


x  co  ■  y-o 
y,  =  0:  (p*=l 
y=oo  ip-0 


providing  (j)  — 


_J, _ 

(9*/n)vi 


(dc*/d^)jj 

(dcydy,)y=0 

(III -82) 


it  can  be  shown  that: 


7 T<fi*'+(3Tr-l)  (/)*'=  0 


(III  —  8  3) 


77=0  :  <!>*•! 

B.  C.  * 

]T  =  60  :  (p  =0 


Eqn.  Ill  -  8  3  admits  the  solution: 


/- 


T'tt/Z) 


•ao  3 

7re"  dr 


(III  -84) 


And 


-c 


* 


% 


N, 


^  re-r 


rWs) 


-77  )  /- 


r(*/3  El) 

r(2/3) 


(III  -  86) 


For  y  i  =  Oj 


(III  -87) 


r* 

-  c  WALL 

% 


\^5 

3x)  6  I 

nJ  T(z7$) 
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* 


f-  (c0-cWALL)(^)  r(2/£) 

For  limiting  conditions,  when  C  wall  =  0,  eqn.  HI-89: 

%r  r 

can  be  compared  with  eqn.  III-76.  If  relations  of  the  type: 

In  (/-  YiL  )=o<(E  -e)  = 


(III  -  88 ) 


(III  -  89) 


(III-90) 


are  descriptive  of  the  polarization  curves,  and  in  fact  represent  In 
(Cw/C0  ):  the  overall  curve  can  be  established  if  (X  and  ij^  are  known. 
However,  since  CX  has  physical  significance,  values  of  iL  (x)  as  defined 
by  eqn.  Ill- 79  are,  in  principle,  sufficient  to  describe  the  current  - 
voltage  relationships. 

In  addition,  wall  concentration  of  oxygen  can  be  calculated  as  a  function 
of  concentration  gradients  near  the  electrode,  i.e.  local  currents. 
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APPENDIX  IV 

IV- 1  MASS  TRANSFER  ANALYSIS 

Analysis  of  the  Possible  Rate  Limiting  Mechanisms  in  the  Diffusion 
of  Oxygen  from  an  Air  Stream  to  the  Surface  of  the  Oxygen  Electrode 
in  a  G.  E.  "GPS"  Type  Ion  Exchange  Membrane  Fuel  Cell 

The  rate  at  which  maaa  may  be  transferred  in  a  fluid  is  determined  by 
the  physical  properties  of  the  fluid,  its  state  of  motion,  and  the 
driving  potential. (44,  45,57)  When  the  fluid  in  which  mass  transfer  is 
occurring  is  stagnant  or  is  in  laminar  motion  normal  to  the  direction  of 
mass  transfer,  the  mechanism  of  transfer  is  that  of  molecular  diffusion 
which  has  been  defined  as  the  spontaneous  intermingling  of  miscible 
fluids  placed  in  mutual  contact  accomplished  without  the  aid  of  mech¬ 
anical  stirring.  The  shear  stresses  in  fluids  in  laminar  motion 
(transfer  of  momentum)  exist  through  the  action  of  viscosity.  At  the 
interface,  the  fluid  velocity  relative  to  the  solid  surface  is  substantially 
zero  while  at  Increasing  distances  from  the  interface  the  velocity  pro¬ 
gressively  increases.  In  effect,  the  fluid  may  be  considered  as  composed 
of  thin  layers  which  follow  the  contours  of  the  surface  without  mixing 
with  the  adjacent  layers  except  as  mingling  results  from  molecular 
diffusion. 

Even  under  the  conditions  of  turbulent  flow  there  is  always  a  layer  of 
fluid  at  the  Interface  in  w’.ich  laminar  flow  is  maintained.  The  thick¬ 
ness  of  this  laminar  layer  or  film  depends  upon  the  conditions  of  flow 
and  becomes  less  as  the  velocity  is  increased.  At  the  low  Reynold's 
numbers  (10  to  150)  observed  in  this  experimental  investigation,  the 
laminar  layer  includes  the  entire  stream. 

On  the  basis  of  the  film  concept  mass  transfer  is  considered  as 
resulting  from  molecular  diffusion  across  the  effective  film  under  the 
influence  of  a  concentration  gradient. 

1  -A.  Countercurrent  Diffusion  of  Oxygen  from  an  Air  Stream  to  the 
Oxygen  Electrode  and  Water  Vapor  from  the  Oxygen  Electrode 
to  the  Air  Stream  in  the  Presence  of  an  Electrochemical 
Reaction.  (45) 

In  order  to  determine  whether  the  diffusion  of  oxygen  from  an  air  stream 
to  the  oxygen  electrode  is  limited  by  a  stagnant  gas  film,  or  by  the 
countercurrent  diffusion  of  water  vapor  into  the  air  stream,  or  by  gas 
diffusion  into  pores,  or  by  a  liquid  film  in  the  electrode  an  analysis  of 
these  mechanisms  was  made  as  follows. 
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A-  1.  Diffusion  in  Complex  Systems  -  Gas  -  Gas 

"In  heterogeneous  chemical  reactions  two  or  more  components  may 
diffuse  to  and  from  an  interface  in  the  presence  or  absence  of  a  film 
of  inert  gas.  The  transfer  may  be  equimolal  in  both  directions  or 
may  correspond  to  an  increase  or  decrease  in  the  number  of 
moles.  "  (45) 

Quoted  from  (45): 

"For  the  general  gaseous  reaction  a  A  +  bB  rR  +  sS  proceeding  in 
the  presence  of  an  inert  gas  I  the  following  equation  which  includes 
the  effect  of  all  components  present  applies.  Thus  for  the  diffusion 
of  A, 


-RT 


i£A  =_Jpa_  cAcB(uA-uB)  **arcacr(ua-ur)  (IV- 1) 

dL  dL  Ab 

**AS  CACS  (“A'  US»  t0*AI  CAC1<UA> 

where  -  dPA  =  partial  pressure  gradient  o(  component  A  in  the 
dD  direction  of  diffusion  L 

°*AB,  AR,  AS,  AI  =  proportionality  factors 

C  -  molal  concentrations  of  A,  B,  R,  S,  and  I,  resp. 

A,  B,  R,  S,  I 

Ua  q  r  S  ”  linear  velocities  of  diffusion  of  A,  B,  R,  S  in  the 
’  '  direction  of  diffusion  of  A 


RT 


An  alternate  equation  is, 


JL _  -0$  =  - L_  <rAaPB'rBaPA)  +  -i— (rAaPR'rRapA) 


dL 


D 


AB 


D 


AR 


Das 
where  Tf 
R 


—  (rAaPS-rSaPA>  +  p—  rAaPl 


AI 


=  total  pressure  (atm) 

f  cm  -atm 

=  gas  constant  ^82.  06  g  mole.oK 


) 
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4.  0 


DAB,  AR,  AS,  AI 
rAa,  Ba,  Ra,  Sa 


PA,  B,  S,  I 
L 


=  absolute  temperature  (°K)^ 

=  diffusion  coefficients  (  ^—/characteristic  of  the 

v  sec. / 

pair  of  components  AB,  AR,  AS,  and  AI 
=  the  average  molal  rate  of  diffusion  of  A,  B,  R,  and  S 

(g  -  moles  \ 
cm^-sec.  J 

=  partial  pressure  of  A,  B,  S,  and  I  (atm) 

=  distance  in  the  direction  of  diffusion  of  A  (CTn) 


The  terms  Uj  and  rj  do  not  appear  since  there  is  no  net  diffusion  of 
inert  gas. 

Because  of  the  uncertainty  of  the  proportionality  factors  and  the 
corresponding  diffusivitie s  in  mixtures  of  several  components  and 
because  of  the  complications  resulting  in  attempting  to  apply  different 
values  for  each  binary  pair,  an  average  diffusivity  DAm  will  be  used 
to  replace  the  separate  values  recorded  in  Equation  IV-2.  The 
average  value  of  the  diffusion  coefficient  will  be  taken  as  the  weighted 
mean  of  the  values  for  each  pair.  Thus  for  component  A, 

U-Vda™  =nbdab  *  nrdar  *  nsdas  *  nida.  "v-3) 

where  NA,  Ng.  . .  =  average  mole  fractions  of  components  A,  B.  .  . 
in  the  diffusional  film. 

From  the  stoichiometry  of  the  general  reaction, 


=  (^a,-r,  =  -(-4-J 


rAa’  r  Sa 


(I  V  -  4) 


Aa 


Ba 

Substituting  Equation  IV-4  in  IV-2  and  replacing  each  diffusion  coeffi¬ 
cient  with  the  average  value  DAmgives  _  (IV-5) 

_T£_  dPA  _  1 


RT 

or,  since  'Tf: 


dL 


D 


Am 


3  (IV-5) 

rAa  PB  +  PR  +  PS  +  pI*PA«!LLp_j 


pA  +  pB  +  PR  +  PS  +  PI 
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r  A  dL  =  - 
Aa 


RT 


Am 


dPA 


/a+b  -  r  •  s  \ 
■"-»a( - a - ) 

I_  /  dpA  \ 

\  ir+  pa'a  J 


(IV- 6) 


whe 


re  6 


RT 

r  +  s  -  a  -  b 


A  " 


Integrating  between  the  film  boundaries  L  (bulk)  and  L  (interface) 
corresponding  to  an  effective  film  thickness  Bq  gives 


rAa  - 


PAm*  /  1  ) 

RT  Bc  i  6A  / 


In  /  IT  +  PA  6  A 
W+  PAi*A 


(IV-7) 


Equation  (IV- 7)  also  may  be  written  in  the  form 


IT 

Am 


(PA  '  PAi> 


Aa  RTB  Gpf 

(•ff^A^)-  (^VAi) 


(IV-8) 


where 


Pf  = 


In 


fTT+<SA  p a 

U+  «A  PAi 


(I  V- 9) 


Bq  =  effective  film  thickness  (  m) 

P£  =  Film  Factor  (atm) 

Equation  IV-9  is  the  logarithmic  mean  value  of| 
boundary  limits  of  the  gas  film.  Thus,  for  diffusion  of  A, 


(Wa*) 


over  the 


^VaI 


lm 


i 

1 


(IV- 10) 
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Similarly,  for  S 


P 


f 


ir-6s 


,  r  _  r  +  s+...-a-b-... 

where  O  e  - 

o  s 

(  TT  +  6a  Pa 

Where  the  ratio  of  l — — - 

\  Tf+  6a  pAi 


(IV- 11) 


is  small, 


for  example  less  than  1.2,  the  arithmetic  mean  is  in  close  agreement 
with  the  logarithmic  mean  and  may  be  used  for  most  purposes. 


For  equimolal  diffusion  it  may  be  recognized  from  Equation 
that  pf  =Tf  and  Equation  (IV- 8)  reduces  to  the  form 


rAa 


rBa 


-  D 


AB 


RT 


d  PA 
L 


(IV- 10) 
(IV-  12) 


which  is  the  fundamental  equation  for  diffusion  of  gases  but  is  directly 
applicable  only  to  the  case  of  equimolal  countercurrent  diffusion  of  a 
binary  mixture  of  A  and  B. 

For  the  diffusion  of  only  one  component  in  the  presence  of  a  stagnant 
gas,  £ A  in  Equation  IV-9  is  equal  to  -1,  and  the  value  of  p^  becomes 
that  given  by  the  following  equation 


p  =  logarithmic  mean  of  p  and  p 
f  B2  Bj 

Subscript  1  =  refers  to  the  bulk  value 
Subscript  2  =  refers  to  the  interface  value 
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Equation  (IV-8)  and  (IV-9)  also  can  be  applied  where  multicomponent 
diffusion  takes  place  without  stoichiometric  relations.  For  example, 
in  the  countercurrent  adsorption  of  one  gas  and  desorption  of  another 
the  constants  a  and  s  may  be.  taken  as  the  actual  molal  rates  o^f 
adsorption  and  desorption,  respectively,  and  the  ratio  - — - 

may  assume  any  vaiue. 

Equation  (IV-  1)  indicates  that  for  a  given  concentration  gradient,  the 
rate  of  diffusion  of  A  is  favored  by  the  presence,  of  another  component 
diffusing  in  the  same  direction  and  retarded  by  another  component 
diffusing  in  the  opposite  direction.  Approximate  values  of  a  and  s 
may  be  obtained  in  such  a  case  by  calculating  the  rate  at  which  each 
component  would  diffuse  under  Vie  conditions  of  the  system  if  the  other 
were  absent.  "  End  of  quote  cf  (45!.  These  principles  will  be  applied  to 
fuel  cell  systems. 

1 A  -  2  Bulk  Diffusion  -  Gas  Film  Controlling  (45). 


System: 


Air  Composition  In 


1.150  cm 


Gas  Composition  Out 


21.0%  02  by  Vol. 
79.  0%  N2  by  Vol. 

=  1.04  atm 
T  =  296.  5°K 
Air  Feed  Rate  = 

5.  32  x  1C4  g-mole 


/  sec. 


////////W 


18.  2%  02  by  Vol. 

77.  2%  N2  by  Vol. 

4.  6%  HzO  by  Vol. 

=  1.  00  atm 
T  =  314.  7°K 
Gas  Film  -  Bq 

•Catalyst  Electrode  at 
319.  3°K 
E.  M. 

Gas  Exit  Rate  =  5.45x  10"4  e-mole/sec. 


Basis:  (1)  G.  E.  I.  E.  M.  Fuel  Cell 

(2)  Current  =  4,  85  amps 

(3)  Current  Density  =  86.  7  ma/cm2  (80.  7  amps/ft2) 

(4)  Geometric.  Catalyst  Area  =  56.  0  cm2  (0.  0602  ft2) 

(5)  Product  Water  Evaporation  Rate  =  2.  51  x  lO"5  g-mole/sec. 

=  100%  at  Indicated  Air  Flow  which  is  8.9  times  the  stoichiometric 
requirement 

(6)  Oxygen  Consumption' Rate  =  1.  254  x  10  ^  g  -mole/sec. 

(7)  Air  Gap  =  0.  1066  c.m 
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21  0+18  2 

%  o, _  =  avg.  Oxygen  Percentage  in  Bulk  (1)  =  - - — - - —  =19.  6% 

L,  3.Vg.  ^  — — 


A  linear  gradient  of  oxygen  concentration  from  the  inlet  to  the  outlet  of 
the  reaction  zone  is  justified  because  of  the  shortness  of  the  path. 


IT 


avg. 


„  .  „  /  1.04  +  1.00  \  , 

=  Avg.  Total  Pressure  =  I - - - )  ~  == 


02  atm 


>2(1) 


Avg.  Oxygen  Partial  Pressure  in  Bulk  =  0.  196  (1.02)  -  0,  200  atm 


p  =  Avg,  Oxygen  Partial  Pressure  at  interface  -  0^  atm 

°2{2) 

The  oxygen  gas  partial  pressure  at  the  catalyst  interface  can  be 
assumed  to  be  negligible  compared  to  the  bulk  oxygen  partial  pressure; 
therefore,  in  order  to  approximate  the  oxygen  partial  pressure  gradient, 
p  is  set  equal  to  zero. 

°2(2) 

p  .  .  =  Avg.  water  vapor  partial  pressure  in  Bulk  = 

H2°U)  [q:+  0.  70(0.0796) 


=  0,  0278  atm. 


Linear  gradient  taken  for  same  reason  as  for  %  C>2  avg. 

p  =  Avg.  water  vapor  partial  pressure  at  interface  =  0.  101  atm 

H  20  ( 2 ) 

lA-2a.  Counterdiffusion  of  Oxygen  and  Water  Vapor  in  Presence  of 
Inert  (45,  46,  47,  48,  56)  Nitrogen  Gas  -  Bulk  Diffusion 
Bulk  Diffusion  Coefficients 


sec. 


1.  75 

.  I- 

-  n 

(  T 

\  /  1.  00  V 

=  0.181  312-4\ 

n2  do2-n2 

l  To 

1  \  1.02  / 

\  273  ) 

312.  4°K  273°K 

1.  02  atm  1. 00  atm 


Do  -H  o  =°-2,0(-p5f)=  0-284  cm2,"c- 


2  2 

312.  4°K  312.  4°K 

1.02  atm  1.  00  atm 
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D„  ~  M  =  0.266  =  0.  260  cm2/sec. 

h2u  "  w2  1.  02 

312.  4°K  312.  4°K 

1.  02  atm  1.  00  atm 

Calculation  of  the  Arithmetic  Mean  Composition  of  the  Gas  Film 


At  Interface  -  =1.02  atm 


p„  =0  atm 

°2(2) 

p.,  =  0.  919  atm 

N2{2) 


p  =0.  101  atm 

H2°(2) 


Mean  Mole  Fraction  of  in  Gas  Film 


In  Bulk  -  If  =  1.  02  atm 


P-.  =  0.  200  atm 

°2(1) 

p.,  =  0.  792  atm 

N2(l) 


p  =  0.  0278  atm 

H2°(l) 


_^0  +  0.20C^/  { 


02  =  0. 0980 


Mean  Mole  Fraction  of  in  Gas  Film  =  1. 02  -  0.  839 

^  A).  101  +  0.0278\/ 

Mean  Mole  Fraction  of  H.,0  vapor  in  Gas  Film  =  ^ - - - - - y  1. 


02  = 
0.0630 


Calculation  of  the  Average  Diffusion  Coefficients  for  Oxygen,  Do2m*  an<* 
Water  Vapor,  Dj^Om 


Using  Equation  (IV-3) 
°02m 


nn2po2-n2  *  nh2odo2-h2o 

1  -  N_ 


_  (0.839)  (0.223)  +  (0.  0630)  (0.  284)  =  0.  227  cm2/sec. 

1  -  0.0980  ■  ■  = 


and  Dn^  = 


NN2  dh2o-n2  +  no2  dh2o  - o2 


1  -  N. 


h2o 


=  (0.  839)  (0.  260)  +  (0.  0980)  (0.  284)  =  0.  262  cm2/ sec. 

1  -  0.0630 
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Calculation  of  Film  Factors  pf 
Using  Equations  (10)  and  (11) 


f  =  K +  ^  n  po 
°2  °2  °2 


where 


o2 


ft  +  <fD  pQ 

°2  °2(  1 ) 


“Tm 


”TT+  (5  o  Pq 
°2  °2 (2) 


In 


TT  +  6  p 

°2  °2(  1) 


TT-r6n  Po 
°2  U2(2) 


r02a  '  rH2Oa  2,  24  x  10~7  -  4.  48  x  10~7 


O^SL 


2.  24  x  10 


-7 


also 


’h2o 


[  1 . 02  +  (-1)  (0.200) 

]  -  [l.02  +  (-1)  (0)] 

then 

Pfo2 

in  J 

*1.02  + 

(-1)  (0.  200)' 

1.02  + 

(-  1)  (0)  J 

9 1 7  atm 


+  ^  h2oph2oJ 


lm 


[TT+  6  H20PH20(1)]  -  ^  +  ^  H20pH20(2j] 


In 


c5  H20PH20(1) 


^  H20PH20(2) 


Where 


6  HP 


rH2°a  *  °2a  _  4.  48  x  10'7  -2.  24  x  10*7 


H  ,Oa 


4.  48  x  10 


~T 


=  0.  5 


then 


pf. 


H20 


fl.  02  +  (0.5)  (0.0278)] 

_ 

[l.  02  +  0.  5  (0.  10lO 

In 

”l.02  +  (0.  5) 

(0.  0278)"" 

_1.  02  +  (0.  5) 

(0. 101) 

1, 058 
atm 
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Calculation  of  the  Rates  of  Diffusion  of  Oxygen  and  Water  Vapor  in  the 
Presence  of  Each  Other 


Using  Equation  (IV-8): 


D0 „m  Tf 


02ac  RTBg  pf 


o2  *o2 


P°2(1)'P°2(2) 


7  (0,227)  (1.02)  (0.200  -  0) 

4  “  (82.06)  (312)  (Bc  )  (0.917) 

UC>2 


thus  bGq2  =  8‘  77  cm 

D 


also  _ r  H2Oac 


pH20(l) 


pH20(2) 


RT  bg  Pf 

vJ-r  x  /-n 


HtO  ‘h,o 


-4.  48  x  10 


.7  _  (0.262)  (1.02)  Eo.  0278  -  0.  lOl] 

(82.06)  (312)  (Be,,  )  (1-058) 


thus  Bq 


h2° 


=  1.  607  cm 


lA-2b.  Estimation  of  the  Percentage  Error  Introduced  by  Calculating 
the  Diffusion  RateB  of  Oxygen  and  Product  Water  Vapor  Neglecting  the 
Effect  of  CounterdiffuBion  of  Each  Component  -  Based  on  Bulk  Gas 
Diffusion  Coefficients,  (45,  46,  47,  48,  56) 

If  water  vapor  and  nitrogen  gas  are  considered  as  a  stagnant  component 
and  the  effect  of  water  vapor  i.B  neglected,  the  rate  of  diffusion  of  oxygen 
from  the  air  stream  to  the  oxygen  electrode  is  obtained  as  follows: 

Dq2  -  N2  IT 

r°2a  ’  RT  PT  BGn  LP°2(1)  ’P°2<21- 

u2 
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and  ro  RTBr- 

2  CL 

a  2 


(0.223)  (1.02)  (0.200  -  0) 

(0.917) 


=  4.95  x  10*2 


Ratio  of  Countercurrent  Diffusion 
No  Countercurrent  Diffusion 


ro2 


RTB 


ac 


GO, 


RTB, 


5-04  *1o4=1.015 

4.  95  x  10-‘  •=— 


Percentage  Error  in  Tq  =  -^^9% 


Calculation  of  the  Rate  of  Diffusion  of  Water  Vapor  Neglecting  the 
Effect  of  the  Counterdiffusion  of  Oxygen  Gas. 

If  oxygen  and  nitrogen  gas  are  now  considered  as  a  stagnant  compon¬ 
ent  and  the  effect  of  oxygen  gas  is  neglected,  the  rate  of  diffusion  of 
water  vapor  from  the  oxygen  electrode  to  the  air  stream  may  be 
obtained  as  follows: 


-rH2Oa 


Dh2o-n2TT  [ph2o(d  ■  pH2o(2)] 

iTTTT 


Br 

1h2o  uh2o 


where 


Pf 


PN2(2)  ’  PN2  (1) 


h2o 


_ _  =  °-v919  -  °r  ?22-  =  Q.  935  atm 

'"[m)  ” 

\pN2(l)  / 


-(0.260)  (1.02)  (0.0278  -  0.  1010)  =  ?  nft  v  HT2 

then  rH  Qa  RTBGh,q  "  (0.635)  ■ 


Ratio  of  Counterdiffusion  Diffusion 
No  Countercurrent  Diffusion 


.  *.  Percentage  Error  in  =  +12I_12% 


rH2Oa  RTBGh2q 

ru  —  RTB  — 
H2°a  GH20 


_  1.855  x  lO'^ 
2.  08  x  10'Z 

=  0.893 
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Re  suit:  Neglecting  the  effect  of  water  vapor  counterdiffusion  on  the 
rate  of  diffusion  of  oxygen  gas  from  the  air  stream  to  the  oxygen 
electrode  causes  an  error  of  -1.  6%  in  the  rate  calculation.  There¬ 
fore,  the  rate  of  diffusion  of  oxygen  gas  from  the  air  stream  to  the 
oxygen  electrode  is  not  appreciably  affected  by  the  rate  of  counter¬ 
diffusion  of  water  vapor  from  the  oxygen  electrode  to  the  air  stream. 

The  error  in  the  rate  of  diffusion  of  water  vapor  from  the  oxygen 
electrode  to  the  air  stream  is  +12.  1%  when  the  effect  of  the  counter- 
diffusion  of  oxygen  gas  is  neglected. 

1 A  -  3  Conclusions: 

1A-3.  1.  Oxygen  Diffusion  from  a  Force -Convected  Air  Stream 
Flowing  at  a  Rate  Sufficient  t?  Remove  100%  of  the  Product  Water 
From  the  Oxygen  Electrode  Considering  the  Counterdiffusion  of 
Water  Vapor  -  Dased  on  Bulk  Gas  Diffusion  Coefficients. 

The  calculated  gas  film  thickness  of  8.77  cm  represents  the  maximum 
distance  that  oxygen  gas  can  diffuse  from  a  force -convected  air  stream 
to  an  oxygen  electrode  against  the  counterdiffusion  of  product  water 
vapor  at  the  100%  removal  rate  with  ar.  oxygen  partial  pressure  gradient 
of  0.  200  atm  and  the  experimental  conditions  of  cell  temperature  and 
total  pressure  and  maintain  a  current  density  of  80.  7  amp/ft  .  Since 
in  the  experimental  hardware  the  average  diffusion  path  or  gas  film 
thickness  is  only  °’  j06  6  =  0.0533  cm  or  1/164  the  maximum  path, 
the  rate  of  oxygen  diffusion  is  r.o1-  limited  by  a  stagnant  gas  film  based 
on  bulk  phase  diffusion  oefficients.  The  effect  of  the  counterdiffusion 
of  water  vapor  on  the  diffusion  cf  oxygen  is  negligible.  Theoretically^, 
on  this  basis  alone  a  limiting  '•urrent  density  of  about  13,  300  amp/ft 
could  be  maintained  under  the  same  operating  conditions.  At  one  time 
the  stoichiometric  air  requiremert  the  limiting  current  density  would 
be  10,  200  amps/ft2.  T'  erefore,  tl  e  re  ently  invented  perforated  air 
electrode  current  colleefcr  is  suitable  for  large  size  cells  and  operation 
at  high  current  densities. 

By  comparison  the  size  of  a  s<=  If- diffusion  type  air  electrode  with  an  air 
gap  of  sufficient  magnitude  (V  0.  007  in)  to  supply  oxygen  at  the  required 
rate  would  be  limited  to  a  round  shape  7.0  in.  in  diameter  t.o  operate 
a  current  density  of  80.7  amps,  ft2  at  the  same  conditions  of  tempera¬ 
ture  and  total  pressure.  The  size  cell  at  162  amps/ft  would  be 
limited  to  a  3.  5  in.  diameter.  The  air  gap  would  have  to  be  greater 
than  0.024".  Much  larger  gaps  are  required  in  practice  to  maintain  the 
average  oxygen  pressure  gradient  a.t  the  0.  200  atm.  level  required  fo 
maintain  an  average  current  density  at  these  levels  (Ref.  Task  III), 

On  the  basis  of  bulk  gas  diffusion  coefficients,  this  type  cell  design  is 
limited  to  small  sizes  and  low  current  densities. 
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1A-3.  2.  Product  Water  Vapor  Diffusion  from  the  Oxygen  Electrode 
to  an  Air  Stream  Flowing  at  a  Sufficient  Rate  to  Remove  100%  of  the 
Product  Water  Against  the  Counterdiffusion  of  Oxygen  -  Based  on 
Bulk  Diffusion  Coefficients. 

Since  the  actual  air  gap  width  is  small  —~2 -  =  0.  0533  cm  compared 

to  the  calculated  gas  film  of  1.  607  cm,  the  rate  of  diffusion  of  water 
vapor  to  the  air  stream  from  the  oxygen  electrode  under  the  experi¬ 
mental  conditions  is  not  limited  by  diffusion  through  a  stagnant  gas_^ 
layer  at  the  electrode  surface  at  a  current  density  of  80.  7  amps/ft  . 
However,  the  rate  of  removal  of  water  vapor  from  the  air  electrode  gas 
compartment  and,  therefore,  the  oxygen  electrode  to  ambient  conditions 
is  dependent  on  the  air  flowrate,  i.e.  it  must  be  sufficient  to  remove 
100%  of  the  product  water  by  maintaining  the  water  vapor  partial 
pressure  gradient.  On  the  basis  of  bulk  gas  diffusion  coefficients 
alone,  the  current  density  at  which  the  stagnant  gas  film  for^water 
vapor  diffusion  would  be  equal  to  0.  0533  cm  is  2430  amps/ft  at  the 
same  experimental  conditions  of  temperature  and  pressure  and  100% 
removal  of  product  water. 

By  contrast,  a  self-diffusion  type  cell  would  be  limited  to  a  maximum 
7.  0  in.  diameter  round  shape  with  an  air  gap^  0.  007  in.  and  ^  0.  632  in. 
in  order  to  avoid  an  oxygen  electrode  product  water  removal  problem 
at  80.  7  amps/ft2.  At  a  current  density  of  162  amps/ft  ,  the  maximum 
diameter  would  be  3.50  in,  and  the  air  gap^  0.024  in.  and  K.  0.316  in. 

1A-3.  3.  Gas  Diffusion  into  Pores  ConUolLng 

In  order  to  determine  whether  the  diffusion  of  oxygen  from  the  air 
stream  into  the  pores  of  the  oxygen  electrode  limits  the  rate  of 
diffusion  of  oxygen  an  analysis  was  made  as  follows: 

Basis: 

(1)  "GPS"  -  Type  I.  E.  M.  Fuel  Cell 

(2)  Standard  Pt  black  Catalyst  Screen  Electrode  Loading  of 
0. 223  g/in. 2 

(3)  Geometric  Catalyst  Area  =  56.0  cir/ 

(4)  Electrode  Thickness  0.0045  in. 

(5)  Catalyst  Porosity  33.  5%  by  the  Mercury  Penetrometer  Method 

(6)  Pore  Size  Distribution  by  the  Mercury  Penetrometer  Method 
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Three  Distinct  Pore  Diameter  Groups: 


Range 

% 

Avg.  Pore  Diameter 

18-90  /JL 

35.6 

35 /K 

0.  15  -  18 /*■ 

27.  4 

1.  5y<i 

0.15  -  0.06/< 

37.  0 

0.08/4 

(7)  Bulk  Density  of  Catalyst  =  3.  40  g/cm 

(8)  Total  Surface  Area  of  Catalyst  =  16  m2/g 

(9)  Diffusion  Tortuosity  Factor  =  1.50 

(10)  Assume  pores  free  of  liquid  water  except  at  the  membrane- 
catalyst  interface  at  the  bottom  of  the  pores 

1A-3.  3a  Diffusion  Tortuosity  Factors  (46,50,  51) 

Experimental  data  and  a  survey  of  the  literature  reported  by  Leonard 
Benjamin  Rothfeld  (51)  indicated  that  diffusion  tortuosity  factors  for 
unconsolidated  media  such  as  catalyst  layers  produced  by  compressing 
powders  are  from  1.4  to  3.0,  for  consolidated  synthetic  media  as  high 
as  5.  3  and  for  naturally  occurring  rocks  up  to  10.  He  concludes  that 
previously  reported  factors  as  large  as  10  for  pelleted  catalyst  were 
probably  based  on  an  incorrectly  assumed  diffusion  mechanism  and 
that  no  sound  experimental  evidence  is  available  to  support  the  conten¬ 
tion  that  the  tortuosity  of  catalyst  particles  is  this  high.  He  further 
states  that  diffusion  tortuosities  are  about  20%  larger  than  flow  tortuosi¬ 
ties. 

Air  permeability  data  reported  by  the  National  Carbon  Co.  (Ref.  p.  155C, 
Chem.  Eng.  Hd.  Bk. )  for  porous  compressed  carbon  and  graphite  slabs 
avg.  pore  sizes  from  5  to  190/C  yield  flow  tortuosities  which  level  off  at 
1.  50  for  very  small  pore  diameters  at  48%  porosity  and  at  1.  23  for  36% 
porosity.  The  corresponding  diffusion  tortuosities  are  estimated  as 
1.80  and  1.48,  respectively. 

The  experimental  facts  are  in  good  agreement  with  a  theoretical  tortuosity 
factor  of  1.41  based  <3n  a  "random"  direction  assumption  by  Ahlborn 
Wheeler  (58). 

Finally,  on  this  basis  a  diffusion  tortuosity  of  1.  50  is  assumed  for  the 
experimental  catalyst.  That  is  the  average  length  of  the  gas  diffusion 
path  be  assumed  to  be  1.  50  times  the  catalyst  electrode  thickness. 
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1A-3.  3b  Gas  Phase  Diffusion  Coefficients  in  Catalyst  Pores  (50,  51) 


( 1)  Laminar  Flow  Regime 

The  magnitude  of  the  gas  bulk  diffusion  coefficient  in  pores  of  10,000  A° 
(ljA  )  radius  or  larger  in  contact  with  one  atmosphere  of  gas  pressure 
is  independent  of  the  pore  radius  and  inversely  proportional  to  the 
total  gas  pressure.  Because  under  these  conditions  the  mean  free  path 
of  a  gas  molecule  is  about  1/10  (oxygen  =  9.940  x  10  ^  cm  =  994  A°  = 

0.  0994^  at  20°C)  the  pore  radius  and  molecules  will  collide  much 
more  often  with  each  other  than  with  the  pore  wall.  Laminar  or  stream¬ 
line  flow  will  predominate.  The  kinetic  theory  gives  the  following 
approximate  formula  for  the  diffusion  coefficient  for  a  mixture  of  gases 
of  similar  mass  and  molecuiar  diameter 

°B  ~T  ■ 

where  Dg  =  bulk  diffusion  coefficient  for  gases,  cm2/sec.  (IV- 13) 

( \jr  -  arithmetic  average  molecular  velocity,  cm/sec. 

\  =  mean  free  path  (cm) 


An  explicit  equation  for  Dp  in  terms  of  molecular  quantities 
and  the  concentration  or  pressure  is  obtained  by  substituting 


0.  707 

tt  a2  CT 


(IV-  14) 


where  Tf  =  total  ga3  pressure,  (atm) 

C  -  molecular  diameter  (cm) 

*2 

CT  =  total  concentration  (molecules/cm-5) 


eqn.  IV- 14  into  eqn.  IV- 13 
0.  707 


whence  = 


or 


3ir 


2  r 

c r  cr 


(IV- 1 5) 


The  range  of  values  for  Dg  for  most  gas  mixtures  at  atmospheric 
pressure  and  a  20°C  temperature  are  0.  1  to  1.0  cm2/sec. 
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(2)  Knudsen  Flow  Regime 


When  a  catalyst  contains  pores  of  1000  A  (0.  1X()  radius  or  less 
Knudsen  or  molecular  flow  rules  the  rate  of  gas  transport  at 
moderate  pressures.  Under  these  conditions,  the  mean  free 
path  between  intermolecular  collisions  will  be  greater  than  the 
pore  diameter.  As  a  result  a  molecule  within  the  pore  structure 
will,  in  general,  strike  a  pore  wall  before  it  strikes  another 
molecule.  Intermolecular  collisions  may  thus  be  neglected, 
and  a  molecule  travels  within  a  pore  by  a  series  of  "random 
flights"  interrupted  by  collisions  with  a  pore  wall.  The  "resist¬ 
ance"  to  diffusion  is  due  to  the  fact  that  after  collision  with  the 
pore  wall,  the  molecule  is  just  as  apt  to  reverse  its  direction  as 
it  is  to  proceed  in  the  direction  of  the  previous  flight.  Knudsen 
showed  that  under  these  conditions  the  flow  of  molecules  through 
a  capillary  of  length  A  X  cm.  is  governed  by  the  equation: 


where  dn 
dt 


gas  flowrate 


molecules 

sec. 


(IV-  1 6) 


m  =  mass  of  a  single  molecule  (g) 

w  =  Boltzmann  Cons*.  ( 1 .  38  x  10' 1  6  ergs /°K/molecule) 

B 

(i.e.  ,  the  gas  divided  by  Avogadro's  Number) 


T  =  absolute  temperature  (°K) 

r  =  radius  of  the  capillary  (cm) 

^  p  =  pressure  difference  between  ends  of  the  capillary 

(dynes/ cm^) 

Knudsen* s  equation  shows  that  flow  or  diffusion  of  a  particular  molecular 
species  in  small  pores  is  independent  of  the  total  gas  pressure,  is  inde¬ 
pendent  of  the  pressure  or  absence  of  other  types  of  molecules,  and 
depends  only  on  the  partial  pressure  gradient  of  that  species.  Because 
each  molecule  travels  completely  independently  of  its  neighbors,  there 
is  no  such  thing  in  small  pores  at  moderate  gas  pressures  as  a  "mass  ' 
or  "hydrodynamic"  flow  in  which  molecules  are  carried  along  in  propor- 
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tion  to  their  concentration,  such  as  occurs  in  ordinary  streamline 
flow.  A  Knudsen  diffusion  coefficient  can  be  defined  by  using  the 
perfect  gas  law 


p  =  c  k  T 
y  B 


and  defining  this  coefficient  for  molecular  species  in  a  pore  by  the 
equation  of  def: 


dn 

dt 


=  A 


(IV-  1 7) 


where  A  is  the  cross-sectional  area  of  the  pore  then 


(IV-  18) 


For  a  molecule  of  average  molecular  weight  A?  is  about  10^  cm/ sec 
an<^  Dfc  has  the  numerical  range  of  values  of  about  0.01  cm^/sec  for 
10a  (0.  00  1  JU^)  to  1,0  cm^/sec  for  10,  000  (1.44,)  pores. 


(3)  Over-All  Equation  for  Gas  Diffusion  Coefficients  in  Pores 


Wheeler  (58)  proposes  two  semi-empirical  over-all  equations  for 
gas  diffusion  coefficients  in  pores  as  follows: 


and 


D  = 


D  = 


4-  AT 


A(-  ) 


—  AT  A( 


'2r/A 

1  +  2r/X 


(IV  - 19) 

(IV -20) 


both  have  the  property  of  showing  the  transition  between  Knudsen  and 
ordinary  diffusion  and  of  reducing  to  the  two  special  cases;  eqn.  IV- 18 
when  r  is  much  less  than  A  and  eqn.  IV- 13  when  r  is  much  greater  than 
^  .  Equation  IV- 19  shows  a  more  gradual  transition  between  the 
two  type  of  diffusion  than  does  IV-20,  yet  he  prefers  equation  IV- 19 
since  the  physical  problem  is  surely  connected  with  the  relative  prob¬ 
ability  of  a  molecule  striking  the  pore  wall  vs.  striking  a  second 
molecule  and  the  bracketed  expression  in  IV- 19  is  exactly  the  probability 
that  a  molecule  will  have  a  collision  with  a  second  molecule  before 
traveling  a  distance  equal  to  the  pore  diameter. 
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Over-all  diffusion  coefficients  calculated  using  equation  IV- 19  in 
the  form  D  =  D]j  (1-e  were  within  10%  and  0.  5%  of 

the  values  determined  experimentally  by  Rothfeld  (51)  on  N2~He 
and  N2  -  A  counterdiffusing  mixtures,  respectively.  It  is  expected 
that  for  mixtures  of  oxygen  and  nitrogen  and  water  vapor  with 
molecular  weights  and  sizes  much  less  different  than  the  pairs 
Rothfeld  reported  on  that  the  use  of  equation  IV- 19  will  permit 
calculation  of  over-all  diffusion  coefficients  within  less  than  10% 
of  the  actual  values. 


(4)  Application  of  Equations 

System:  Same  as  (1)  Bulk  Diffusion  -  Gas  Film  Controlling 


(a)  Calculation  of  Over -All  Diffusion  Coefficients  for  Each  of  the 
Three  Average  Pore  Size  Ranges  (45,  46,  47,  50,  55,  56) 

-2: 


now  Dn  =  DQ  _ 

°2m  (0.  A)  U2m 


and 


312.  4°K 
Tf  =  1.  02  atm 

<V~  Oz  =44 
312.  4°K 
Tf  =  1.02  atm 


r\rcx^W~ .  . 


^  .  -2r??/3D02l^ 


L.  dUU  - 

•000  vte)  fsH)' 


312.  4°K 
1.  02  atm 


=  45,  000  cm/ sec. 


AT 


fl/"  =  47,  100  (1.022)  =  48,200  cm/sec. 

ATH  o  =  58,  700  (1.022)  =  60,000  cm/sec. 


avg. 


=  AT^ 

°2  °2 


+  WV  +  N' 


n2-n2 


or 

h2°  h2o 


=  (0.0980)  (45,000)  +  (0. 839)  (48, 200)  +  0.0630  (60,000)  =  48, 69_0  cm/ sec. 


(b)  then  for  35.M.  avg.  pore  diameter 


®2  m  (O.A.)  =  0.227 
312.  4°K 
TT=  1.02  atm 


-2(1.75  x  10~3)  (4,869  x  104) 
3  (0.  227) 


1  -X 


D02m 


0.  227  cmVsec. 
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(c)  also  for  1.  5>Uavg.  pore  diameter 


(0.  A. ) 
312.  4°K 

7f  =  1.02  atm  =  D( 


=  0.227 

_  =  0.  227  cm2/sec. 

02m  - 


- 2(7,  5  x  10~5)  (4.869  x  104) 
3  (0.227) 

1  -  JL 


(d)  Finally  for  a  0.08M,avg.  pore  diameter 


Do.  in  a  \  =0.  227 

O2m{0.  A.  ) 


1  -  JL 


-2(4  x  IQ'6)  (4,869  x  104) 
3  (0.227) 


=  0.  0992  cm2/sec.  ^  D°2m  =  — 227  cm  'Bec* 


and  D  =  0.262  (0.391)  =  0.  1024  cm2/D..  ~  =0.262  cm2 

H2Om(0.A.)  ==sTc.N  H2°m  =  «c. 


Note: 

For  the  0.  08j^avg.  pore  diameter,  the  over-all  pore  diffusion  coefficient 
is  56.  5%  less  than  the  bulk  diffusion  value  for  oxygen  and  60.9  less  for 
water  vapor. 

1 A  -  4  Maximum  Current  Density  Limitations  Due  to  Gas  Pore  Diffusion 
Mechanism  (45,  50,  51) 

The  use  of  the  over-all  gas  pore  diffusion  coefficients  for  oxygen  and 
water  vapor  to  estimate  the  maximum  current  density  possible  for 
"GPS"  I.  E.  M.  fuel  cell  with  the  standard  catalyst  electrode  operating 
on  force-convected  air  under  the  same  conditions  stated  for  the  bulk 
diffusion  case  is  as  follows: 

Since  the  screen  of  the  catalyst-electrode  occupies  about  16%  of  the 
geometric  surface  area  of  the  electrode-membrane  interface  and  pene¬ 
trates  to  about  1/3  the  thickness  of  the  catalyst,  the  average  pore  length 
will  be  somewhat  longer  0.  01728  cm  vs.  0.  01712  cm  to  the  reaction 
zone  which  is  assumed  to  occur  in  the  limit  at  the  bottom  of  the  pores 
at  the  catalyst-electrolyte  interface.  The  maximum  rate  of  diffusion 
and  of  current  density  will  then  be  determined  by  a  weighted  average 
of  the  rates  of  diffusion  through  the  two  sections  of  different  length  and 
area. 
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Now  the  total  gas  flowrate is  the  same  for  both  sections 

hence  \  1.0/  \0.84  /  „  . 

O2a(0.  A.  )  avg.  -  -  —  r -  rO2a(0.A.  )1 

\~T  ) 


=  1.  127  r 


O 2*  (0 • A. )  j 


g-moles/cm2-sec. 


al8°  rH2Oa  (O.A.)  avg.  =  1.127  A  ^ 

°02m  (0.  A. )  avgTT 

•  _  _ 

•  •  r02a  (O.A. )  avg. 


2(1) 


’2(2)] 


RT  ",  Bc, 


'2  (O.A.  )  avg. 


>re  size 


treating  the  Dq  ^(O.A.)  for  the  three  differcnt  average  poi 
as  conductance  £  in  parallel. 

DO2ro(0.A.Uv8.  =  (por.^action)  Dq^o.  A. ,  + (por.  fraction)  X 

DO2m(0.A.)  4  DO2m(0.A.) 

=  0.  356  (0.  227)  +  0.  274  (0.  227)  +  0.  370  (0.  0992) 

=  0.  1797  cm2/  sec. 


and  D, 


H^Om  (0.  A. )  avg. 


(0.  63)  (0.  262)  +  (0.  37)  (0.  1024)  =  0.  1688  cm2/sec. 


_  (0.  1797)  (1.  02)  (0.  200  -  0)  _ =2.  0u 

(82.  06)  (312)  (0.  917)  (2.  24  x  IQ-7)  (1.  127)  Cm 


then  B 

°2(0.  A, )  avg, 
at  80.  7  amps /ft 2 


0.  335 
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(0.  1688)  (1.02)  (0.  101  -  0.0278) _ 

also  (o.A.  )__._  (82.06)  (312)  (1.052)  (4.48  x  10  )  (1.  127) 

£  g'  0.335 


=  0.  310  cm 


Since  the  length  of  the  diffusional  path  is  only  0.01728  cm  in  the  experi¬ 
mental  hardware  and  the  calculated  stagnant  gas  film  thickness  is  2.06 
cm  for  oxygen,  the  rate  of  diffusion  of  oxygen  to  the  bottom  of  the 
catalyst  pores  is  not  limited  by  the  rate  of  diffusion  through  this  film 
and  neither  is  the  current  density  of  80.7  amps/ft2  at  the  experimental 
conditions  of  100%  product  water  removal,  temperature,  and  pressure. 

In  fact,  on  the  basis  of  oxygen  pore  diffusion  coefficients  alone,  the 
theoretical  limiting  current  density  would  be  9,600  amps/ft  .  The  rate 
of  product  water  vapor  diffusion  from  the  bottom  of  the  catalyst  pores  to 
the  air  stream  is  not  limited  by  diffusion  through  a  stagnant  gas  film 
in  the  pores  at  a  current  density  of  80.7  ampB/ft2.  Theoretically,  the 
current  density  at  which  water  vapor  pore  diffusion  does  become  limiting 
is  1,  450  amps/ft2  with  100%  product  water  removal  and  the  same  condi¬ 
tions  of  temperature  and  pressure.  However,  the  rate  of  removal  of 
water  vapor  from  the  oxygen  electrode  and  from  the  air  electrode  gas 
compartment  does  depend  on  the  air  flowrate,  i.  e.  it  must  be  sufficient 
to  remove  100%  of  the  product  water  under  the  experimental  conditions 
of  temperature,  water  vapor  partial  pressure  gradient,  and  total 
pressure. 


1-B.  Oxygen  Diffusion  Rate  From  An  Air  Stream  to  the  Oxygen 
Electrode  Limited  by  Gas  Diffusion  Across  a  Liquid  Film 
at  This  Surface.  (45,  46,  47) 


The  basic  equation  for  gas  diffusion  through  a  liquid  is  given  by  Fick  s 
law 


rAa 


AB 


(IV-21) 


where  r^^  =  rate  of  diffusion  of  A  (g-moles/ sec-cm  ) 


CA  =  concentration  of  A  (g-moles/cm3) 
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D  .  =  diffusion  coefficient  of  gas  A  in  the  liquid  phase  (cm^/sec. ) 

AB 


since  C, 


Xa 


"A  flT m 

where  XA  =  Mole  fraction  of  A  in  the  liquid  phase  (g-moles/g-moles) 

=  average  molal  volume  of  liquid  phase  (cm^/g-mole) 
d^  =  effective  liquid  film  thickness  (cm) 

Subscript  i  =  refers  to  properties  at  the  interface 

Equation  21  becomes: 

n  H  v  . 

(IV-22) 


Aa 


AB 

flAn 


__A 

dL 


Equation  IV-22  may  be  integrated  over  a  film  thickness  BL  to  give: 


AB 


Aa 


^mBL 


(Xa  ■  -*Ai) 


(IV- 2  3) 


This  equation  may  be  applied  to  the  diffusion  of  oxygen  through  a  liquid 
film  using  the  same  air  data  as  in  the  bulk  gas  diffusion  case 


D~  uo  =  1.  98  x  10'5  SELL  +  (28.  3°C)  (0.03  *  )  (1.98  x  10’5) 

-  « 2'-'  sec.  C 

(at  3 19.  3°K)  (at  291°K) 

=  3.  66  x  10'5  •  (Ref.  Int.  Crit.  Tables  Vol.  V) 

sec. 
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X02  at  the  aye  rage  conditions  of  the  surface  of  the  liquid  film  is 
estimated  using  Henry's  law  which  is  followed  by  oxygen  at  low 
and  moderate  pressures. 


h2o  Xo2 


(IV-24) 


where  K  depends  only  on  the  solute -solvent  pair  and  the  tempera¬ 
ture.  It  is  independent  of  the  pressure  up  to  5  atms. 

0.  200 

5. 70  x  10 

(at  319.  3°K) 
also  X  o2.  =0 


uiereiure 


Xo, 


4  =  3.  51  x  10“^  (g-moles/g-moles) 


The  oxygen  gas  mole  fraction  at  the  catalyst-electrolyte  interface  can 
be  safely  assumed  to  be  negligible  compared  to  the  bulk  mole  fraction; 
therefore,  to  approximate  the  oxygen  concentration  gradient  X02^ 
is  set  equal  to  zero. 

cm  ^ 

and  Ar  »  AT  =  18.  20 - —  the  molal  volume  of  water  since 

fv*  m  mHzo  g-mole 

the  molal  volumes  of  oxygen  and  nitrogen  in  the  liquid  phase  can  be 
neglected. 

then  r~  ln  .  ,  (based  on  a  33.  5%  catalyst  porosity  and 

02am  (0.  A. )  avg. 

correction  for  screen) 

(1.  127)  (2.  24  x  10'7)  =  (3.  66  x  10~5)  (3.  51  x  10~b  -0) 

(0.  335)  18.  20  Bl 

.  Bt  =  8. 69  x  10~k  cm 

«  •  L  - 

at  a  current  density  of  80,  7  amps/ft^. 
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The  significance  of  thi6  liquid  film  thickness  may  be  estimated  as  follows: 


Avg.  Pore  Diameter 

%  of  Total 
Pores 

Total  Pore 
Geom. 

Surface  Area 

%  of 
Total 

Ratio  213^ 
Pore  Diam 

Area  of 
Reaction 
Zone,  cm2 

Length  of 
Reaction 
Zone, cm 

35X\=  3.  5x  10'3ctti 

35.  6 

1. 108xl02crti2 

0.  2 

0.00487 

7.  12 

1. 108x10' 

1.  5*^=  1.  5x  10'4cm 

27.  4 

1.  993xl03cm2 

3.  7 

0. 1158 

5.  47 

4. 75x10' 3 

0.08Ul^=  8.  Ox  10'6cm 

37.0 

5.  05x  104cm2 

96.  1 

2.00 

3.  21 

1. 10x10' 6 

The 

liquid  film 

thickness  represents 

0.487%,  11. 

58%,  and  200%  of 

the  diameters  of  the  35>t(,  1.  5A^,  and  0.08JL^  average  pore  sizes, 
respectively.  Also,  the  geometric  surface  area  of  this  film  must 
be  equal  to  15.80  cm2  from  the  diffusion  equation.  It  will  be  instruc¬ 
tive  then  to  compare  this  value  first  with  the  total  pore  geometric 
surface  area  and  the  reaction  zone  area  for  each  pore  size  with  the 
total  area  for  each  size  and  second  the  length  of  the  reaction  zone 
with  the  pore  length  for  each  pore  size. 


the  number  of  35_A^pore  s  /rL  = 


(0.  356)  (0. 335)  (56. 0)  (0. 84)  (0.  01728) 
~jr  (3.  50  x  10-3)2  (0.01728) 

4 


=  5. 84  x  10  pores 


and  the  number  of  1.  5/*  pores 


=  (0.274)  (0.335)  (56.0)  (0.  84)  =  *  £(  1.  5  x  10 


V 


Or 

*.  4S' k 

*  *nI.5><-  =  2.  45  x  108  pores 


finally/\o^ 


(0.  370)  (0.335)  (56.0)  (0.84)  =  1#  l62  x  ioUpores 
(0.  785)  (8.0  x  10-b)  i  — 


The  total  geometric  pore  surface  area  is 

=  A35^  +  A1.5>4  +  A0.08X{ 

=  (5.  84xl05)  (If)  (3.5xl0'3)  (1.728X10-2)  +  (2.  44xl08)Iff.)  (1.5X10'4) 
(1.  728  x  102) 

+  (1.  162x1 01 !)  Ctf)  (8.  0x10'°)  (1.  728xl0'4)  =  S.  25x10  cm 
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A  comparison  between  this  number  and  the  liquid  film  area  shows  that 
the  entire  diffusion  of  oxygen  and  reaction  occurs  at 

*5,  80  _  =  0.0301%  of  the  total  geometric  pore  surface  area. 

5.  25  x  104  ■■ 


The  liquid  film  surface  area  for  each  average  pore  size  is  estimated 
from  both  the  individual  cross-sectional  pore  area  fractions  and  the 
pore  diffusion  coefficients.  Whence 


Liquid  Film  Area  for  35  jA  pore  size 

1.5M  pore  size 
0.  08 JA  pore  size 


=  7.  12  cm2 
=  5.  47  cm2 
s  3.  21  cm2 


Thus  the  percentage  total  pore  geometric  surface  area  occupied  by 
the  liquid  film  reaction  zone  is  6.  4%,  0.  27%,  and  0.  006%  for  the 
35A(,  1.57^,  and  0.  08L(  average  pore  sizes,  respectively. 

The  length  of  the  reaction  zone  is  estimated  as  follows: 

35  pore  diameter  m.7Tbli=  7.  12 

_ _  7.  12 _  =  l.  108  x  10'3  cm 

(5.84  x  10=)  (3.  14)  (3.  5  x  10-3)  ========= 


5.  47 

(2.45  x  108)  (3.  14)  (1.  5  x  10-4) 


4.  75  x  10~  ^cm 


3.21 


0.  (1.  162  x  101 1)  (3. 14)  (8. 0  x  10'6) 


=  1.  10  x  10'^cm 


These  values  represent  6.42%,  0.275%,  and  0.00638%  of  the 
1.5^(,  and  0.  08t^  pore  lengths,  respectively. 


Conclusions: 

_  f. 

The  reaction  zone  varies  from  1.  108  x  10  cm  to  1.  10  x  10  cm 
in  length  and  covers  from  6.  4%  to  0.  006%  of  the  individual  pore  areas 
or  about  0.03%  of  the  total  geometric  pore  surface  area. 
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Theoretically,  the  minimum  thickness  for  thisgfilm  is  the  thickness 
of  a  monomolecular  layer  of  water  =  2.  7  x  10  cm  but  data  in  the 
literature  indicate  that  water  vapor  adsorption  on  platinum  is 
multimolecular  and  consists  of  from  51  to  380  layers  or  from 
1.  375  x  10“°  cm  to  1.  025  x  10'5  cm  at  the  saturation  pressure  for 
water  vapor.  This  would  appear  to  be  the  range  of  thicknesses  for 
the  liquid  film  in  the  catalyst  pores  at  the  reaction  zone.  ^Hence  the 
limiting  current  density  could  range  from  about  68  amp/ft  for  a  non 
wet-proofed  or  non-uniformly  wet-proofed  catalyst  to  as  high  as 
510  amps/ft2  in  the  ideal  wetproofed  catalyst  case  both  operating 
at  100%  removal  of  product  water.  Therefore,  the  maximum  rate 
of  oxygen  diffusion  from  an  air  stream  to  the  oxygen  electrode 
catalyst-electrolyte  interface  and  the  maximum  current  density 
possible  at  100%  product  water  removal  on  an  ion  exchange  membrane 
fuel  cell  are  limited  by  the  diffusion  of  the  oxygen  through  a  thin 
liquid  film  on  the  walls  of  the  pores  at  the  catalyst-electrolyte  inter¬ 
face. 
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APPENDIX  V  1 


Heat  and  Mass  Transfer  in  an 


Air-Breathing  Fuel  Cell  System 
with  Air  Recirculation 


V  - 1  Introduction 

In  the  analysis  which  follows,  reference  will  be  made  to  various 
flow  stations  which  are  shown  in  Figure  V-l  which  is  a  schematic 
drawing  of  the  system.  The  reader  is  referred  to  Section  4.  4.  2 
for  a  description  of  the  basic  system. 


V-2  Air  Composition  and  Properties 

The  following  composition  of  dry  atmospheric  air  was  used 


Table  V-l  -  Standard  Composition  of  Dry  Air 


Constituent  Notation 


Molecular  Weight,  m  Mole  Fraction,  X 


Nitrogen 

N2 

28. 016 

Oxygen 

°2 

32. 000 

Argon 

Ar 

39. 944 

Carbon 

co2 

44. 010 

Dioxide 


Dry  Air  28.966 


0.  7809 
0.2095 
0. 0093 
0.  0003 


1.  000 


Since  nitrogen.,  argon,  and  carbon  dioxide  do  not  enter  into  the 
reactions,  they  may  be  considered  diluent  gases  and  lumped  to 
gether  for  mathematical  convenience. 


V-3  Mass  Flow  Balance 

Assuming  that  the  flow  is  one  dimensional  at  any  flow  station  in  the 
system,  the  mass  rate  of  flow  of  any  constituent  can  be  represented 
by  the  following  expression 

W.  =  mj  Pj  (lb.  m/sec.)  I1) 

J  1545.  32  Tk 
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where:  W-  =  mass  flow  rate  of  any  constituent 

J  (j  =  Oz,H2  or  D) 

m.  =  molecular  weight  of  the  constituent 

p.  =  partial  pressure  of  the  constituent,  psia 

A^  =  flow  c ross-sectional  area,  sq.  in. 

(k  =  station  1,2.3,  etc .  ) 


=  velocity,  fps 

T  =  absolute  temperature,  deg.  R 
k 


V-4  Product  Water  Mass  Flow  Balance 


The  increase  in  mass  rate  of  flow  of  water  flowing  through  the 
system  rate  of  generation  of  product  water  if  the  bleed  flow  rate 
through  the  system  is  controlled  such  that  product  water  is 
carried  off  in  the  bleed  system  at  the  same  rate  at  which  it  is 
generated.  The  product  water  generation  rate  is  a  function  of 
the  oxygen  consumption  rate  and  since  2H2  +  -+■  2H20 


and  substituting  Equation  ( V -  1 )  into  Equation  (V-Z) 


A6V6  (  \ 

V  H2°'  6 


!A6V6  .  ,  A’V7 

=  2Nrip°2v— 


(V-3) 


and  solving  Equation  V~3  for  6 ,-^J  gives. 

WV7'Tfc 

K(V6(TT  _  ( PHZo)  7  +  Z\P°Z h 

'  -  V)4  t2fp°2)6 


(V-4N, 
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V-5  Diluent  Mass  Flow  Balance 

The  mass  flow  rate  of  the  diluent  gases  in  and  out  of  the  system 
must  equal  in  order  for  equilibrium  conditions  to  prevail  or: 

<WD>7  =  <WD>6 


(PD>7  =  P7-(PH2<^  7{P02)7 

(V-6) 

(pD)6  =  P6-(PH2o)6  ~(P0^}b 

where  and  p^  represent  the  summation  of  all  the  static  partial 
pressures  at  the  two  stations;  then  from  Equations  V-l  and  V-5: 

C^TJ££\.  ip°*h 

^-(PH2^6fP02)6  ' 

If  Equations  V-4  and  V-7  are  equated  and  expanded,  several  terms 
cancel  out  giving: 

[2P’-(PH2^7J  = 

P6fPH2()7  =  P02  7  [2p6  '  fVj  6 1 

and  solving  Equation  V-8  for  pQ^  ^  gives.  (V-8) 


But  for  the  atmospheric  air  entering  the  system: 

(po2y°-20^tH2°)\i 


(V-10) 
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and  Equation  (V-9)  becomes: 


0.  2095 

Ip 7  -  °-  5| 

t 

(pH24][»6  <»H2o)6 

-  0.  5 

r 

?6(Ph2^7  '  P?l 

K°)( 

L 

[P6  ‘  °,5fPH2o)( 

s] 

Equation  (V - 1 1)  is  a  quasi-rigorous  expression  for  the  relationship 
showing  that  the  oxygen  partial  pressure  at  the  discharge  of  a  fuel 
cell  stack  is  essentially  a  function  of  the  fuel  cell  operating 
temperature  and  the  atmospheric  conditions  when  the  amount  of 
water  leaving  the  system  is  equal  to  the  amount  entering  plus  the 
amount  generated. 


A  simplified  version  of  Equation  ( V - 1 1)  which  is  more  suitable  for 
computation  in  a  control  system  may  be  developed  by  multiplying 
and  dividing  Equation  ( V - 1 1 )  byj^p^  +  0.  ^(PH^o)^]  ’  P6  =  P7’ 

and  neglecting  all  terms  where  the  producf  of  two  values  of  vapor 
pressure  occurs  gives: 

(pq?)7  ■  0-20,5[p<> '  °-  5K°^[P6  -  K°)6][p<.  *  °-  5(M^ 


-  0.  5 

[p&  (PP 

[2o)7  1 

[Ph2<^ 

fP6  '  0 

'  1 

^p6  +  °*  5  fPH2d  6 

which  first  reduces  to: 


(  P°2)  *  2095  P6  '  °'  60475''Ph2^  7  +  39525fPH2°.  6 


(V-  12) 


A  test  of  the  accuracy  of  Equation  (V-12)  may  be  made  by  comparing  the 
results  of  an  extreme  example  with  Equation  ( V - 1 1 ) .  Let  p^  =  p7 
=  14.  696,  t^  =  120  F,  and  1 7  =  130  F  where  for  saturated  conditions 
fu  \  =  1.  6924  psia  and  p„  ^  =  2.  2225  psia. 

(PH2C)6  H  VH29? 


Relative  Vapor  Press.  Oxygen  Partial  Press. 


Humidity 

at  Inlet 

at  Outlet 

,p^'  ,  psia  Error-Per  cent 

"V 

(M  6’  Pi“ 

Eq.  (11) 

eI\u) 

0 

0 

1. 7348 

1.  7348  0.0 

1.0 

1.  6924 

2. 3909 

2.4037  0.53 
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Equation  (V-12)  can  be  made  still  more  accurate  by  adjusting  the 
last  constant  from  0.  3952  to  0.  3877  or: 


V  -6 


Po2')^0.2°95p6  -  0.6048^^*  0.3877^^ 

Oxygen  Mass  Flow  Rate 


(V  - 1 3) 


It  is  instructive  to  compare  the  oxygen  mass  flow  rate  through  the 
system  with  the  rate  of  consumption.  Let  the  through-flow  oxygen 
equivalence  ratio  o-  be  the  ratio  of  the  mass  flow  rates  of  oxygen 
entering  the  system  and  oxygen  consumed  within  the  system.  Then 


(wo2)6  -(-<*), 


and  using  Equation  V-l). 


A6  v6 


A7V, 


’,(.P°?)7  T7 


(V- 1 5) 


Rearranging  Equation  (V-15)  and  equating  to  Equation  V-4)  gives. 


(PH2o)  7  *  j[p°z)7  /  -,(^7 

iS4t  *  Ni  ^ « •  ‘  ''  (P°8l« 


(V  - 1 6) 


\y6 

When  solved  for/p  \  ,  it  becomes: 

V  7 

(PoN 


*  t 


fpo2.6 


(V  - 1 7) 
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Equations  which  are  similar  to  Equations  (V-16)  and  (V-17)  can  be 
developed  from  Equations  (V-7)  and  (V- 15);  thus- 


(V-18) 


(V-19) 


Setting  Equation  (V-17)  equal  to  Equation  (V-19)  and  eliminating 
common  factors  gives: 

_  .  LP7  '  CPiWt  j 

<r.  f  6  4  f  'Oj)6  -(PH2^fc] 

and  solving  for  gives 


P6(PH2^7'P7  (^H2c)6 

and  since  =  0.2095^-^. 


(V-21) 


Equation  (V-21)  becomes 


p6(pH20|7  ’ 

\  ,* 

Equation  (V-22)  shows  that  the  ratio  of  through-flow  rate  through 
the  system  to  the  stoichiometric  oxygen  flow  rate  requirement  is 
entirely  a  function  of  the  static  pressures  and  water  vapor  pressures 
in  and  out  of  the  fuel  cell  system 


’  [W  *  ( PH  Q 

7  b  \  h2,  6 


pvv7 


(V -22) 
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V-7  Heat  and  Mass  Balance  of  Flow  Through  the  Fuel  Cell  Stack 

The  heat  and  mass  balance  of  the  flow  through  the  fuel  cell  stack  is  an 
iterative  type  of  solution.  Certain  assumptions  can,  however,  be 
made  that  give  a  reasonably  accurate  initial  solution.  A  second 
solution  can  correct  for  any  discrepancies,  the  assumptions 
used  in  the  following  analysis  are: 

1.  One -dimensional,  steady-state  flow  through  the  air  passages. 

2.  The  initial  temperature  of  the  water  evaporated  into  the  air 
stream  is  equal  to  the  inlet  gas  temperature. 

3.  All  of  the  stack  heat  losses  are  rejected  to  the  gas  stream  by 
means  of  sensible  heat  losses  and  evaporation. 

4.  The  gas  is  saturated  with  water  vapor  at  the  entrance  and  exit 
of  the  stack. 

Regarding  the  last  assumption,  when  the  recirculating  air  flow 
rate  is  .  100  times  that  necessary  to  provide  the  stoichiometric 

oxygen  requirement,  approximately  72  per  cent  of  the  heat  rejection 
is  by  means  of  evaporation.  A  small  error  in  the  final  gas 
temperature  would  therefore  have  a  small  effect  on  determining 
the  amount  of  water  needed  to  provide  the  necessary  heat  rejection. 

V-7-1  Diluent  Mass  Flow  Balance 


Equation  (V-7)  can  also  be  applied  to  the  conditions  entering  and 
leaving  the  fuel  cell  stack,  or  for  constant  area  flow. 


V-7-2  Oxygen  Mass  Flow  Balance 

The  amount  of  recirculation  is  related  to  the  stoichiometric  oxygen 
requirement  by  the  ratio  as  follows: 


-  184  - 


4.  0  FACTUAL,  DATA  (Cont'd) 


V.7.2  Oxygen  Mass  Flow  Balance 

The  amount  of  recirculation  is  related  to  the  stoichiometric  oxygen 
requirement  by  the  ratio  as  follows: 


r°j2 

(W°A  •  (*°A 

where:  /W  \  =  oxygen  flow  rate  leaving  the 

'  ^2*2  fuel  cell  stack 

(\y  \  -  fw  \  =  oxygen  consumption  in 

°2'1  '  the  fuel  cell  stack 


(V-24) 


or  rewriting  Equation  (V-24): 


(V-25) 


and  using  Equation  (V- l)  for  the  continuity  expression  for  constant 


area  flow: 


(V-26) 


Rewriting  Equation  (V-26)  gives 


(V-27) 
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The  partial  pressure  of  the  oxygen  at  the  stack  entrance  can  be  found  by 
setting  Equation  (V-23)  equal  to  Equation  (V-27),  or: 


P2  | 

(Ph2°)2  ‘  1 

IP°z1 

Iz  . 

i! 

w 

^  (V-28) 

Pl'l 

("Hz0)!  ‘  1 

(P°z) 

i,  '  v  ■« ' 

'i 

:p°zi 

1, 

which  reduces  to  the  following  expression  when  solved  for 


(ffc  +  [P1  ’  (PH2o)1 

(P°z)z 

•c  [pj  -  ev)2j  * 

(V-29) 


Substituting  Equation  (V-29)  for 


in  Equation_(V-27)  gives: 


<T 

C 

(V)J+( 

P°z)z 

(T 

C 

P1  -  (pH 

Z°)l 

(V-30) 


V.  fft  Temperature  Change  with  Saturated  Conditions 

If  the  gas  is  saturated,  the  partial  pressure  of  the  water  vapor  is 
only  a  function  of  gas  temperature.  A  change  in  temperature  is 
approximately  proportional  to  the  change  in  the  partial  pressure  of 
the  water  vapor  for  small  temperature  changes.  The  approximation 
can  be  expressed  as  follows: 
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T  — 


(PhA  '  (V)i 


8  Pi 


(V-31) 


H2° 


8T 


sat 


The  evaluation  of  (8  p„  _/8T)  from  Keenan  and  Keyes  is  as  follows 

H^O  sat 

for  the  temperature  interval  from  118  F  to  120  F: 


sat 


1.  6924  -  1.  6006  _  0.  0918 
120.  0  -  118.  0  '  2.  0 


0.0459  VE£W  (V-  32) 
deg.  F 


and  from  Equation  (V-31): 


(V-  33) 


Substituting  for  T/T-  in  Equation  (V-30)  using  Equation  (V-33)  gives: 
1  ^ 


f 

(T 

C 

f2' 

(v)2. 

+ 

,  (Ph2°)2 ' 

(PH2°)i 

%E 

>1  -  (pH2 

°\] 

> 

'ra 

lt  x2 

(V-  34) 


V,  9  Fuel  Cell  Stack  Heat  Balance 

At  120  deg.  F  the  heat  of  formation  of  condensed  water  is  =  68,  204 
cal/gm-mole  which  corresponds  to  an  ideal  voltage  of  E.  =  1.478  volts 
using  the  Faraday  equation.  The  fuel  cell  thermodynamic  efficiency  T)p£ 
can  be  expressed  as  the  ratio  of  the  output  voltage  E  and  Ei  or 
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E  _  E 

^FC  E.  1.478 


(V-35) 


The  amount  of  heat  rejected  by  the  fuel  cell  stack  qR  can  be  related 
to  the  amount  of  oxygen  consumption  by  following  equation  using  the 


reaction  2H^  +  =  21^0. 


qR  =  0  -  nFC)  68,204 


cal  ( 

2  gm-mol  H^O 

gm-mol  H?0\ 

32  gm  O 

*  (1.  798823) 


BTU/lb. 


cal/gm 


(V-36) 


or: 

qR  =  7660  (1  - 

The  heat  rejected  from  the  stack  is  assumed  to  be  absorbed  by 
evaporation  of  water  (product  water  plus  additional  water  wicked  from 
the  condenser)  and  by  increasing  the  enthalpy  of  the  circulating  gas, 
or: 


BTU/sec. 


(V-37) 
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Assuming  a  constant  flow  area,  rewriting  Equation  (V-38)  by  using 
the  equation  for  continuity,  Equation  (V-  1)  gives: 


(V-39) 


Multiplying  Equation  (V-39)  by  T^/V^  gives 
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The  terms  in  Equation  (V-40)  containing  the  oxygen  partial  pressure 
can  be  simplified  by  using  Equation  (V-27)  or. 


(V-40A) 


Equations  (V-27).  (V-29),  (V-30)  and  (V-34)  are  then  used  to 
eliminate  all  unknowns  in  Equation  (V-40)  except  (p„  Q)  assuming 

V  1 

that  the  conditions  at  station  2,  the  exit  of  the  stack,  are  known. 
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Equation  (V-41)  is  quite  complex  to  attempt  to  solve  analytically;  therefore, 

the  solution  was  sought  by  making  an  initial  estimate  ^'r 

partial  pressure  of  the  water  vapor  at  the  exit  of  the  fuel  cell  stack  and 

then  determine  a  more  refined  value  (p  JN  ^  by  the  Newton-Raphson 

2  4 

method  which  is  expressed  as  follows- 
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(2)  = 

1 


(V- 42) 


where:  (l),  (2)  indicate  values  based  on 

successive  estimates 

y  =  residual  equation  (residual  is  essentially 
zero  when  a  satisfactory  solution  is 
obtained) 


To  use  Equation  (V-42),  Equation  (V-41)  must  be  different  with  respect 


T 


2 


IP 


1 
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V.  10  Example  of  Solution  to  Heat.  Balance  of  Flow  through  the  Fuel  Cell  Stack 

The  following  example  is  used  *o  illustrate  the  solution  to  the  heat 
balance  equations  for  the  fuel  cell  stack 


Given: 


P 


I 


p^  =  14.  71 1  psia 


7  042 


B 


lb  -mol 
m 


t  =  12o'V  =  579.  7  R 


C  6.9i< 


2 


1 . 6924  psia 


C  =  8.045 

V 


<*  =  ioo  =  o.5 

c  FC 

E  -  0.7  39  volt s 

Dry  air  enters  system 


>V, 


Initial  assumption 


1 


o 

118  F, 


1.  6006 


The  gas  composition  at  station  7  is  the  same  as  at  station  2,  therefore, 
P6-  P2»"“(PH20)7  -  (pH2o)2  Since  the  enenng  air 

/  n  \  =0  The  partial  pressure  of  the  oxygen  leaving  the  fuel 

(pH2o)6 

cell  stack  from  Equation  (V-12}  is  found  as  follows 


lpcO 


=  0.2095(14.711}  -  0.6048(1.6924)  =  2.  058  psia 


2'2 


Equation  (V-40A)  is  used  to  determine  y.  Miscellaneous  terms  wi*hin 


Equation  (V-40A)  are  determined  as  follows 
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_  459.7  +  1  18  _  1 

T^  '  459. 7  +  120  "  1.  00346 


P1  ’  (ph  o)  =  19.711  -  1.6924  13.018b  psia 


P1  *  (PH  o')  =  14.71  1  -  1. 6006  =  1  3.  1  104  psia 


101(13.  1 104,  (2.  0580; 
100(13.  0186)  t  2.  0580 

=  2. 0900  psia 


=  1.01553(2.058) 


(V-29) 


V  /  T 


( l.oi 

VtJ  1.0155 


1 


0155J  1,005475 


(V  -  27) 


"  1 . 005475i 1 . 00418;  1.00968 


1.01 

1 . 00346 


r  1.00651 


-  H  711  -  1.6924  2.0580  =  10.961  psia 


°  00S4T1  0.  00545 

1  V  1.005475 


(ph2°)2‘  (ph2o')1(v2/\t1 


V,VT2)  •  1  6924  -  A  1>00-6-  =  0,1004 

1  -  1. lt 005475 


(PH20)2‘  (PH2o)1l  V 


Ii\  ,  1.6924-  -±4^  =  0.1061 


1. 00895 
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Substituting  in  Equation  (V-40A): 


y  =  18.  016(102  3)(0.  1004) 

+  579.  7  (j.  042(-  0.  OOb51)(2.  058)  +6.  934(10.  961)  •  (0.  00545)  +8.  045(0.  1072j] 

7660(0.  5)(32)(2.  058) 

100 

y  =  1856  +  681  -  2520  =  +17 


Equation  (V-43)  is  now  evaluaied  in  order  to  correct  the  assumed  value  of 

( n  \  =  1. 6006  psia,  or 

' "  H2°)l 


'(%o\ 


7.042(1. 01)(2.  058) 
0.0459 


100(1  3.  0186)  +  2.058 


100(13.  1094) 


71  1(18.  016)11026.  3) 
+  6.  934(579.  7)i  10.  961)  • 

+  8. 045  79  7) 


1  + 


13.0186 


0. 0459(579 


J 


l" 


14.711(1. 6924)  -  2(14.71  1  -  0.  8003)(l .  6006) 
71  1  *  0.  0459(579.  7) 


8y 


#(PH20)l 


-  318.6  - 


272, 000  +  65, 600  +  72, 050 
13.  2 


=  -  318.  6  -  31, 000  =  -  31,  319 
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4.  0  FACTUAL  DATA  (Cont;d) 

From  the  Newton-Raphson  relationship  equation 

(»h,o),<2)  -  ‘  ^-(oTJT 9)  -  1-60,>  ps,“ 

A  repetition  of  the  above  series  of  calculations  does  not  produce  a 
significant  improvement  in  the  value  for  ^p^  q\  . 

From  Keenan  and  Keyes  the  temperature  corresponding  to  ^p^  = 

1.6011  psia  is  tj  -  118.01  deg.  F. 

The  first  term  in  Equation  (V-40A;  is  proportional  to  the  amount  of  heat 
rejected  by  the  evaporation  of  water;  the  term  is  evaluated  as  follows 


'■•(Mi " 


6011  =  13.1  100  psia 


fo  x  101(13.  1 100/(2.  0580 L  =  1.01548(2.0580) 
(p02)i  100(13.  0186)  t  2  058 


( V - 29) 


=  2.  0897  psia 


Ii  (h) .  _i 

v.  It  I  1. 


01 


1 


01548  1.005425 


(qR)  _  <x  18.  016(1026.  3) 
H2° 


1. 6924  - 


1.6011 
1. 005425 


.1  . 


1828 


(V-27) 
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4.0  FACTUAL  DATA  (Cont’d) 


The  total  heat  rejection  is  proportional  to  the  last  term  in 
Equation  (V-40A)  which  was  evaluated  previously  as  (qR)tot  tt  2520. 
The  percent  heat  rejection  by  evaporation  of  water  is  then  found  as 
follows: 


Percent  heat  rejection  by  evaporation 


(qR) 


H2° 


(qR) 


tot 


x  100 


e  -1818(100)  =  n.  ,  percent 


/ 
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CONCLUSIONS 


In  order  to  satisfy  the  multiple  requirements  necessary  to  be 
fulfilled  in  a  broad  approach  at  IEM  fuel  cell  battery  designs, 
systematic  investigations  have  been  undertaken,  emphasizing 
determinations  of  rates  as  related  to  the  electrochemical 
reaction,  the  transport  of  water  in  ion  exchange  membranes, 
detailed  heat  and  mass  transport  across  interfaces  and  overall 
systems  analysis. 

The  program  included  applied  research  as  well  as  development 
activities  in  an  attempt  to  apply  experimental  results  and  con¬ 
cepts  immediately  to  single  cell  investigations.  Appreciable 
improvements  in  air-operated  IEM  fuel  cell  capabilities  have 
results  from  this  type  of  approach. 

Properties  of  the  electrolyte  (IEM) /electrode  interface  have 
become  important  parameters  in  establishing  an  analytical 
model  descriptive  of  interfacial  processes  and  describing 
the  nature  and  quality  of  the  2-phase  bond.  Because  of  systems 
complexity  and  the  limited  amount  of  published  work  in  the 
literature,  novel  measurement  techniques  need  investigating 
and  improving  before  definite  conclusions  will  become  available 

Water  transport  in  ion  exchange  membranes  car.  be  induced  by 
various  gradients,  of  which  pressure  gradient  effects  have  been 
determined.  Classical  equations  for  transport  in  porous  media, 
accounting  for  membrane  properties  have  yielded  a  general 
correlation  to  describe  transfer  rates.  This  result  will  suggest 
structural  membrane  modifications  in  order  to  control  water 
transport  rates  compatible  with  overall  heat  and  mass  transfer 
processes. 

Detailed  rate-controlling  mechanisms  for  the  air-breathing 
IEM/Pt  black  electrodes  still  remain  to  be  elucidated  and  will 
require  more  experimental  work  beyond  the  scope  of  the  present 
program.  However,  the  semi-empirical  relationships  derived 
for  self-breathing  and  forced-flow  conditions  in  air-channels 
have  yielded  parametric  influences  of  gas  properties  (temperature, 
partial  pressures  of  oxygen,  water  and  nitrogen)  operating  con¬ 
ditions  (electrode  surface  temperature,  air  flow  rate)  and  sys¬ 
tems  geometry  (channel  length  and  height).  The  results  suggested 
new  designs  2. 11  owing  for  quasi-uniformity  of  the  main  parameters, 

i.  e.  ,  electrode  temperature,  partial  pressures  of  oxygen  and  water. 

I 

A  broad  approach  to  battery  design  aimed  at  the  multitude  of  air- 
operated  applications  is  now  possible,  as  a  result  of  the  establishment 
of  the  importance  of  major  parameters. 
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A  novel  current  collector  design,  allowing  for  short  air-channels, 
resulting  in  better  surface  distribution  of  the  oxidant  (air)  and 
more  uniform  surface  and  channel  properties  (temperature,  water 
removal,  thermal  conductivity,  etc.  .  .)  has  resulted  in  systems 
stability  and  high  power  capabilities  (50  watts/ft  ).  Availability 
of  experimental  results,  for  steady  state  operation,  will  suggest 
further  modifications  and  applicability  of  cell  voltage  control 
devices. 

Fast  electrode  surface  temperature  and  cell  resistance  response 
to  varying  load  conditions  and  air  flow  rates  has  allowed  the 
design  of  a  control  device  which  can  maintain  cell  voltage  quasi 
time-independent.  Since,  specific  fuel  cell  batteries  may  require 
controlling  devices  for  unattended  operation,  the  designed  system 
is  now  becoming  the  object  for  such  investigations. 

Preliminary  designs  new  available  for  a  system  based  on  the  con¬ 
trol  of  water  partial  pressure,  according  to  parametric  influences 
determined  in  single  channel  studies,  display  over-all  simplicity 
and  expected  system  stability  without  appreciable  etlects  oi 
variations  of  environmental  conditions.  Such  a  system  shows 
promising  applicability  to  multiple  air-operation  applications 

Presently,  sufficient  experimental  data,  information  and  knowledge 
are  available  to  allow  for  confidcr.ce  in  the  successful  design  ot 
air-operated  IEM  Fuel  Cell  Batteries. 
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PROGRAM  FOR  NEXT  INTERVAL 


6.  0 


Planned  work  up  to  the  Final  Technical  Report  (covering  the  period 
from  October  1,  1962  to  September  30,  1963)  will  remain  in  agree¬ 
ment  with  the  Technical  Guidelines,  EPP  No.  58-238  of  July  17,  1962, 
and  include: 


Task  I:  Interfacial  Membrane/Electrode  Properties 

Limited  additional  work  is  planned  for  the  purpose  of 
testing  conclusions  with  hydrophobic  binders  possessing 
critical  surface  tension  values  between  those  of  TFE  and 
Kel-F.  Additional  work  will  be  conducted  on  sample 
preparation  techniques  for  optical  and  electron  micro¬ 
scopy  observations. 

Task  II:  Water  Transport  through  Membranes 

Additional  work  on  hydrodynamic  as  well  as  vapor  perme¬ 
abilities  will  be  conducted  on  other  membrane  systems. 
Work  on  elcctro-osmotic  rates  of  water  transport  will 
be  continued. 

Task  III:  Electrochemical  Oxygen  Reduction 

No  additional  experimental  work  is  planned  Limited 
data  analysis,  as  required  to  interpret  results,  will  be 
conducted. 


Task  IV:  Heat  and  Mass  Transfer  Investigations 

Limited  additional  experimental  work  is  planned.  Emphasis 
will  be  placed  or.  analysis  of  heat  transfer  data  and  inter¬ 
pretation  of  experimental  results  obtained  on  single  cells. 

Task  V:  Design  Parameters  Affecting  Air-Operated  Fuel  Cells 

An  extensive  experimental  program  is  planned,  associated 
with  testing  and  evaluation  of  the  preliminary  design, 
described  in  detail  in  this  Report. 


Took  VI:  Ion  Exchange  Membranes 

Sporadic  supporting  activities  will  be  associated  with 
experimental  work  of  Tasks  I  and  II.  New  membranes 
will  be  investigated  as  electrolyte,  as  suggested  by  Contract 
Officers  of  the  USAELRDL. 
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7.0 


IDENTIFICATION  OF  KEY  PERSONNEL 


During  the  period  covering  January  1,  1963  to  June  30,  1963,  the 
following  principal  investigators  have  been  associated  with  the 
program  tasks: 


Dr.  J.R.  Boyack 
Mr.  D.  W.  Craft 
Mr.  J.  Dankese 
Dr.  S.  J.  Krumbein  - 
Dr.  H.J.  Maget 


Water  Transport  in  Ion  Exchange  Membranes 
Design  Parameters  -  Air-Operated  IEM  Fuel  Cells 
Heat  and  Mass  Transfer  -  Air-Operated  IEM  Ce'l 
Interface  Structure  and  Imerfacial  Properties 
Kinetics  of  the  Electrochemical  Oxygen  Reduction 


Detailed  Time  Distribution  (January  1  to  June  30,  1963) 


A.  Principal  Investigators 


Dr.  J.R.  Boyack,  Physical  Cncmist  &  Electrolytes 
Mr.  D.  Y-  Craft,  Senior  Preliminary  Design  Engineer 
Mr.  J.  Dankese,  Laboratory  Chemical  Engineer 
Dr.  S.J.  Krumbein,  Physical  Chemist-Electrochemistry 
Dr.  H.J.  Maget,  Senior  Chemical  Engineer 


511  hr s. 
543  hrs. 
918  hrs. 
481  hrs. 
73 2  hrs. 
3  185  hrs. 


B.  Supporting  Laboratory  Specialists  and  Technicians 

Mr.  J.  Amore,  Specialist  -  Fuel  Cell  Electrode 
Mr.  W.  A-mstrong,  Specialist  -  Instrumentation 
Mr.  R.  Donovan,  Technician  -  Laboratory 
Mr.  E.  Glosl,  Technician  -  Polye’.ectrolyte  Physics 
Mr.  J.  Hubs,  Fuel  Cell  Evaluation  Lab.  Technician 
Others  with  individual  contribution  less  than  100  hours 


909  hrs. 
132  hrs. 
47  1  hrs. 
287  hrs. 
1  022  hrs. 

473  hrs. 
3  294  hrs. 


C.  Reports  and  Drafting 

Mr.  E.  P.  Breau  and  others  4311 

*  Data  available  up  to  June  16  and  extrapolated  to  June  30,  1963. 

The  technical  background  of  principal  investigators  is  presented  in  the 
following  resumes: 
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7.  0 


IDENTIFICATION  OF  KEY  PERSONNEL  (Cont'd) 

James  R.  Boyack,  Ph.  D.  -  Electrolyte  Physical  Chemist 

Dr.  Boyack  completed  all  his  work  for  a  Ph.  D.  from  the  University 
of  Utah  in  November  1961.  (degree  June  1962)  and  joined  the  Fuel  Cell 
Laboratory  of  General  Electric's  Aircraft  Accessory  Turbine 
Department  in  December  1961.  He  received  his  B.S.  in  Chemistry 
in  1958  at  the  University  of  Utah.  His  doctoral  thesis  was  concerned 
with  zone  diffusion  and  mobility  in  paper  electrophoresis.  Included 
in  his  background  is  a  summer  as  research  assistant  at  Eastman 
Kodak  Company. 

Since  joining  General  Electric.  Dr.  Boyack  has  been  working  with 
electrolyte  transport  mechanisms  in  various  systems  such  as  lon- 
cxchangc  polymers,  gels,  etc. 


D.  W,  Craft  -  Senior  Preliminary  Design  Engineer 

Mr.  Craft  holds  B.S.  M.E.  andM.S.M.E.  degrees  from  Purdue 
University,  1945  and  1949  respectively.  His  principal  fields  are 
thermodynamics  and  fluid  mechanics. 


From  1949  to  1954  he  was  on  the  teaching  and  research  staff  at  Purdue 
University  where  he  instructed  a  senior  aeronautical  engineering 
course  in  gas  turbines  and  jet  propulsion  and  was  engaged  in  gas  turbine 
and  rocket  engine  research. 


He  joined  General  Electric  in  1954  and  was  originally  concerned  with 
the  field  of  turbomachinerv.  He  conceived  and  developed  a  high 
speed  dynamometer  which  uses  a  turbine  as  its  primary  clement  -  now 
used  as  standard  for  small  turbine  development.  He  also  developed 
a  theory  for  the  compressible  fiow  s  through  drag  turbines  and  com¬ 
pressors.  He  designed  and  developed  a  high  suction  specific  speed  fuel 
pump.  Other  important  contributions  were  in  design  of  supersonic 
nozzles,  mass  flow  meters  and  the  preparation  of  digital  computer 
programs  for  thermodynamic  and  fluid  dynamic  analyses. 


His  fuel  cell  experience,  since  1961.  has  been  primarily  in  the  field  of 
mass  transfer.  His  efforts  have  been  instrumental  in  the  development 
of  a  passive  system  utilizing  a  fixed  pressure  gradient  imposed  across 
a  porous  plate,  for  the  delivery  of  water  from  a  condensing  system 
under  zero  gravity  conditions.  Under  his  leadership  this  principle  was 
successfully  applied  to  a  production  version  of  a  hydrogen-oxygen  fuel 
cell  system  in  which  product  water  is  recovered  from  the  fuel  cell  for 
the  potable  water  supply  in  a  manned  spacecraft.  In  addition,  a  patent 
has  been  applied  for  in  his  name  for  a  method  removing  entrained  liquid 

from  a  gas  stream. 
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IDENTIFICATION  OF  KEY  PERSONNEL  (Cont'd) 


Currently,  Mr.  Craft  is  applying  heat  and  mass  transfer  theory  to 
air-breathing  solid  electrolyte  fuel  cells.  He  is  investigating 
methods  of  maintaining  an  optimum  water  balance  in  the  electrolyte 
while  operating  at  high  current  density. 

Mr.  Craft  is  a  member  of  ASME,  Tau  Beta  Pi,  Pi  Tau  Sigma,  and 
Sigma  Xi. 

Joseph  P.  Dankese  -  Laboratory  Chemical  Engineer 

Mr.  Dankese  received  his  B  S.  in  Chemical  Engineering  from  the 
Massachusetts  Institute  of  Technology  in  1954,  and  has  done  graduate 
work  at  Northeastern  University. 

Before  joining  General  Electric  in  early  1961.  he  accumulated  a 
wide  background  in  chemical  engineering  work  at  Microwave 
Associates,  Ionics,  Inc.  ,  Dewey  and  Almy  Chemical  Co.  ,  and 
A.D.  Little,  Inc.  This  included  studies  of  conductive  coatings,  glass 
to  metal  seals,  and  high  temperature  sintering  of  metals;  design  of 
the  electrical  and  hydraulic  systems  for  ion-exchange  membrane 
electrodialysis  equipment  and  field  testing  of  the  final  product,  pro¬ 
cess  design  and  development  for  the  vinyl  chloride  and  butadiene 
type  copolymer  plants;  developed  a  battery  separator  for  electric 
storage  batteries  (U.S.  Patent  No.  3,  026,  366-March  20,  1963), 
analysis  and  process  design  of  the  liquid  oxygen  and  hydrocarbon 
fueling  system  for  the  Tital  missile  -  this  work  involved  heat  transfer, 
mass  transfer  and  thermodynamics;  conducted  experimental  studies 
on  the  effect  of  radiation  on  the  reactivity  of  gaseous  hydrocarb  r  ir 
fuel  cells;  surveyed  various  power  sources  including  batteries,  fuel 
cells,  thermoelectric  and  thermionic  devices. 

In  General  Electric's  Direct  Energy  Conversion  Laboratory,  he  had 
conducted  mass  and  heat  transfer  studies  and  authored  a  G  E 
Technical  Report  No.  62DE-2  on  Mass  Transfer  and  Voltage  Control 
in  Fuel  Cells,  investigated  new  fuels  and  submitted  two  joint  patent 
disclosures  on  borohydride  and  methanol  liquid  fuels,  improved 
electrode  structures,  advanced  ion  exchange  membranes  for  fuel 
cells  (GPS  Membranes),  and  developed  a  modified  air  -  electrode 
current  collector  for  use  with  I.  E.  M.  fuel  cells  for  which  a  paten* 
disclosure  will  be  submitted. 

Mr.  Dankese  is  a  registered  professional  engineer  in  Massachusetts 
and  a  member  of  the  Boston  Chapter  of  the  American  Institute  of 
Chemical  Engineering. 
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IDENTIFICATION  OF  KEY  PERSONNEL  (Cont'd) 


Simeon  J.  Krumbein,  Ph.  D.  -  Physical  Chemist  -  Electrochemistry 

Dr.  Krumbein,  received  his  Ph.  D.  in  Physical  Chemistry  from 
New  York  University  in  October  1961  and  joined  the  Fuel  Cell 
Laboratory  of  General  Electric's  Direct  Energy  Conversion  Operation 
immediately  thereafter.  In  addition  to  being  a  research  assistant 
at  New  York  University  from  1957  to  1959,  he  had  been  working  for 
the  past  two  years  under  a  National  Science  Foundation  Fellowship. 

His  doctoral  thesis  dealt  with  the  effects  of  proteins  on  the  electro¬ 
metric  behavior  of  simple  inorganic  ions. 

Prior  to  this,  he  worked  in  the  laboratory  of  the  Technical  Research 
Group  (TRG),  Inc.  ,  where  he  conducted  analysis  of  the  properties  of 
materials . 

Since  coming  to  the  Direct  Energy  Conversion  Operation,  Dr.  Krumbein 
has  been  primarily  concerned  with  fundamental  electrochemical 
mechanisms  as  they  exist  in  fuel  cell  applications. 

Henri  J.R.  Maget,  D.  Sc.  -  Senior  Chemical  Engineer 

Dr.  Maget  has  a  background  combining  organic  chemistry  and  chemical 
engineering.  He  holds  degrees:  B.S.  and  Licence  es  Sciences  from 
University  of  Strasbourg,  France,  1950  where  he  majored  in  Chemistry; 
M.  S.  in  Petroleum  Engineering  from  Petroleum  School,  Paris,  1953, 
M.S.  in  Chemical  Engineering  from  University  of  Michigan,  1956  ;  and 
D.  Sc.  from  University  of  Toulouse,  France,  1959  He  joined  the 
Aircraft  Accessory  Turbine  Department  in  August  I960. 

His  background  includes  two  years  of  research  at  the  University  of 
Michigan  as  an  exchange  scientist.  This  was  followed  by  two  years 
of  research  work  at  Institut  du  Genie  Chimique  in  France  and  more 
recently  as  research  scientist  at  Yale  University.  His  specialization 
leading  to  his  doctorate  was  chemical  kinetics  and  polymerization  of 
epoxides. 

Since  he  joined  the  General  Electric  Company,  Dr.  Maget  has  been 
working  on  physico-chemical  properties  of  electrochemically  active 
catalysts,  and  is  investigating  the  relationship  between  these  properties 
and  the  performance  of  fuel  cells. 
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